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Abstract 



We propose the dual gravity of a non conformal gauge theory which has logarithmic 
running of couplings in the IR but becomes almost conformal in the far UV. The 
theory has matter in fundamental representation, non-zero temperature and under a 
cascade of Seiberg dualities, can be described in terms of gauge groups of lower and 
lower rank. We outline the procedure of holographic renormalization and propose 
a mechanism to UV complete the gauge theory by modifying the dual geometry 
at large radial distances. As an example, we construct the brane configuration and 
sources required to attach a Klebanov-Witten type geometry at large r to a Klcbanov- 
Strassler type geometry at small r. Using the supergravity description for the dual 
geometry, we compute thermal mass of a fundamental 'quark' in our theory along 
with drag and diffusion coefficients of the gauge theory plasma. We compute the 
stress tensor of the gauge theory and formulate the wake a probe leaves behind as it 
traverses the medium. Transport coefficient shear viscosity rj and its ratio to entropy 
rj/s are calculated and finally we show how confinement of 'quarks' at large separation 
can occur at low temperatures. We classify the most general dual geometry that gives 
rise to linear confinement at low temperatures and show how quarkonium states can 
melt at high temperatures to hberate 'quarks'. 
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Resume 



Nous proposons une theorie de la gravite qui correspond a une theorie de jauge avec 
des constantes de couplage a comportement logarithmique dans I'infrarouge et de- 
venant quasi-conforme dans I'ultraviolet profond. La theorie contient la maticrc dans 
la representation fondamentale, a une temperature non-nuUe et peut-etre decrite en 
terme de groupes de jauge de rangs de plus en plus petits par une cascade de dualites 
de Seiberg. Nous discutons brievement du processus de renormalisation holographique 
et proposons un mecanisme pour completer la theorie de jauge dans I'ultraviolet. 
Comme exemple, nous construisons une configuration de membranes et de sources qui 
joignent une geometric de Klebanov-Witten a grand r a une geometric de Klebanov- 
Strassler a petit r. En utilisant la supergravite pour decrire la geometric duale, nous 
calculous la masse thermique d'un quark fondamental ainsi que les coefficients de 
trainee et de diffusion du plasma de la theorie de jauge. Nous evaluons le tenseur de 
stress de la theorie de jauge et quantifions le sillon laisse par une sonde traversant le 
milieu. Nous calculous le rapport de la viscosite de cisaillement a la densite d'entropie, 
77/ s. Finalement, nous montrons comment le geometric duale la plus generale possible 
peut donner lieu a la fois au confinement lineaire de paires quark- antiquark a basse 
temperature et a la dissolution du quarkonium a haute temperature. 
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Chapter 1 
Introduction 



Since the discovery of a unified description of strong, weak and electromagnetic in- 
teraction back in the 1970's, much activity was centered around its consequences. 
The foundation of the unification hes in the principle of local gauge symmetry where 
interactions between particles are mediated by gauge fields. Through the mechanism 
of spontaneous symmetry breaking as in the Higgs model Weinberg [5] and 

Salam [6J showed that at low energies only photons and neutrinos remain massless 
while vector bosons responsible for weak interaction acquire mass. In the far UV 
the vector bosons become massless and strong and weak interactions become long 
range- just like electromagnetic interaction. Thus at high energies symmetry between 
strong, weak and electromagnetic force gets restored [6] [5]. 

On the other hand phase transition is a process which changes the symmetry of 
the system and before or after the phase transition, symmetry is broken [7]. This 
suggests the existence of various phases in the theory of elementary particles. In par- 
ticular, highly dense matter under extreme temperature and pressure should exhibit 
various phases as temperature and density is altered. By heating nuclear matter up 
to extreme temperature and pressure, one expects to create a new phase of matter, 
a plasma of quarks and gluons where quarks become massless with strong, weak and 
electromagnetic interaction becoming long range, restoring the symmetry between 
them. The estimates for the critical temperature Tc for phase transition in renormal- 
izable gauge theories were first made by [5]-[in] and there is an extensive literature 
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2 Chapter 1 Introduction 

on the phase structure of quark matter along with studies of critical temperature 

m-m- 

These theoretical expectations for a new phase of matter known as the Quark 
Gluon Plasma (QGP) led to the Relativistic Heavy Ion Collider (RHIC) program 
where heavy nuclei are made to collide at relativistic velocities. The matter formed 
in the early stages of the collision has indeed been identified not as a collection of 
color neutral hadrons, but a new state of dense matter [I5]-[T8]. The relativistic fluid 
created at RHIC cannot be described in terms of hadronic degrees of freedom and one 
may conclude that the energies reached by the experiments give rise to a plasma of 
quarks and gluons i.e. QGP. But what are the properties of this fluid? In particular, 
is the coupling between the constituents of the fluid strong or weak? The answer to 
this is crucial as it will determine the theoretical tools to study the fluid. At weak 
coupling, perturbative methods can be useful whereas at strong coupling, effective or 
lattice field theory need to be applied. 

To address this issue, first observe that one of the key characteristics of a fluid is 
how it responds to pressure variation and flows to equilibrium state. At the heavy 
ion collider, the overlapping region between two heavy nuclei is of elliptic shape with 
one axis longer than the other- thus creating anisotropic pressure. The fluid will of 
course try to equilibrate to spherical symmetric configuration and this leads to elliptic 
flow. This flow as observed in the experiments at RHIC [IS]- [22] is well described 
by ideal hydrodynamics which assumes that the QGP fluid is strongly coupled. This 
means perturbative field theory techniques are not useful to analyze this strongly 
coupled system and one usually studies the theory on the lattice or using effective 
field theory . While the former becomes quite challenging with lattice simulations 
limited by computational ability of the numerical method, the latter relies on effective 
Lagrangians which are only approximations and could lead to incomplete results. 

But these are not the only tools to describe a strongly coupled field theory. One 
may apply the principle of holography which relates the Hilbert space of a gauge 
theory with that of a theory of gravity, to study strongly coupled field theories. The 
key observation was made by Maldacena in the late 1990 's while studying anti de 
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Sitter black hole solutions and how they arise from D brane configurations [23]. He 
conjectured that strongly coupled A/" = 4 supersymmetric conformal field theory is 
dual to a theory of weakly coupled gravitons describing AdS^ x 5*5 geometry, where 
AdS^ is the five dimensional anti de Sitter space and 5*5 is the five sphere. This 
duality known as the AdS / CFT correspondence opened a whole new avenue to study 
strongly coupled quantum field theories using weakly coupled classical gravity and is 
only a part of a more general correspondence between gauge theory and geometry. 

In this thesis we extend the AdS/CFT correspondence to incorporate renormal- 
ization group fiow in the gauge theory and study the system at strong coupling using 
weakly coupled dual gravity. We propose the dual geometry for a strongly coupled 
gauge theory which has logarithmic running of couplings in the far IR but becomes 
almost conformal in the far UV. The gauge theory has matter in fundamental repre- 
sentation and has non-zero temperature- which is achieved by introducing a black hole 
in the dual geometry. Using the weakly coupled gravity description for the strongly 
coupled field theory, one can easily compute expectation values of various gauge the- 
ory operators that are extremely difficult to obtain using conventional techniques. By 
computing quantum correlation functions using dual classical action, one can analyze 
the thermodynamic properties and learn about the kinematics of strongly coupled 
gauge theory plasma. But how does the rich structure of quantum field theory get 
exact description in terms of a simple classical theory of gravity? To be more precise, 
what are the limits of validity of the correspondence? 

It turns out that for a field theory at strong coupling to have a description in 
terms of weakly coupled gravity, the number of colors N in the gauge theory must 
be very large. In the large limit, it seems that the field theory has a classical de- 
scription in terms of dual gravity where quantum correlation functions can be exactly 
calculated. At first glance this simplification might appear to be quite 'accidental' 
but a careful analysis can shed light on the issue. The answer may be linked to an 
observation made by 't Hooft who suggested that the number of colors A^ in a gauge 
theory can be thought of as an expansion parameter [2l]. There are planar and non 
planar diagrams that contribute to the propagator and in the limit A^ — )■ oo, only 
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planar diagrams survive. This is because non planar diagrams are suppressed by 
0{1/N^),k > 1, relative to planar ones [2l][25], resulting in a drastic simplification 
of the propagator. While at finite propagators and subsequently S matrix elements 
are extremely difficult to compute due to the large number of non planar diagrams at 
various loop order, in the large N limit, the theory has a rather simple description. 
This simplification from a field theory analysis is completely consistent with our dual 
gravity description where theory can be described by classical gravity. 

The thesis is organized as follows: In chapter two, after a brief discussion on 
AdS / CFT correspondence, we describe in some detail gauge theories that arise from 
various brane configurations and then present their dual geometries. We outline the 
procedure of holographic renormalization of non conformal gauge theories which have 
dual gravity description and then discuss how to UV complete the theory by attaching 
geometries. In the final section of chapter two, we give a specific example of such 
UV completion by considering localized sources which allow one to attach asymptotic 
AdS geometry at large r, to a Klebanov-Strassler type geometry at small r. Using 
the dual geometries of chapter two, in chapter three we compute various gauge theory 
quantities crucial in analyzing a plasma. We compute the thermal mass of 'quark' in 
fundamental representation, the drag it experiences as it moves through the plasma 
and the wake it leaves behind as it traverses the medium. We also compute the 
momentum broadening of a fast moving jet along with transport coefficients such as 
shear viscosity and its ratio to entropy. Finally we show how confinement can be 
achieved using dual gravity and propose the most general dual geometry that can 
realize linear confinement at large distances. We conclude the thesis with a brief 
summary of our results along with future directions that can be explored. 



Chapter 2 

The Gauge/ Gravity 
Correspondence 

Quantum field theory with gauge symmetry has been extremely successful in describ- 
ing the dynamics and collective excitations of highly energetic particles at very short 
distance scales. On the other hand, Einstein's theory of gravity has been crucial in 
understanding the dynamics of massive objects separated by large distances under 
the force of gravity. But what is the connection between these two seemingly distinct 
theories at opposite distance scales? General relativity dictates that any stress energy 
will couple to gravitons and thus given a field theory one can always compute the 
gravitons sourced by the stress tensor of the field theory. Thus given a gauge theory, 
we can always compute the corresponding gravitons and there is a natural equivalence 
between energies of the gravitons with that of the fields. Note the coupling of field 
theory with gravity is very weak compared to other gauge couplings (such as strong 
interaction coupling ccg or the electromagnetic or weak interaction coupling) up to 
the energy scale relevant for current experiments. Nevertheless there is a naturalness 
in this correspondence between gauge theory and the geometry it sources, all because 
of the coupling of the two. 

In string theory, a more remarkable and a little less obvious duality arises. Open 
string excitations are described by a quantum field theory which has scalar and vector 
fields. In addition, closed string excitations are described by quantum fields which 
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6 Chapter 2 The Gauge/Gravity Correspondence 

are tensors of rank two and higher [26] [27] [28] . Thus open string quantization can 
give rise to gauge theories with fermions and spin one bosons whereas closed string 
quantization can incorporate gravitons. In general there will be interactions between 
open and closed string modes but it turns out that when the interaction is weak, 
there is indication that the Hilbert space of the open string modes and the Hilbert 
space of the closed string modes are identical. Thus the conjecture [23] that open 
string field theory is dual to closed string field theory in the limit where the modes 
decouple! This remarkable duality is one the primary motivations for proposing a 
general correspondence between a gauge theory and theory of gravity - in the limit 
they decouple from each other. In the following sections we will explore this gauge 
gravity duality, first for conformal field theories and then for gauge theories with 
running couplings. 

2.1 Conformal Field Theory and AdS Geometry 

Gauge theories naturally arise when one studies the excitations of branes. In particu- 
lar excitations of open strings ending on D branes can be described by supersymmetric 
Yang- Mills multiplet with vector, spinor and scalar fields [22]. If we consider strings 
ending on a single D3 brane where the D3 brane is embedded in ten dimensional 
flat space time, then the gauge group is U{1). The D3 brane fllls up four dimen- 
sional Minkowski space and can move in six dimensional space with coordinates 
x\j = 4, ..,9. its world volume has coordinates = 0,1,2,3 and the massless 
gauge bosons Ai{x'^),i,l = 0, ..,3 live on the D3 brane world volume. There are six 
scalar flelds (f)j{x'^),j = 4, .., 9 and they describe the oscillation of the position of the 
D3 brane. 

On the other hand if we consider N number of parallel D3 branes, then two ends 
of a string can lie on different branes or on the same brane. There are A^^ choices to 
place the endpoints which gives N"^ number of vectors. Thus the excitations of the 
strings stretching between and ending on the branes is described by a U{N) gauge 
theory [2^] [2^] with vector, scalars and spinor fields. Note that one takes the massless 
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modes of the string excitations which means strings have vanishing length and imply 
that we have coincident branes. Thus excitations of coincident D3 branes give rise 
to U{N) gauge theory living in four dimensional flat space time. It also turns out 
that the gauge coupling does not run and we have a conformal field theory with U{N) 
gauge group in four dimensions. 

The total action for the collection of branes embedded in the ten dimensional 
space time can be written as [30] 



^ i^branes ~l~ i^gravity ~l~ ^int (^■-^) 

where iSbranes is the action for theory on the branes, iSgravity is the action for the back- 
ground geometry on which the branes have been embedded and iSmt is the interaction 
of the brane theory with gravity i.e the the interactions of the gauge fields with gravi- 
tons. At low energy, the coupling of the gauge theory with the gravitons is negligible 
and we can ignore the interaction term to obtain a gauge theory living in fiat space 
and a theory of free low energy gravitons. 

Let us analyze the system of branes embedded in this geometry using supergravity 
which is the low energy limit of string theory. The supergravity solution incorporating 
the back reaction of the flux sourced by the branes was computed back in the 90's 
1. The metric reads [29] 



ds^ = -^(-Gir)dt^ + dl^')+Z^dr^ + r^VHdnl 
VH^ ^ G{r) ^ 



~4 



H = l + L7r^ ^^ = 1--^, L^=47rgsNa 



(2.2) 



where dQf is the metric of five dimensional sphere S^, Qs is the string coupling and 
a' = 11 is the string scale. Here horizon is located at r = fo and extremal limit is 
achieved by taking fo = 0. In the extremal limit, we see that we have an horizon at 
r = and for the near horizon limit, r <^ L, the above metric takes the following 
form 

ds^ = ^^[-de + dl^^)+^dr^ + LHnl (2.3) 
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which is the metric of AdS^ x 5*5, where the AdS throat radius given by L. For an 
observer located at the boundary r = oo, gravitons from the near horizon region 
appear to be of low energy as they are red-shifted. As the energy of the gravitons 
should be measured by an observer at the boundary, low energy modes consist of 
the gravitons coming from the near horizon region i.e. AdS^ x 5*5 geometry and low 
energy gravitons from the bulk. 

Thus we have two descriptions of the same system of brane configuration at low 
energy; one in terms of branes and low energy gravitons, the other in terms of AdS^ x 
5*5 geometry and low energy gravitons. We can identify the two descriptions and 
conclude that the brane theory described by iSbranes has an equivalent description in 
terms of gravitons of AdS^ x geometry. Thus the gauge theory on the branes 
which is a Conformal field theory in flat four dimensional space time is dual to Anti 
de-Sitter geometry pS]. This duality is known as the AdS/CFT correspondence. 



Figure 2.1: Open string gauge theory from D3 branes and its dual geometry. 





D3 Branes 



AdS geometry 



Gauge Theory 



Dual Gravity 



One may ask at what regime of couplings constant space is this duality valid? To 
answer this, first note that the string coupling Qs which describes the coupling between 
the strings and thus the coupling between the gauge fields, can be identified with the 
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Yang-Mills coupling gyu = 4:T^gs- On the other hand the t' Hooft coupling for U{N) 
gauge theory is A = Qym^ ^^us the AdS throat radius L"^ = g^y^Na'"^ = Xa'^ in 
a' units is the 't Hooft coupling. The supergravity solution in fl2.2f2.3p is only valid 
when curvature of space is small. Observing that the Ricci scalar for AdS^ and 5*5 
are of 0{L~'^)., we conclude that for the solution in f l2.2|2.3p to be valid, we need L 
to be very large. Keeping in mind that length is measured in string units a', this 
means gsN = X/Att ^ 1. We also want to keep Qs small so we can ignore effects 
of string loops in the gravity action. The only way to keep (7^ <^ 1 and still have 
QsN ^ 1 is to take — 00. Thus in these limits, we have large strongly coupled 
gauge theory with 't Hooft coupling A ^ 1 dual to Anti de-Sitter geometry with small 
curvature. Hence strongly coupled Conformal field theory is dual to weakly coupled 
Anti de-Sitter gravity. 

2.2 Non Conformal Field Theory and Dual Geom- 
etry 

The duality between conformal field theory and AdS geometry is a particular example 
of the Holographic Principle mapping Hilbert spaces of gauge theories with that of 
gravity. By studying brane excitations in various background geometries one can 
attempt to generalize the AdS / CFT correspondence to include theories with running 
couplings. In the following sections we first analyze the gauge theories that arise from 
excitations of branes placed in geometries with conical singularity and then describe 
their weakly coupled dual gravity. 

2.2.1 Gauge Theory from Brane Configuration 

We will consider branes placed in geometries with conifolds and study their excita- 
tions. Before going into the brane setup, we briefly discuss the conifold geometry. 
For details, consult [32]-[38] 

Consider a six dimensional cone with base T^'^ and radial coordinate r. This cone 
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is a manifold with conic singularity i.e a conifold which is a solution to the Einstein 
equations in vacuum and it has the metric 



dsl = dr"^ + r'^ds%i_,i 



(2.4) 



The metric of the base T ' is given by 



ds. 



dip + cos9id(j)i 



i=l 



1 ^ 



(2.5) 



where ip is an angular coordinate ranging from to Att and {6i, 0i) , {62, ^2) param- 
eterizes the two two-spheres S'^'s. The form of the metric makes it clear that T^'^ is 
a f/(l) bundle over S*^ x S"^ and has the topology of S'^ x S*^. 

Recall that a two dimensional cone embedded in three dimensional space has the 
familiar embedding equation 



2,2 2 
X + y = z 



where the base of the cone is a circle with radius z"^ = r^. The six dimensional 
conifold is a generalization to higher dimensions and has the embedding 



Z1Z2 — ^3^4 — 



in C^, where Zi are complex coordinates given by 



(2.6) 



Z2 
Z3 
Z4 



^3/2gi/2(V'-0l-02) 
^3/2gj/2(V'+0i+9i2) 



sin {9i/2) sin {62/2) 
cos (^1/2) cos (^2/2) 



^3/2gi/2(^+<^i-<^2)^Qg (0i/2) sin (62/2) 



^3/2gi/2(V'-</'l+02) 



sin {9i/2) cos {62/2) 



(2.7) 



To see the symmetries of the space, we can parameterize the conifold with another 
set of complex coordinates Wi given by 



Zl = Wi+ iW2, Z2 = Wi - iW2 

Z3 = —W3 + iw4, Z4 = —W3 — iwi 



(2.8) 



2.2 Non Conformal Field Theory and Dual Geometry 
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In terms of Wi coordinates, the conifold equation (12. 6 p becomes 

E^' = (2.9) 
On the other hand, the base T^'^ is the intersection of the cone with surface 

Ekil' = ^' (2.10) 

We note from (12. 9p . (12.1 OP that T^'^ is invariant under rotation of the four Wi coordi- 
nates, that is under the group 5*0(4) ~ SU{2) x SU{2) and an overall phase rotation. 
Thus the symmetry group of T^'^ is SU{2) x SU{2) x U{1). This will come in handy 
shortly. 

With the understanding of the symmetries of our conifold geometry, consider 
embedding D3 branes in ten dimensional manifold with the metric 

dslo = -dt^ + d-^^ + dsl (2.11) 

where ds\ is given by (12. 4p . That is we have four dimensional Minkowski space along 
with the six dimensional conifold. The D3 branes live in the flat four dimensional 
space and and are placed at the tip of the conifold at fixed radial location r = 0, as 
shown in Fig. 2.2. 




D3 Branes 



Figure 2.2: D3 branes placed on conifold singularity. 

The excitations of the massless open strings ending on these D3 branes are de- 
scribed by gauge fields and complex matter fields Ai,Bi, i = 1,2 which transform as 
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bi-fundamental fields under the gauge group SU{N) x SU{N). Note that the matter 
fields Ai,A2 transform under global SU(2) and so do -81,-82 under another SU{2) 
and we also have global U{1) phase rotation. Thus we have SU{2) x SU{2) x U{1) 
global symmetry, which is also the symmetry of the conifold! This is not surprising as 
these fields describe motion of the D3 branes [26] and the branes move in the conifold 
direction. Thus the fields Ai, Bi are really coordinates-ordinates of the conifold and 
can be written as 



zi = AiBi, Z2 = A2B2 

= A1B2, Zi = A2B1 



(2.12) 



Having analyzed the global symmetry, we would like to understand the origin of 
the gauge group SU (N) x SU (N) and the bi-fundamental nature of the matter fields. 
The nature of the gauge group becomes clear once we analyze the T dual setup of 
the Type IIB brane configuration in conifold geometry. We will be brief on what 
follows and for detailed discussion on T dual of Type IIB which is Type IIA brane 
constructions, please consult [39] -j 



strings giving 
matter in (N,N) of 
SU(N) X SU(N) 




Figure 2.3: T dual of KW model. 



Under T duality the cone becomes two intersecting NS5 branes extending along 
the 012345 and 012389 directions, while the D3 branes become D4 branes along 01236 
directions [10]. Here 0123 are the four Minkowski directions , 6 is the radial direction 
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and 45789 are the angular directions. The T dual setup is shown in Fig 2.3 where 
we have suppressed 0123 and the 7 direction. As D4 branes can end on NS5 branes, 
they get divided into two branches between the NS5 branes along the 6 direction as 
shows in Fig 2.3. On each branch there are N D4 branes giving rise to SU{N) gauge 
group. But there are also strings with one end on a D4 brane on one side of the NS5 
and the other on the other side NS5 brane, as in Fig 2.3. These strings give rise 
to matter fields which transform as (A^, N) of SU (N) x SU {N) gauge group. Thus 
D3 branes at the tip of the conifold gives SU{N) x SU (N) gauge theory and this is 
known as the Klebanov-Witten (KW) model [55] . 

With the understanding of the bi-fundamental nature of the gauge group, we now 
analyze the various couplings of the KW model. There are three couplings namely 
the gauge couplings gi, g2 corresponding to the gauge group SU (N) x SU{N) and the 
coupling h to quartic super potential [551 US]- The beta functions for the couplings 
which govern how they change with energy scale /i is given by 



where r] = hfi is the dimensionless coupling of the theory. Observe that all the fields 
in the theory have the same anomalous dimension 70 ((^i, 92, h) [33| |46] and the three 
beta functions vanish exactly when 



which is one equation for the three couplings. The couplings fiow with scale /i and 
once they reach a value such that {gi, g^, h) satisfy equation fl2.14p . there is no fiow in 
the theory. The theory reaches a fixed point and all the couplings stay at that value as 
energy scale is changed thereafter i.e. the field theory becomes conformal. As equation 
fl2.14p is an equation of a surface in three dimensions, the fixed points in this theory 
form a two-dimensional surface in the three-dimensional space of couplings (see Fig 
2.4 for details). The fiow of the couplings with scale is of course the Renormalization 
Group (RG) fiow as illustrated in Fig 2.5 below. Since the sign of the two beta 




'io{gi,g2,h) 



1 



(2.14) 
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Figure 2.4: The two-dimensional RG surface in the Klebanov-Witten theory. 

functions are negative any arbitrary flow in the coupling constant space brings us to 
the fixed point surface. Thus the brane setup in Fig. 2.2 leads to a conformal field 
theory with gauge group SU{N) x SU{N). 

Now consider embedding M D5 branes in the conifold geometry with D3 branes 
where the D5 branes wrap the two cycle 5^ of the conifold base and extend in four 
Minkowski directions. The D5 branes will slide down to the tip at r=0 and we have 
the brane setup of Fig. 2. 6. We again have bi-fundamental matter fields Ai, Bi with 
global symmetry group SU{2) x SU{2) x U{1). But now due to the additional D5 
branes wrapping one of the S^'s, the gauge group is SU{N + M) x SU{N). 

The bi-fundamental nature of the fields and the different size the gauge groups 
become clear when we go to the T dual picture as shown in fig. 2.7. The M D5 
branes which wrap only one of the S'^'s (albeit of vanishing radius) become D4 branes 
which stretch between the NS5 branes only on one branch and thus contribute to 
only one of the gauge groups [H]. This effectively gives N + M D4 branes stacked in 
one branch between the NS5 branes where the other branch has D4 branes. This 
gives rise to the group SU{N + M) x SU{N). 

As M 7^ we can define kM = N + M and then the gauge theory is SU{kM) x 
SU{{k — 1)M). The gauge couplings are now gk,gk-i for the two gauge groups 
respectively and we also have rj to be the dimensionless coupling of the quartic super 
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Two dimensional RG 
fixed points surface 



Figure 2.5: The typical RG flows in the Klebanov-Witten theory. 



potential. The three beta functions now are ^51 H6] : 



IGtt^ 



;i + 27o) + f(l-7o) 



1 - 



kM 



/3, = r7(l + 27o) 



gU{k-l)M 
167r2 



:i + 27o)-^(l-7o) 



1 - 



(2.15) 



from which we see that they differ from (12.131) by 0{l/k) factors. Note that there 
is no point in the coupling constant space where all the three beta functions vanish 
exactly and thus the theory has no conformal fixed points (except when all couplings 
are zero and we have free theory at all scales). 

Solutions to (I2.15P along with boundary conditions determine the RG flow of 
theory and the flow in the space of coupling constants is depicted by the arrows in 
Fig 2.8. From (12.151) we see that the gauge couplings run logarithmically with scale 
fi and when both gk, Qk-i are not zero, the difference between them grows with scale. 
This is the Klebanov-Strassler (KS) model |15]. As the QCD coupling as also runs 
logarithmically with scale, the gauge theory we obtained above could have common 
features with QCD. Furthermore, as we shall see in the next section, this gauge 
theory has a dual gravity description and we may learn about strongly coupled QCD 
by analyzing the dual geometry of the brane setup of Fig 2.6. But before making 
connections to QCD, lets flrst analyze the RG flow of KS model. 



16 



Chapter 2 The Gauge/ Gravity Correspondence 




D5 Branes 



D3 Branes 



Figure 2.6: D3-D5 system placed at conifold singularity. 



The RG flow incorporates a Seiberg duality cascade, where under a series of du- 
alities higher rank gauge groups have equivalent description in terms of lower rank 
groups. Here we will briefly review the Seiberg duality cascade and its realization in 
the KS model. For a comprehensive review please consult HT]. 

Seiberg duality states that strongly coupled SU{M) gauge theory with Nf flavors 
is dual to weakly coupled SU {Nf — Af) gauge theory with Nj flavors. Recall that 
we have SU{N + M) x SU{N) gauge group and the SU{N + M) branch has 2N 
effective flavors while the SU{N) branch has 2(A^ + M) flavors [HI HB]. To see this, 
consider the fields Ai,A2 and note that each of them have color indices under two 
color groups SU{N + M) and SU{N). If one fixes the color under one of the group, 
say SU{N + M), then the color indices for the group SU{N) can be thought of as 
flavor indices. Then for the SU{N + M) color symmetry group, the SU (N) group 
appears as flavor symmetry group. As we have two fields Ai and A2, this means that 
the SU{N + M) branch has 2N effective flavors. The same argument shows that 
SU{N) group has 2{N + M) effective flavors. 

Now consider the flow on the rj = plane given by curve 2 of Fig 2.8. The 
coupling gk-i corresponding to gauge group SU{N) shrinks while gk corresponding 
to SU{N + M) grows as scale is changed from UV to JR. . At the end point of curve 
2, the SU{N + M) gauge group is strongly coupled with 2N effective flavors and thus 
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strings giving 
SU(N) 



string giving 
SU(N+M) 



strings giving matter 
in (N, N+M) of 
SU(N)X SU(N+M) 




Figure 2.7: T dual of KS model. 



it is dual to weakly coupled SU{2N - {N + M)) = SU{N - M) gauge theory under 
Seiberg duality. This means the end point of curve 2 has two equivalent descriptions; 
one in terms of SU {N + M) strongly coupled gauge theory with 2N flavors and the 
other in terms of SU{N — M) weakly coupled gauge theory with 2N flavors. Thus 
SU{N + M) X SU{N) gauge theory is dual to SU{N - M) x SU{N) gauge theory. 
But N — M = {k — 2)M and thus we can draw the endpoint of curve 2 on another 
coupling constant space with couplings gk-i, gk-2 as shown in Fig 2.9. In fact the RG 
flow of curve 3 can be identified with that of curve 1 in the space of new couplings 
gk-i,gk-2 and new quartic coupling fj. The RG flow of the new couplings gk-i, gk-2, V 
can be obtained by replacing k with k — 1 and rj hj fj in fl2.15p . Solving the RG 
equations, one observe that when gk of SU{N + M) sector grows, gk-2 of SU{N — M) 
shrinks, so the flows of curve 3 in Fig 2.9(a) and of curve 1 in Fig 2.9(b) are in 
opposite direction. This also means 77 ~ 1/f/ and the identiflcation of curve 3 with 
curve 1 is shown in Fig 2.9 is justifled. 

Denoting N — M = N, we now have SU {N + M) x SU (N) and repeating the 
same arguments as before we can perform another Seiberg duality to obtain the 
SU{N — M) X SU{N) gauge group. At each step of the duality, the effective N 
(deflned from gauge group SU{N + M) x SU{N) is reduced by M units and we have 
a cascade of dualities known as the Seiberg duality cascade. The reduction continues 
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Figure 2.8: RG flow of KS model. 



until we reach a point when one of the groups has zero size and we end up with 
SU{0) X SU{M) = SU{M) gauge group. 

By suppressing the quartic couphng ri,fj, we can draw the RG flows of the couphngs 
{gk, gk-i, Qk-n) on the same plot as depicted on Fig 2.10. Each hne of the boundary 
is the axis for couphng gj^ and its Seiberg dual 5'^_2- This depiction is indeed valid 
as when g^ grows and becomes strong, its dual gj^_2 shrinks and becomes weak and 
vice-versa. Intersecting lines describe gauge couplings {gyg%_]) and we move from 
one pair of couplings to the other with the flow. The cut in the Fig 2.10 indicates 
that the flow do not take us back to the original coupling gk-i-, rather the flow takes 
us to another coupling gk-n and we move from one sheet of coupling constant space 
to another. 

Thus we see that in the KS model there are multiple equivalent descriptions of 
the same gauge theory and each description is related to the other by Seiberg duality. 
As we go from UV to IR, we can describe the gauge theory with lower and lower rank 
groups. In the language of Wilsonian RG flow, we are integrating out UV modes 
and obtaining effective Lagrangians for lower and lower rank gauge groups in the 
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IR. If we keep all the relevant, irrelevant and marginal operators as we integrate out, 
physical observables of course should not change with the RG flow. But if we leave out 
certain operators to obtain the effective Lagrangians, the flow will result in changes 
in physical observables. We will have more to say about this subtlety in the coming 
sections. 

In summary, the gauge theory in the KS model has a very rich structure with the 
rank of the gauge group growing in the UV. In fact the far UV of gauge theory is 
Seiberg dual to an infinite rank gauge group while the far IR is a rather simple SU{M) 
gauge theory. But we still do not have matter in the fundamental representation and 
to make connection with QCD, we need to amend the KS brane setup. 

In order to obtain fundamental matter, one introduces D7 branes. The D7 branes 
can be embedded in various ways [381 HHl HSl EOl EH [521 153] and in particular the 
D7 brane world volume will depend on the background geometry along with the 
boundary conditions. In Ouyang's model [3H] seven branes are embedded via the 
following equation (see also [5^): 



3 

z = raexp 



- 01 - 02) 



6l 02 

sin — sin 1^ — 1^ (2-16) 
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Figure 2.10: RG flow of tlie gauge tlieory as it cascades under Seiberg duality. 

where /i is a complex quantity. In the limit where — )■ 0, the seven branes are 
oriented along two branches: 

Branch 1 : = 0, 0i = 

Branch 2 : 02 = 0, 02 = (2.17) 

From the above observe that the seven branes in branch 1 wrap a four cycle {62, 02 ) 
and {ip, r) in the internal space and is stretched along the space time directions 
{t, X, y, z). Similarly in seven branes on branch 2 would wrap a four-cycle (^1, 0i, r, -0). 

With the above embedding one needs to check whether Gauss's law is violated. 
As the seven branes wrap a non-compact four cycle filling the entire r direction, the 
field lines sourced by the branes extend only in the compact directions. If we draw 
a Gaussian surface, all the field lines that go out of the surface come back into the 
surface as the space is compact. This means there cannot be any net charge due 
to the seven branes. This paradox can be resolved by allowing the seven brane to 
wrap a topologically trivial cycle so that it can end abruptly at some r = rmin when 
the embedding is (12.161) . Thus the D7 brane extend in the radial direction with 
'^min < r < 00 which is similar to the seven brane configuration of [IHl ESj • Of course 
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there are Nf D7 branes which can overlap and the final configuration of branes is 
sketched in Fig 2.11. In section 2.4 we will discuss modification of the embedding 
f l2.16p which will be more relevant in describing large QCD. 




D5 Branes 
D3 Branes 



Figure 2.11: D3-D5 branes along with D7 branes embedded in conifold geometry. 



In addition to the bi-fundamental fields Ai,Bi, introduction of the D7 branes give 
rise to fiavor symmetry group SU{Nf) x SU{Nf) and matter fields q,q,Q,Q which 
transform as fundamental under the gauge group SU{N + M) x SU{N) [38]. The bi- 
fundamental nature of the fields follows from the same T duality argument as before 
and the beta functions now take the form: 



glkM 



(l + 27o) + f(l-7o) , Njil 



167r2 



1 - 
DM 



9lt 



+ 



27o 



87r2 



kM 



167r2 



:i + 27o) - - 7o) ^Nf{l- 27,) 



gl_^{k-l)M 
87r2 



{k-l)M 



/3, = r7(l + 27o) 

(2.18) 



where 7, is the anomalous dimension of the q field. Again we have logarithmic running 
of the couplings but also matter which transforms as fundamental of the gauge group. 
In Table 2.1, we list the various matter fields and their representation under local 
and global symmetry groups for the brane setup in Fig 2.11. 
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Field 


SU{N+M)xSU{N) 


SU{Nf) X SU{Nf) 


SU{2) X SU{2) 


Q 


(N + M, 1) 


(Nf, 1) 


(1,1) 


Q 


(N + M,l) 


(l,Nf) 


(1,1) 


Q 


(1,N + M) 


(N7,l) 


(1,1) 


Q 


(1,N + M) 


(1,N7) 


(1,1) 




(N + M,N + M ) 


(N;,Nf) 


(2,1) 


-Bl,2 


(N + M,N + M ) 


(Nf,N7) 


(1,2) 



Table 2.1: The field content and their representation under symmetry groups. 

The field theory here also cascades to lower and lower rank gauge groups under 
Seiberg duality. But now due to the presence of fundamental fiavor, the effective 
flavor for SU {N + M) strongly coupled gauge group is 2N + Nf and its Seiberg dual 
weakly coupled theory is SU{2N+Nf-N -M) ^ SU{N-M + Nf). This means that 
at the ith step of the cascade, the effective N (again defined from the gauge group 
SU{N + M) X SU(N)) is reduced by (M-kNf) units, with k a natural number. This 
reduction continues until we end up with just one group SU (0) x SU{M — jNf) = 
SU{M ~ jNf). Here j is the number of dualities performed starting with gauge group 
SU{N + M) X SU{N) at some UV scale. This brane setup in Fig 2.11 is the Ouyang 
embedding of D7 branes on Klebanov-Strassler model and we will refer to it as the 
Ouyang-Klebanov-Strassler (OKS) model. 

In summary, for OKS model we have a field theory which is non conformal with 
matter in fundamental representation with Nf fiavors and the gauge couplings run 
logarithmically with scale. At some UV scale the gauge group has description in terms 
of SU{N + M) X SU{N) group while at the far IR we end up with SU{M - jNf) 
gauge group. 

For the gauge theories arising from branc configurations in Fig 2.2, 2.6 and 2.11, 
we would like to compute matrix elements as they are physical quantities of interest. 
Note that it is the t'Hooft couphngs = (A^ + M)gl{K) and Afe_i = Ngl_^{K) that 
are relevant for computing propagators and scattering amplitudes. At a scale when 
either of the two t'Hooft couplings A^, \k-i become large, perturbative methods fail 
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and we need to resort to other techniques. One approach is to study the theory on 
the lattice while the other is to find the dual gravity. We will discuss the latter in the 
following section. 

2.2.2 The Dual Gravity of the Brane Theory 

Having discussed in some detail the properties of gauge theories that arise from branes 
in conifold geometry, we will now analyze the same brane configurations at low energy 
using supergravity (SUGRA) approximation of string theory. Type IIB supergravity 
action including local sources in ten dimensions is [55] [56] : 

+ I C4 A A + -V / ^A^i^ + (2.19) 
Jsg SzKfo J Imr 

in Einstein frame, where r = Cq + ie'"^ is the axio-dilaton with Co being the axion 
and the dilaton field and is the five-form fiux sourced by the D3 branes. Here 
G3 = F3 — tH^ with F3 the RR three form fiux sourced by D5 branes and = dB2 
the NS-NS three form fiux with B2 being the NS-NS two form. We also have C4 the 
four-form potential, G = ^Jdet~G^ with G^,y being the metric, -R(2) the curvature 
two-form, and S'loc is the action for localized sources in the system (i.e. D7 branes 
and other local sources that we may consider). The above action fl2.19p is the most 
general supergravity action one gets by placing branes in various geometries and we 
will consider the relevant terms for the the brane setups of Fig 2.2, 2.6 and 2.11 
separately. 

First consider the brane setup of Fig 2.2. The D3 branes sources Ramond- 
Ramond (RR) five form fiux in the supergravity action but do not source three form 
fiux and neither the axio-dilaton field. Minimizing the action f l2.19p with the only 
F5 7^ being the source, we get the following metric [55] 




'1,1 
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H = 1 + h, h = —, L^ = AnQsNa'^ 
QsF^ = d^x A dh-^ + * (rf^x A dh'^) (2.20) 

where (is^i,i is the metric of T^'^ given by (12. 5p . * is the hodge star operator and we 
have taken the extremal hmit.The horizon is located at r = and the near horizon 
limit (r ^ L) of the above metric (I2.20p is that of AdS^ x T^'^ . Using the same 
arguments as in section 2.1, the near horizon gravitons appear to be of low energy as 
measured by an observer at the boundary r = oo. On the other hand at low energy, 
the SU{N) X SU {N) gauge theory that lives on D3 brane world volume i.e on four 
dimensional space, decouples from gravity. Thus we conclude that SU{N) x SU{N) 
gauge theory in four dimensional fiat space is dual to AdS^ x T^'^ geometry. 

But in what regime is this duality valid? For our supergravity solution to hold, 
we need L ^ 1 in a' units, which means that QsN ^ 1. As we want to ignore 
effects of string loops in supergravity action, we need Qs 1 and combined with the 
requirement QgN ^ 1, this means — ?■ oo. On the other hand the gauge couplings 
91,92 corresponding to the gauge group SU{N) x SU{N) is related to the dilaton and 
i?2 appearing on SUGRA action by [53] [S] 




(2.21) 



For brane configuration in Fig 2.2, the dual geometry has no B2 and dilaton is 
constant with its value set by the string coupling i.e. e""^ = I/qs which gives = 
gl = Sirgg. Thus from dual supergravity we conclude that the gauge couplings do 
not run with scale and we have conformal field theory with large t'Hooft coupling 
Ai = A2 = SirggN ^ 1. This is exactly consistent with our analysis of the gauge 
theory of KW model in the previous section where we found that the gauge couplings 
always reach conformal fixed point surface of Fig 2.5. 

The field theories arising from the brane setups in Fig 2.2, 2.6 and 2.11 are at zero 
temperature. Introducing temperature in the field theory boils down to introducing 
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black holes in the dual geometry. For the KW model, metric for the dual geometry 
incorporating finite temperature of the field theory is given by 

ds"^ = — ^df + \dlt'^ + —dr'^ + r'^\/hdsii i 
^/h h g 

9 = l-'f, (2.22) 

where r/j is the horizon of the black-hole. Note that the above metric f l2.22p can be 
obtained by minimizing the SUGRA action f l2.19p with only F5 7^ and every other 
source zero. Thus metrics in (I2.20p and fl2.22p can be obtained from the same action 
with the same sources. 

Next for the brane setups in Fig 2.6 and 2.11, the most general SUGRA solution 
takes the form [57] 

ds^ = (-9idt^ + dx^ + dy'^ + dz^) + v^f^^Mr^ + r'^dM'l\ (2.23) 

where Qi are functionj^ that determine the presence of the black hole, h is the lOc? 
warp factor that could be a function of all the internal coordinates and dAi\ is given 
by: 

dMl = hi{d'ip + cos 61 d(j)i + cos 62 d(j)2f + h2{del + sin^6'i d(j)l) + 

+hi{h^dO\ + sin^6'2 dcf)^ + cos ip {d9id92 — sin 6'isin 92d(f)id(f)2) + 

+ /i5 sin (sin 9i d92d(f)i - sin 6*2 d9id(l)2) (2.24) 

with hi being the six- dimensional warp factors. The advantage of writing the back- 
ground in the above form is that it includes all possible deformations in the presence 
of seven branes, fiuxes and other localized sources in the theory. The difficulty how- 
ever is that the equations for the warp factors hi are coupled higher order differential 
equations which do not have simple analytical solutions. The original KS solution 
without black hole and seven branes is obtained in the limit 

hs = 9i = 1, hi = fixed (2.25) 



^They would in general be functions of {r,9i). We will discuss this later. 
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In the presence of seven branes we obtain: 

= 0, hs = 1, /i4 - /i2 = a, gi = I (2.26) 

which puts a seven brane in a resolved conifold background with a = constant [5^ 
using the so-called Ouyang embedding [SBj- For a = and the following choice of hi, 

hi = 1, h2 = hi = l, hs = l (2.27) 

the metric takes the following form [SS] 



ds^ = (^-dt^ + dx^ + dy^ + dz^) + Vh\dr^ + r^ds^i,i 



\ , 39sM\ j,^^9sNff^ ^'^\_^9sNf^ (.91.92 

1 H logr < 1 H logr H — H ^log sm — sm — 

2'kN ^1 27r V 2/ 47r ^ V 2 2 



(2.28) 



which is Ouyang's solution [38]. In the limit Nj = 0, the above metric f l2.28p reduces 
to the metric of Klebanov-Strassler model [15] . 

As mentioned above, a black hole could be inserted in this background by switching 
on a non-trivial Qi. However a naive choice of fluxes in 02.261) will break supersymme- 
try [M]. In general we will not restrict to dual geometries which are supersymmetric 
even at zero temperature and in particular supersymmetry will be explicitly broken 
by the introduction of black hole. 

Once we introduce black hole, Qi 1 and we do not expect hi to remain constant 
anymore. We also expect M and Nf in fl2.28p to be given by some M^s and N^^ 
respectively. Our first approximation would then be to make the following ansatz for 
the hi,M^s and Nf: 

hi = l + Oig,), h2 = h^ = l + 0{gs), h3 = l + 0{gs) (2.29) 
9 6 

Meff = M + 5] amn{9sNfr{gsMr, Nf = Nf + h^nigsNfTigsMf 

ra>n m>n 

with amn,bmn could in principle be functions of the internal coordinates {4':'Piy9i). 
Note that we have made m > n in the above expansions because the precise limits 
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for which our supergravity solution would be valid are: 

gs, gsNf, glMNf, ^ ] ^ 0, {gsN, gM) ^ oo (2.30) 



These limits of the variables bring us closer to the Ouyang solution with little squash- 
ing of the two-spheres. This also means that the warp factor h in fl2.23l) can be written 
as [571: 




Figure 2.12: Dual geometry of Ouyang-Klebanov-Strassler model with black hole. 



Now for the black hole factors, what are the choices for gi that are consistent with 
taking back reactions of all the sources in our setup? The Einstein's equations derived 
from the action (12.191) for the KS and OKS model (where sources are functions of 
internal coordinates) fixes gis to be functions of {r, 61,62) ■ Our ansatz therefore is 

4 4 

g,{r, 6u ^2) = 1 - ^ + OiglMNf), g^ir, 6^,62) = 1-^^ + 0{glMNj) (2.32) 

where is the horizon, and the (6*1, 6*2) dependences come from the 0{glMN f) 
corrections. The resolution parameter a is no longer a constant but a function of 
horizon radius and number of D5 and D7 branes i.e. a = a{rh) + 0{glMNf). In Fig 
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2.12 we sketched the dual geometry of OKS model with a black hole. The geometry 
is a resolved deformed conifold with squashing of the two S^'s and as the local action 
for D7 branes appear explicitly in fl2.19p . we have sketched their embedding in the 
dual geometry. Note that D3 and D5 branes contribute to the geometry through the 
fluxes but do not explicitly enter into the action. Hence the these branes are not 
sketched in the dual geometry and only D7 branes form the localized sources. 

But what are the fluxes in fl2.19p that give rise to a metric of the form (12.231) 
with Qi as in (I2.32p ? The background RR three and five- form fluxes can be succinctly 
written as: 

= dr Ae^ A (ci d9i + C2 d62) + dr A (ca sin 6i dOi A d(f)i — C4 sin 62 d62 A d(f)2) 

CT^ -\- Qo? T \ 
ci sin 6*2 d(t)2 C2 sin 61 dcpi A dOi A d62 , 
2r 2 J 

F3 = dr A A (ci sin 61 d(j)i + C2 sin 62 d(j)2) 

9s 

H — A (cs sin 61 dOi A dcpi — cq sin 62 d92 A d(j)2) 
9s 

-— sin 01 sin 62 (- C2 dOi - ^ ci dOA A d(pi A d(j)2 ■ (2.33) 

9s V2 2r ; 

where H-^ is closed and = — CqH^^, Cq being the ten dimensional axion. The 
derivations of the coefficients appearing in (I2.33P are rather involved and can be found 
in [57] y. Here we just quote the results for a constant resolution parameter a: 



c, = ^y^^ (72a4 - 3r^ - 56a logr + a^Hogir^ + 9a^)) cot ^ 

4:7rr{r^ + Qa'^Y 2 

C2 = ^^l^^f (r^ - 9a^ log(r^ + 9a^)) cot % (2.34) 



^Observe that the background EOMs cannot be trivially worked out by solving SUGRA EOMs 
with fluxes and seven branes sources. This is because, even if we know the energy momentum tensors 
for the fluxes, the energy momentum tensors for Nf coincident seven branes are not known in the 
literature. In particular non-abelian Born-Infeld action for Nf seven branes on a curved background 
is unknown. In the absence of such direct approach, we use an alternative method to derive the 
EOMs. This method uses the ISD (imaginary self-duality) properties of the background fluxes and 
flelds. Details on this appears in [38] [54] for the case without black hole . For the geometry with a 
black hole we found [57] that one could find consistent solutions to EOMs using similar arguments. 
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— 36a^ — 36r^ log a + 34r^ log r 
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^2MiV. 
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247r(r2 + 6a2) 
glMN 



ISa^ - 36(r2 + 60^) loga + S{2r'^ - 90^) logr 

+ (10r2 + 63a2)log(r2 + 9a2) 



247rr2 



- 36a2 - 36r2 log a + IGr^ log r + (lOr^ + Sla^) log(r^ + 9a^ 



with K 



'^2"!"a°2 ■ All the above coefficients have further corrections that we will 



discuss later. Finally, this allows us to write the NS 2-form potential: 
B2 = Ih (r) cot ^ dOi + 62 (r) cot y rf^s ) A 



(2.35) 



+ 



47r ^ ^ y 2 2 y J 

ilMI^ (-36a2 + 9r2 + IGr^ logr + log(r^ + 90^)) log fsin — sin — j + bi{r) 
X sin 6'2 o?^2 A d02 



with the r-dependent functions 



6i(r) 
&2(r) 



^?2MiV^ 



247r(r2 + 6a2) 



(ISa^ + (16r2 - 72a2) logr + (r^ + 90^) log(r2 + 90^)) 



57rr^ 



r^ + 9a^) log(r^ + 9a^ 



(2.36) 



and b^ir) and &4(r) are given by the first order differential equations 



Kir) 



+ 



glMNf 



SgsMr 

r2 + 9a2 ' 87rr(r2 + 9a2) 



-36a^-36a^loga + 34r2logr (2.37) 



+ (10r2 + Sla^) log(r2 + 90^ 



3gsM{r^ + 6a^) g^MNf 



„2 I Q/5„2^ 



18a2-36(r2 + 6a2)loga 



(2.38) 



+ (34r2 + 36a2) logr + (lOr^ + 630^) log(r^ + 9a^ 

Putting back the forms of Cj in fl2.33p we can see how exactly the fluxes change 
with conifold coordinates, but being quite involved, it is hard to compare with the 
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Ouyang model. However there exist an alternative way to rewrite the fluxes which 
would tell us exactly how the black hole modifies the original Ouyang setup. This 
can be presented in the following way: 

^gsNf 



F3 = 2MAi (1 + log r) 64, A- (sin 9^ dOi A dcjyi - Bi sin 62 (102 A 

3gsMNf dr / ^2 . , ^1 . , \ 

: -A2 — A A cot — sin 62 d(p2 — B2 cot — sin 9i dcpi 

A-K r ^ \ 2 2 / 



^2 



^QsMN 



Stt 



( Q 9 \ 

^ A3 sin Oi sin 62 I cot y dOi + B3 cot y I A rf^i A d(f)2 (2.39) 

/ 9gsNf , QsNf , 9i 92\dr 

Hs = 6(7,A4M 1 + ^^logr + ^logsm-^sm^ — 
\^ 47r 2tt 2 2 J r 

A^ (sin 6*1 dOi A d(f)i - B4 sin 6*2 rf6'2 A rf02 ) + ^^iMI^l^^ [ ^ A - ^c?e^ ) 
2y y87ryr2y 

A^cot y rf^2 - B5 cot y c/^i^ 

where we see that the background is exactly of the form presented in [38] except that 
there are asymmetry factors Ai, Bi. These asymmetry factors contain all the infor- 
mations of the black hole etc in our backgrounqj. To order 0{gsNf) these asymmetry 
factors are given by: 

Ai = l+^-^-^-{2-3\ogr) + 0{aYsN',) 

B2 = 1+ ,^^fo^2\' +0{a'g!N') (2.40) 
r^ + lSaVlogr ^ ^ ^ 

= 1 + ^ ■ log r + 0{a'g',Nf) 

- ^ + y ■ Anr- + 9gsNfa\2 - 3 log r) + ^^^^^ 
A3 = 1 - ^ ■ log r + Oia'glN'f) 



„ ^ 36a^log r n 9,^0, 

= 1+ 2 ig 21 + a2^?'iVj 

— ISa^log r ■' 

- + C(a^^^A^I), B4 - ^ ^ ^2 _ 3^2 ^ 



A4 = 1 - ^ + 0(a^^,^iV^), B4 = 1 + + Oia'gjN'^ 



■^One can easily see from these asymmetry factors that one of the two spheres is squashed. This 
squashing factor is of order 0{gsNf) and therefore could have a perturbative expansion. Note also 
that although the resolution factor in the metric is hidden behind the horizon of the black hole the 
effect of this shows up in the fluxes. As far as we know, these details were first considered in [57] 
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A5 = 1 + ^ + Oia'g!N]), B5 = 1 + — T^TT^T + O a^^^^ivj 

r r + 3oa^iog r 

These asymmetry factors tell us that corrections to the Ouyang background [38] come 
from (9(a^/r^) onwards. Thus to complete the picture all we now need are the values 
for the axio-dilaton r and the five form F5. U zi = fi gives the location of a single D7 
brane as in fl2.16p . then from F theory one obtains r ~ log(zi — /i) |38][58] near the 
D7 brane. In section 2.4 we will discuss in some detail how one obtains this form of 
T. As r = Co + ie"*^ we get the following form for the dilaton and axion: 
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+ *{d'^x A dh'^) 






(2.41) 




9s 









where 0{fi, a) denotes all orders in fi and the resolution parameter a, h is the ten 
dimensional warp factor discussed above. Thus combining fl2.96p and fl2.4ip our 
background can be written almost like the Ouyang background [38j with deviations 
given by (12.401) . 

In order to extract temperature from the geometry, we look at the metric in 
(I2.23P in the near horizon limit r ^ r^. To be exact, we really need to start from 
the ten dimensional supergravity action and then integrate out the internal directions 
to obtain a five dimensional effective action. Minimization of that five dimensional 
effective action will give the five dimensional effective metric which is the same as 
integrating the five dimensional metric over the internal directions. That is the five 
dimensional metric is given by g^^ = J d9id(f)idipGfj,u{0i, 4>i, ip) where Gfj,^, /i, = 0, .., 4 
is the metric as in (I2.28P and we only integrate over the internal coordinates that Gf^,y 
is a function of. We will denote the resulting warp factor for the five dimensional 
metric g^i, as h{r) for the ensuing analysis of temperature. 

Now, looking at the r, t direction of the metric g^i, and by change of variable. 
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under the assumption that / dOkdcpjdip Qi ^ g{r), k = 1,2, we can define as: 



4:\/h(rh)q(r) 

P" = \ u 12 (2-42) 

so that the near horizon hmit of five dimensional effective metric takes the following 
Rindler form: 

ds' = -P \'}"f, J tl + dp- (2.43) 
where prime denotes differentiation with respect to r and we only wrote the r, t part of 



the metric in terms of new variable p and tc = y 9i''"c)t- The reason behind rescaling 
time at fixed is that with this time coordinate tc, the five dimensional metric 
induces a four dimensional Minkowski metric at every re- 

Now the temperature observed by the field theory with time coordinate tc can be 
extracted by writing the metric in fl2.43p in the following form 

ds'- = -Att'-T^p'- dtl + dp- (2.44) 

Thus comparing fl2.43p and fl2.44p . we obtain the temperature as: 

= f^'^''^ (2.45) 

In the limit where we have (^i = (72 = fl' = 1 ~ 7!, we can easily compute the 
corresponding temperature using the above formula fl2.45p . This is given by: 

7; = ^ + + E CmnJ^^^^ ^T, + 0{l/rc) (2.46) 

where L is defined earlier, Cmnp is in general functions of {gs,M,N,Nf), and > 
Tdeconf (where Tdeconf is the deconfinement temperature) is the temperature at rc 00 
i.e 

n = Tboundary = ^i^L, = F{T,) = T (2.47) 

47c^h{rh) 

where F{Ti,) can be obtained by inverting the first equation above. We can do this 
exactly once we know the black hole factor g{rh) as well as the warp factor h{rh) to 
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all orders in gsNj, ggM . Here we identify the black hole horizon radius Vh with what 
characteristic temperature T which is the only scale in the theory when there are no 
D7 branes. 

With the knowledge of the dual gravity for the branes setups of KS and OKS mod- 
els with temperature, one may ask whether Seiberg duality cascade can be realized 
from supergravity. The answer lies in identifying the gauge theory effective degrees 
of freedom using gravity. For AdS/CFT correspondence this identification is rather 
simple; the AdS throat radius L is directly related to the number of colors of the 
gauge theory by A^ = ^^g^^,2 as can be seen from (I2.2p . Generalizing this result for 
the non- AdS /non- CFT models, one can immediately identify the effective degrees of 
freedom of the gauge theory with the throat radius Z/'^(r) = r'^h where h is the warp 
factor appearing in (I2.23p . This gives 



^eff(r) = #4 = (2-48) 

where is a constant which depends on the five dimensional compact manifold M.^. 
For example, if = T^'\ then = f tt [33]. 

Now with the warp factor h given by fl2.3ip . we see that Aefr grows with r and 
this gives a flow of effective degrees of freedom with changing r of the dual gravity. 
More precisely one can introduce a cutoff Vc in the geometry and evaluate the four 
dimensional dual field theory on the boundary r = re- The bulk geometry obtained 
this way extends from to and the dual gauge theory has degrees of freedom 
A'ofr(rc). Note that the radial coordinate of dual gravity is identified as the energy 
scale of the gauge theory, that is r = A and the geometry with cutoff is dual to 
effective field theory at scale Ac- Of course one has to be specially careful in obtaining 
effective gauge theory Lagrangian from dual gravity. In particular one has to attach 
geometries from r = rc to r = oo in order to account for appropriate irrelevant and 
marginal operators that have been integrated out. This procedure of attaching UV 
caps to geometries will be discussed in detail in the comings sections. 

As Aeff grows with r, we conclude from gravity that effective degrees of freedom 
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PARENT CASCADING THEORY 




Figure 2.13: Seiberg duality cascade captured by dual geometry. 



grow in the UV and shrink in the IR, which is exactly consistent with Seiberg duality 
cascade. At some UV scale Auv, we start with SU{N + M) x SU{N) and as we go 
to lower and lower energy scales, we can describe the theory with lower and lower 
rank gauge groups using a cascade of Seiberg dualities. Thus supergravity captures 
the generic feature of Seiberg duality cascade and N^s at some value of of r, that is 
at some value of scale A describes the overall degrees of freedom of the gauge theory. 
Fig 2.13 shows a simplified depiction of how an effective field theory at scale A^ can 
be mapped to a geometry with cutoff r^. For more details consult \57\ . 

Although generic features of the duality cascade is captured by gravity, there are 
crucial distinctions between the cascade and the flow obtained from dual gravity. To 
describe a theory with its Seiberg dual, a key requirement is that the coupling has 
to be either very strong or very weak. In the space of coupling constants, one must 
be at the boundary curve of Fig 2.14. For A^, M ~ the t'Hooft couplings are 

roughly equal to the gauge couplings, that is ~ gk so we can treat gk as the t'Hooft 
couplings. Then Seiberg duality can be performed at a scale when t'Hooft coupling 
Afc of SU{N + M) is very large and Xk-i of SU{N) is very small. On the other hand 
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Figure 2.14: Effective field theories at various scales along with their dual geometry 
with cutoffs. 

supergravity solution is valid when N, M ^ 1 and using fl2.2ip one obtains that both 
t'Hooft couplings Xk,Xk-i are very large. This means the SUGRA solution lies in 
the interior of coupling constant surface of Fig 2.14 and not on the boundary curve 
where Seiberg duahty holds. The situation is depicted in Fig 2.14. 

The regime where supergravity solution holds is precisely beyond the phase space 
where Seiberg duality holds. The flow captured by supergravity is a smooth RG flow 
with continuously changing degrees of freedom rather than a step by step reduction 
of gauge group by some units (for KS gauge theory by M units) as described by 
the cascade. However supergravity may capture crucial features of the gauge theory 
which are independent of smooth RG flow or patchy Seiberg duality cascade. 

With a clear understanding of the gauge theories arising from the branes and 
their dual gravity, we would like to know how to compute relevant gauge theory 
observables. In the following section we discuss the precise mapping between gauge 
theory correlation functions and geometry. 
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2.3 Gauge Theory Observables from Dual Gravity 



In order to compute expectation values of certain observables of the gauge theory, one 
first computes the partition function. For a system at thermal equilibrium, various 
thermodynamic quantities such as free energy, pressure, entropy etc. can easily be 
obtained from the partition function. On the other hand propagators of quantum 
fields and subsequently scattering amplitudes of particles can be easily evaluated 
by taking the functional derivative of the partition function with respect to some 
appropriate sources. Thus once we know the partition function of a quantum field 
theory, we can essentially determine its thermodynamic properties. 

As the Hilbert space of certain quantum field theories arising from brane excita- 
tions is contained in the Hilbert space of the geometry sourced by the branes, there 
should be a one-to-one correspondence between the partition functions of the gauge 
theory and the dual gravity. A precise mapping was proposed by Witten [59] and sub- 
sequently by the authors of |;60j . Witten's proposal states that the partition function 
of the strongly coupled quantum field theory should be identified with the partition 
function of the weakly coupled classical gravity. In particular if we are interested in 
computing the expectation value of an operator < O > with a source (po, one can 
make the following identification of the partition function as a functional of the source 

-Z^gaugef^o] = (exp / (f)oO) = 2^gravity [0o] 

= exp(S'suGRA[0o] + 5'gh[0o] + 5'counterterm[0o]) (2.49) 

where is a Minkowski manifold, S'gh is the Gibbons-Hawking boundary term 
[6Tj . 00 should be understood as a fluctuation over a given configuration of field and 
^counterterm is the couutcr-term action added to renormalize the action. We will briefly 
describe how each term is obtained from supergravity action and their necessity. 

First observe that is a bulk field, which is in general a function of the coordinates 
of the geometry and 0o is the boundary value of 0. One integrates the classical 
ten dimensional gravity action (12.191) over the compact manifold A^s to obtain an 
effective action Sf^ for a five dimensional manifold. For AdS/CFT correspondence 
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the five dimensional manifold one obtains after integration over 5 = 6*5 is of course 
AdS^ geometry whereas for non-Ads geometries described in the previous section, 
one obtains modified AdS space. The bulk five dimensional geometry has a boundary 
at r = 00 and the boundary describes a four dimensional manifold with an induced 
Minkowski flat metric. Thus the boundary value of 

(j)o{t, X, y, z) = 0(r = 00, t, X, y, z) (2.50) 

is a field which lives in four dimensional fiat space and plays the role of source in the 
four dimensional gauge theory. 

Now to obtain the gravity action as a functional of 0O) one integrates the five di- 
mensional effective action Sf^ over the radial coordinate r and we are left with gravity 
action only being function of 0o- Given a boundary value 0o, one can uniquely build 
the bulk field and write the bulk action for as a function of the boundary value 
00- This was shown first by Witten [59j for AdS bulk fields and for the geometries in 
[571 162] which are perturbations on AdS space, the same argument holds. 

It turns out that in general for the gravity action S'[0] to be stationary under 
perturbation — t- 0+50, one needs to add surface terms which are known as Gibbons- 
Hawking terms [6IJ. On the other hand, after the radial integral is done, for some 
sources 0, S'[0o] becomes infinite and one needs to regularize the action to obtain 
finite expectation values for operators. The regularization requires addition of extra 
terms which cancel the infinities that appear in the gravity action and are denoted 
by Scountertemi- Keeping all this in mind, one takes the functional derivative of (I2.49P 
with respect to 0o to obtain 



^ Q ^ — ^^gaugc _ ^ ^gravity 

500 500 

~ exp (5totai[0oJ) 7- Uo=o (2.51) 

Ocpo 

where S'total[0o] = 'S'suGRa[0o] + 'S'gh[0o] + •S'countortorm [0o] • 

Observe that in the usual AdS/CFT case we consider the action at the boundary 
to map it directly to the dual gauge theory side. For general gauge/gravity duali- 
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ties, there are many possibilities of defining different gauge theories at the boundary 
depending on how we cut-off the geometry and add UV caps. We will elaborate the 
addition of UV caps in section 2.3.2 and in the following section 2.3.1, we briefly 
describe renormalization of the supergravity action. 

2.3.1 Holographic Renormalization 

As already mentioned, supergravity action S'sugra[0o]) for some sources (^q, becomes 
infinity and one needs to regularize the action. In particular if we want to compute 
the stress tensor T^"?, p,q = 0, 1, 2, 3 of a field theory in four dimensional fiat space 
time, then the source 0o is the metric ripg, i.e the Minkowski metric. It turns out that 
if we naively integrate over the radial coordinate r of S^^, the resulting action as a 
function of rjpg blows up. 

First observe that it is indeed possible to write the five dimensional supergravity 
action for AdS or modified AdS space as a functional of fiat four dimensional metric. 
This is because for both AdS and non-AdS space (KS and OKS background with 
black hole), the five dimensional metric has the following form: 




= g^Jx^'dx" (2.52) 

At the boundary r = oo, (^i — )■ 1, and the induced metric is rjpq/ ^ h{oo) which is 
a constant factor times the Minkowski metric. Thus Sf^[gfj,,y] can be written as a 
function of rjpg once the radial integral is done. 

Now for AdS space, starting with f l2.19p with only ^ being the source, 
integrating over the internal coordinates give the following action: 

Sfis, = I d'x r dr^ (i?5 + A) (2.53) 

where g = detf^^^y, A is the cosmo logical constant and — )■ oo. Note that the above 
action is not only the action for AdS space but also for asymptotic AdS spaces i.e. 
geometries such that limr_j,oo G^j, = g^j^^ with g^j^^, given by (I2.52p where h = L'^/r^ 
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After the r integral is done, the result contains terms which diverge as Vc — > 
oo. As the AdS or asymptotic AdS metric is a power series in r, the divergent 
terms are precisely the ones containing positive powers of Vc and by subtracting 
these terms we get a finite result for the action. This subtraction procedure can be 
done by introducing counter term S'counter to the action fl2.53p such that they cancel 
the infinities coming from it. The renormalization scheme mentioned here was first 
proposed by Skenderis and for more details please consult [63]-[67]. 

With a finite renormalized action S^^^g^^ [rjpq] , and using the equation of motion for 
the metric g^,^, one can take the derivative with respect to rjpg to obtain the following 
result for the stress tensor [65] 

T^'^ = (^(4)pg - ^9{o)pg [(Tr^(2))^ - ^^9(2)^] - ^(^(^2))pg + ^c/(2)pgTr^(2)) (2.54) 

where is the five dimensional Newton's constant and is defined through the 
bulk metric G^u the following way: 

G.^dz'^dx" = ^ + -g„Jxf',p)dxPdx'^ 
4p^ p 

9{x,p) = 9io)ix) + ... +p'^/2^(d) + .. (2.55) 

where p = Thus are Taylor coefficients in the expansion of the four dimen- 

sional metric Qpg^x^.p) which induces the four dimensional boundary metric g{Q)pq at 
the boundary p = 0. 

The story gets somewhat more involved for non-AdS geometries and we will briefly 
review the key points of holographic renormalization of dual supergravity for non 
conformal field theories. Details of the analysis can be found in our work [SZl 

The divergent terms appearing in the supergravity action are uniquely determined 
once we know the form of the warp factor h in (12.231) . Observe that the logarithms 
appearing in h fl2.31l) results from the logarithmic running of the dilaton field fl2.41l) 
and that of NS-NS field i?2- The dilaton only behaves logarithmically near the location 
of the D7 brane and its global behavior is determined by F theory. For large r and 
away from any D7 branes, we expect the dilaton to behave as a constant. In fact in 
section 2.4 we will analyze the precise running of the axio-dilaton field for all r where 



40 Chapter 2 The Gauge/ Gravity Correspondence 



we will also see how the logarithmic running of NS-NS B2 field can be modified to 
give finite value for B2 at large r. With B2 and axio-dilaton r behaving as logarithms 
in the local neighborhood of some r but approaching constant finite values for large 
r, in [57] we proposed the following form of the warp factor: 



^ - r^-ei + r^-2e, r^-e, + ^4_^.i/2 - A. ^4^^ (2-56) 



9sM 2Nf 



TV \ N 



where ei,r(a) etc are defined as: 

ei = — - H log sm — sm — , eo ■ 

27rN 87rm SvriV ^\ 2 2 / ' 

2 

1-e, ^ ^1 ^2 ^2 ^2 

r(,) = r e(i) = -, e(2) = -, e(3) = -, e(4) = - 

r(±„)=r^^^w, L(i) = L(2) = L(4) = L^ L(3) = -2L^ (2.57) 

which makes sense because we can make to be very small. Note that the choice of 
ei doesn't require us to have gsNf small (although we consider it here). In fact we 
can have all (A^, M, Nf) large but small. A simple way to achieve this would be to 
have the following scaling behaviors of {gs, N, M, Nf): 

gs e", M ^ e-^, Nf ^ e"^ (2.58) 

where e — )■ is the tunable parameter. Therefore all we require to achieve that is to 
allow: 

a + 7>2/3 + K, a > K, ^ > a (2.59) 
where the last inequality can keep gsNf small. Thus gsN, ggM are very large, but 



9s, ^^,gsNf are a. 
gravity background 



1 very small to justify our expansions (and the choice of super- 



''For example we can have gs going to zero as gs — > e^/^ and {N,M,Nf ) going to infinities as 
{e^^,e~^,e~^) respectively. This means {gsN, gsM) go to infinities as (e""'^/^, e""'/^) respectively, 
and {gsNf, g^MNf, gsAP /N) go to zero as (e^^^, e, e^^^) respectively. This is one limit where we can 
have well defined UV completed gauge theories. Note however that for the kind of background that 
we have been studying one cannot make Nf large because of the underlying F-theory constraints 
[SH [58] [82] . Since we only require gsNf small, large or small Nf choices do not change any of our 
results. 
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The warp factor f l2.56p has a good behavior at infinity and reproduces the 0{gsNf) 
result locally. Our conjecture then would be the complete form of the warp factor at 
large r will be given by sum over a as in f l2.57p but now a can take values 1 < a < oo. 
This conjecture in fact justifies the holographic renormalizability of our boundary 
theory as we will find out shortly. 

With the warp factor given by fl2.56p . we introduce perturbation Z^,^ to background 
five dimensional metric Qf^u given by f l2.52p . Recall that for non-AdS space the five 
dimensional effective action S^^ derived from ten dimensional action fl2.19p takes the 
form 

Sf = J d^xVG {R + A(G'^.)) (2.60) 

where G^,^ = (7^^ + l^i, and the source A(G^,^) is in general a function of the metric. 
As l/j^u is a perturbation, one can write the action fl2.60p as a Taylor series in the 
perturbation. Considering terms only up to quadratic order in 1^^,^ one obtains: 



5(1) [$] = / y [ dr {^Ar(r,g)[$!ikr,g)<W(r,-g) + $:;l^kr,g)$W(r,-g) 

+i?r(r, g)$l'l(r, g)$f (r, -q) + ^Cr(r, q) [<f'W(r, q)^l\r, -q) + $ W(^, q)^f\r. -q) 



+Dr{r,q)^^^{r,q)^^:;\r,-q) + Tr{r,q)^^^{r,q) + ET^'^^\r^ 



where m,n = 1,...,5, prime denotes differentiation with respect to r, the script [1] 
denote the total background to 0{gsNf,gsM'^/N); and the explicit expressions for 
j\rnn^ j^mn ^ (jmn ^ pm ^ specific casc are givcu in Appendix D of [57]. We have also 



(2.61) 



defined $W(r, g) in the following way (with go = ^y9{f^c) as before): 

$W(r,g) = / ^-^= e'^'^'^'-'^^^-'^"^'n^m{t,r,x,y,z) (2.62) 

J (27r)4y^ 

We will see that the effective four dimensional boundary action is independent of 
Z)]"" and and hence their explicit expressions do not appear in the appendix of 
[57] . Furthermore, note that the derivative terms in f l2.6ip all come exclusively from 
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\/—GR^ whereas fluxes contribute powers of $1^' but no derivative interactions. In 
tfie following we keep up to quadratic orders, and therefore the contributions from 
the fluxes will appear in and -E™. 

The equation of motion for $ W (r, —q) is given by: 



Dr{r,q)^l\r,-q) + ^^rir, g)<W(r, -g) - hor{r,qW^\T,-q) 



+ mr^q) - E'r{r,q) + Fr(r,g) = 



(2.63) 



The next few steps are rather standard and so we will quote the results. The variation 
of the action (I2.6ip can be written in terms of the variations 5$ W (^f^ and 5$ W (r, —q) 
in the following wajQ: 



2^ (271^ ^g^)Jr, 



dr< 



[1] 



[1] 

m 



(2.64) 



which includes the equations of motion as well as the boundary term. We can then 
write the variation of the action SS^^^ in the following way: 



1 r 
+ 2 



dr 



[1] 



^ , , ^, (EOMfor$W)5<l>W + (EOMfor$W)5<l>„ 

(25r - Ar)(<^'l''5$L'' + -^-f 5$„) + (Cr - Ap)($W(5<l>W + $W5$W) 



boundary 



(2.65) 



^Henceforth, unless mentioned otherwise, $m , will always mean $[„ (?", <z) and (r, — g) 
respectively. Similar definitions go for the variations 6^m and 6^n^ . 
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where by an abuse of notation by the "boundary" here, and the next couple of pages 
(unless mentioned otherwise), we mean that the functions are all measured at Vh and 
Tc i.e the horizon and the cut-off respect ivel}o. It is now easy to see why D^^, 7^"*" 
and i?^" etc do not appear in the boundary action. Finally, we need to add another 
boundary term to fl2.65p to cancel of the term proportional to This is precisely 
the Gibbons-Hawking term 
1 r d^q 



4™"<|)W<T)'[1] I 4mn<T,[l]<T)'[l] I 2F-,"$' 

1 rri n '1 n m ' 1 r 



[1] 



(2.66) 



boundary 



Taking the variation of f l2.66p 5/Ci we get terms proportional to 5$'^^' as well as 5$'^'. 
Adding 5/Ci to 5S^^^ we can get rid of all the 5$''^] terms from fl2.65p . This means we 
can alternately state that the boundary theory should have the following constraint 

(5$lil(r„g) = 5<l>'W(r„-g) 







(2.67) 



With all the above considerations we can present our final result for the boundary 
action. Putting the equations of motion constraints on fl2.65p . as well as the derivative 
constraints fl2.67p . we can show that the variation f l2.65p can come from the following 
boundary 3 + 1 dimensional action: 



5(1) 



+ 



, , , , , ,C^r(r,g) -^r(r,g)]$S(r,g)$W(r,-g) 

(27r)Y(7(r^ax) I' 
Sr(r, q) - Ar{r, q)] [<l>'W(r, g)$W(r, -q) + $W(r, q)^i'\T, -q) 



+ (E^-F[^)^^^{T,q) 



(2.68) 



However the above action diverges, as one can easily check from the explicit expres- 
sions for A^^", 5]"", C{"", E"^ and for the specific case worked in Appendix D 



^Note that we have not carefully described the degrees of freedom at the boundary as yet. 
For large enough we expect large degrees of freedom at the UV. This would mean that the 
contributions to various gauge theories from these degrees of freedom would go like e^^"" , which 
would be negligible. Thus unless we cut-off the geometry dX r — rc and add UV caps with specified 
degrees of freedom we are in principle only describing the parent cascading theory. For this theory 
of course A/'ofr is infinite at the boundary, which amounts to saying that UV degrees of freedom don't 
contribute anything here. We will, however, give a more precise description a little later. 

^One can impose similar constraints at the horizon also. 
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of [57]. Indeed, comparing it to the known AdS results, we observe that from the 
boundary there are terms proportional to in each of C{"'^, A']^"", E'^ and -F{™ and 
proportional to in A'™,!?]"". As Tc — )■ oo the action diverges so as it stands Vc 
cannot completely specify the UV degrees of freedom at the scale Ac. Thus we need 
to regularize/renormalize it before taking functional derivative of it. This renormal- 
ization procedure will give us a finite boundary theory from which one could get 
meaningful results of the dual gauge theory. 

Once we express the warp factors in terms of power series in r(Q,) fl2.56p the renor- 
malizability procedure becomes much simpler. Note however that this renormaliza- 
tion is only in classical sense, as the procedure will involve removing the infinities in 
(I2.68P by adding counter-terms to it. Comparing with the known AdS results, and 
the specific example presented in Appendix C of [57] , one can argue that the infinities 
in fl2.68p arise from the following three sources: 

1. Cr{r.,q)-A'r{r,,q) = T.H^\i<l)<ia) + finite terms 

a 

2. 5r(rc,g)-Ar(^c,g) = E^lT(?)^c(a) + finite terms 

a 

3. Er(r„g)-Fr(r„g) = E^H + finite terms (2.69) 

a 

where in the above expressions we are keeping a arbitrary so that it can in general 
take both positive and negative values; and the finite terms above are of the form 
^c(a) with n>l. Therefore to regularize, first we write the metric perturbation also 
as a series in l/r(^a)'- 

^1^^ = EE7^ (2-ro) 

fc=0 a ' (a) 



where the above relation could be easily derived using (I2.62p . taking the background 
warp factor correctly. 

Plugging in fl2.70p and f l2.69p in (I2.68P we can easily extract the divergent parts 
of it . Thus the counter-terms are given by: 



c(l) _ / " V r rrm, 

'-'counter — / — ^ / ^ = / . 



I o{0)[/3]„(l)[7V3 \ , / (2)[/3] (0)[7]»,2 , (0)[/3] (2)[7U2 , „(l)[/3] „(1)[7]^2 ^ 
' mm nn (02) / \ mm nn (03) mm nn (04) mm nn (0^5) , 
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^(^mmJOM I „(0)[/3] (3)[7]^^ , s(2)[/3]„(l)[7]^, , , s(2)[7] V, 

'y'^m.m ^nn ' ("6} ~ ^mm "^nn ' (ar) ~ ■^mm ^nn ' (as) ~ ^mm "^nn J' (019) 

^ |a| K^mm '^nn (oio) ' '^mm ^nn (qh)/ \ '^mm ^nn (ai2) 



I Q,(0)[/3] J3)[7]^^ , I Os(l)[/3]s(2)[7]^, , , ,(2)[/3] „(1)[7]^, , , Q s(0)[/3] J3)[7U, 
^'-'"'mm "'nn ' (an) ^ '^•^nn "^mm ' (ais,) ^ ^nn "^mm ' ("19) ^ '-'"'nn "'mm ' ("20)^ 

^-'la^VO '^Jl^'Smm ^(a) + "^mm ^(ai) + "^mm ^(az)+^m'm ^{ae))} l^- ' -^J 

with an equal set of terms with ?"(_„.). In the above expression r(Q.) = r^^^("«'; and 
as before, the integrand is defined at the horizon and the cutoff Vc- The other 
variables namely, sl'^}ll^\ H^J^ , K^" and are independent of r but functions of q\ 
For one specific case their values are given in Appendix C of [57]. Finally the e(a^) 
can be defined by the following procedure. Lets start with the expression: 

4 

o-mn (a)[/3] (fe)[7] p ^ Trrnn (a)m (b)['y] ^(a) 
-'-'lal "^mm *nn ' (a^) — -'-'|Q|'^mm "^nn a 

5 

Tv-mn (c)[/3] (d)[7] q ^ jv-mn (c) (d) [7] ^(«) 79^ 
-''■a *mm *nn '(a,) — -"^ a "'mm ^nn n „rf l^' ' 



from where one can easily infer: 



p = 4-a-b, e(„^) 



4e{a) — ae(/3) — &e(7) 



4 - a - 6 

g - 5-c-d, - 5-c-d ^ ^ 



Using this procedure we can determine all the r(Q-) in the counterterm expression 

At this point the analysis of the theory falls into two possible classes. 

• The first class is to analyze the theory right at the usual boundary where — )■ 
00. This is the standard picture where there are infinite degrees of freedom at the 
boundary, and the theory has a smooth RG flow from UV to IR till it confines (at 
least from the weakly coupled gravity dual). The action obtained by setting Vc = oo 
describes the 'parent cascading theory'. 

• The second class is to analyze the theory by specifying the degrees of freedom at 
generic energy scale given by Vc = Ac and then defining the theories at the boundary 
by adding appropriate irrelevant and marginal operators for scales greater than A^. 
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By adding different operator for scales A > Ac, one gets different theories. All these 
different branches meet the parent cascading theory at some scale Vc = Ac, as depicted 
in Fig 2.13. The gravity duals of these theories are the usual deformed conifold 
geometries cutoff at various Vc with appropriate UV caps added (of course for r < Vc 
the geometries change accordingly). More details on UV caps to geometries will be 
discussed in the next subsection. 

The first class of theories is more relevant for the pure AdS / CFT case whereas the 
latter is more relevant for the non-AdS casqj.For the pure AdS/CFT case without 
flavors e(a^) = (so that the subscript a^'s can be ignored from all variables), we 
can subtract the counter-terms f l2.7ip from the action fl2.68p to get the following 
renormalized action: 

c(i) = c(i) _ c(i) 

^rcn ^ "-^counter 



(2vr) 



rrmnf U) (0) , (3) (1) , (2) (2) , (1) (3) 
y mm nn ' mm nn ' mm nn ' mm nn 

+ 0(0) - /r™"f4s(0) + Ss^^) ^(3) +4s(2) -U s(3) ,.(1) 

' mm nn J y mm nn ' mm nn ' mm nn ' mm nn 

4.(0)5(4) , 3s(i)c(3) I 5(3)5(1) + /"^5(4) 1 (2 74) 

"nn"mm ' ""nn mm ' "nn mm I ' "mm \ ' J 



where we have made all the 0{l/rc) terms vanishing, and in the limit small the 
small shifts to s^{| given by s^^\ + 0{rf/j can also be ignored. Observe that we can 
reinterpret the renormalized action fl2.74p as the following new action: 

Sgl = /^[z-$„(g)$„(-g) + f/-'^($^(g)$;(-g) + <|.;,(g)<|.„(-g) 



(2vr) 



(2.75) 



where $m(g) = ^m{q,r = oo) , $'„(g) = d^rn{q,r)/dr\r=oo and X,Y,Z,U,V could 
be functions of r and if but evaluated for fixed r = oo. We can determine their 
functional form by comparing (I2.75P with f l2.74p . For us however the most relevant 



^In both cases of course we need to add appropriate number of seven branes to get the finite 
F-thcory picture. The holographic renormahzation procedure remains unchanged and the far IR 
physics remains unaltered. The UV caps affect mostly geometries close to Tc, as expected. 
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part is the energy momentum tensors which we could determine from 12. 75] by finding 
the coefficients and F". One can easily show that, up to a possible additive 
constant, Y"^, Y'"" are given by: 

= ff^'^oW _ 4 I^n"^„(4) yn _ TTmn U) _ 4 Z>'mn (4) 

nn ^-L^ '^nn^ mm ^-i-^ ^mm 

ym ^ j^mn^i^)^ = K^^ S^^^m (2-76) 

Now let us come to second class of theories wherein we take any arbitrary r = 
Tc, with appropriate UV degrees of freedom such that they have good boundary 
descriptions satisfying all the necessary constraints. For these cases, once we subtract 
the counter-terms (I2.7ip . the renormalized action (specified by r) takes the following 
form: 

c(i) = c(i) _ c(i) 

ren '-^counter 




I ^(3)[/3]„(l)[7]u4e(^)+e(^)-5e(„) 1 ^3e(^)+2e(^) -5e(a) t\ , rm „(4)[/3] 5e(^)-5e(„) 1 77) 
^^mm ^nn V \) ~ '^\a\'^mm ' j \'^-''J 

evaluated at the cut-off and the horizon radii Vh as usual. Notice now the appear- 
ance of r™''^(")"'""'^('')"'"^'^(T) factors. One can easily show that: 

Since the warp factor h is deffned only for small values of QsNf, QgM'^/N, g'^NfM'^/N 
we don't know the background (and hence the warp factor) for finite values of these 
quantities. Therefore for our case we can put 0[e(^a,f3,'y)]^ to zero so that the value in 
fl2.78p is identically 1. For finite values of these quantities both the warp factor and 
the background would change drastically and so new analysis need to be performed to 
holographically renormalize the theory. Our conjecture would be that once we know 
the background for finite values of gsNf,gsM'^/N, the terms like fl2.77p would come 
out automatically renormalized by choice of our counterterms. 
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Once this is settled, the renormahzed action at the cut-off radius Tc would only 
go as powers of ^c(a)- Thus we can express the total action as: 



5(1) 




j=0 'c(a) 



J 



M-"(<|.^(g)<|.;(-g) 



I T^mn 



1 - 



(2.79) 



r-^^^ |r=rc- Here we have ignored the terms evaluated 



where are as before and r^j-^^ 
at the horizon as they are not proportional to the sources $„ and hence will not 
contribute to the expectation value. The explicit expressions for the other coefficients 
listed above, namely, G"^", M™-", S^^^-^-), J™, e'/nnij) ^^"^ t)e worked out easily 

from our earlier analysis (see Appendix C of [57] for one specific example). Note that 
X[r(Q)] is a function independent of and appears for generic renormalized action. 

Now the generic form for the energy momentum tensor is evident from looking at 
the linear terms in the above action (12. 79 p . This is then given by: 



rpmm 



Tjrnn , rrnmN (4)[/3] _ A(L-rnn , LCJimy (4)1(5] , f ly-mn , /^-^^UCS)!/?] 



1°l 



1°l 



PS 



(2.80) 



\j=0 ' c{a) , 

at r = Tc (we ignore the result at the horizon) and sum over (a, /3) is implied. This 
result should be compared to the ones derived in [63] [M] [65] [66] [67] [68] [69] which 
don't have any Tc dependence. The AdS/CFT result is of course the first line of 
the above result. Notice that the second line also has a Vc independent additive 
constant which are irrelevant for our purpose because the energy-momentum can 
always be shifted by a constant to absorb this factor. Now for non-AdS geometries, 
Tc determines the number of degrees of freedom at scale Ac with the relation (13.331) and 
inverting it with r'^h ~ log(r) we get ~ e^'^"' . Thus once we specify the effective 
degrees of freedom at scale A^, then our result shows that the energy- momentum 
tensor not only inherits the universal behavior of the parent cascading theory but 
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there are additional corrections of Ole'^^"^) . But energy scale should eventually be 
set to infinity, so how is keeping C(l/rc) terms in f l2.80p justified? 

To answer this, first consider a daughter gauge theory for which number of degrees 
of freedom do not change for scales A > Ac and then N^fi^Ac) = A^uv(oo)- Taking 
the energy scale to infinity means replacing in fl2.80p with e^u^''°°-* as only gauge 
theory observables such as number of degrees of freedom should appear in the final 
result for stress tensor. But as Ncsi-^c) = N\jy{oo), for this particular theory, we can 
keep the corrections coming from 0{l/rc) as in fl2.80p . 

This argument of replacing Vc — ?■ 0{e^^^) is somewhat naive and needs elabora- 
tion. If we cut off the geometry at r = r^, and do not add any geometry from r = 
to r = oo, we are essentially ignoring gavitons that are present in the region r > r^. 
This will in turn mean we are ignoring UV modes in the dual field theory that have 
energies A > Ac. Thus fl2.80p is incomplete. What we mean here is once we add 
gravitons from region r > r^, their contribution to the final stress energy tensor will 
can be accounted for by replacing in fl2.80p with (9(e^u^). 

Note that for an AdS geometry, the dual gauge theory is conformal and do not 
change degrees of freedom with scale. Hence treating Ncs{Ac) = A^uv(oo) amounts 
to adding an AdS UV cap to our dual geometry from r = Tc to r = oo. The final 
result for stress tensor of a dual gauge theory with an AdS UV geometry will be of 
the form fl2.80p with Tc — )■ (9(e^u^) which may not be infinity. We will discuss this in 
detail in the following section. 

In general, various gauge theories which are quite distinct at far UV scale can look 
identical in the IR. One can start with a given number of degrees of freedom at a scale 
Ac and demand that it is the infrared limit of several different UV theories. Finally 
specifying the UV degrees of freedom distinguishes the theories. From the gravity side 
we can replace Tc in fl2.80p with diffent values of e^^v obtain distinct UV gauge 
theories which are identical at the infrared scale Ac. In fact a more careful analysis 
shows that one can cut a bulk geometry at some Vc and attach various UV geometries 
from Tc to oo to obtain distinct gauge theories at the boundary. We discuss this in 
some detail in the following section. 
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2.3.2 UV Completion from Dual Gravity 

The first important issue here is that we can study infinite number of UV completed 
theories in full F-theory setup. All of these theories have good boundary descriptions 
and have same degrees of freedom as the parent cascading theory at certain specified 
scales. The simplest UV complete theory is of course the parent cascading theory of 
OKS model. Using (13.331) with warp factor given by fl2.31l) one easily gets 

iVeff(A)~L^(l + 0(log(A))) (2.81) 

which grows with scale A and diverges for A — )■ oo. The question now is how to 
construct other possible theories by defining the degrees of freedom at scale A — )■ oo 
i.e. at the boundary of the geometry r — i- oo. The action for the boundary theory 
should be identified as: 

[S^ll]Z = i^^llYrl + [S^X (2-82) 

where the boundary is at r — )■ oo. For the boundary cascading theory the above 
expression simply means that 



+ ^^-"^-^'^^'^'^^ factors (2.83) 
oo 

where the sign is crucial and sum over a is again implied (note (a) the cut-off de- 
pendence, and (b) rc(^) is some function of rc(a) that one can determine easily). We 
now see that the contributions from the UV cap give the following values for Bj for 
parent cascading theory: 



-°n(0) ~ ^n(l) — -°n(2) " -°n(3) ~ ^ 



B. 



(a) 



1 



n(4) 2 



b%)j"nro - re) + o(i/rjr, i^rjT, 4i 



(2.84) 



The choice above is made precisely to exactly cancel the 0{l/rc) contributions 
coming in from the action measured from Th <t < Tc- Finally the boundary energy- 
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momentum tensor reads 



4E 



(2.85) 



which is the result derived in [63]- [67] and [69] . 

The above way of reinterpreting the boundary contribution should tell us precisely 
how we could modify the boundary degrees of freedom to construct distinct UV 
completed theories. There are two possible ways we can achieve this: 

• From the geometrical perspective we can cutoff the deformed conifold background 
at r = Tc and attach an appropriate UV "cap" from r = Tc to r — )■ 00 by carefully 
modifying the geometry at the neighborhood of the junction point. As an example, 
this UV cap could as well be another AdS background from rc to r 00. There are 
of course numerous other choices available from the F-theory limit. Each of these 
caps would give rise to distinct UV completed gauge theories. 

• From the action perspective we could specify the value of the action measured 
from Tc to r — )■ 00, i.e [5^^en]^- The simplest case where this is zero gives rise to a 
boundary theory for which A^eff(Ac) = A'uv(oo) as discussed before. To study more 
generic cases, we need to see how much constraints we can put on our integral. One 
immediate constraint is the holographic renormalizability of our theory. This tells us 
that the value of the integral can only go as powers of l/rc(a) otherwise we will not 
have finite actions. This in turn implies 



Vi=0 ^c{a) J \j=0 ^c{a) 



00 



B 



(a) 



where by specifying the coefficients ^^^(j), -^mn(i) ^^"^ ^mlj) ^^"^ specify the precise 
UV degrees of freedom! The finite terms are independent and therefore would only 
provide finite shifts to our observables. They could therefore be scaled to zero. Notice 
also that the contributions from (I2.86P only renormalizes the coefficients a^^ri{j)i ^mn(j) 
and in (12.791) . and therefore the final expressions for all the physical variables 
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for various UV completed theories could be written directly from fl2.79p simply by 
replacing the l/vc dependent coefficients by their renormalized values. This is thus 
our precise description of how to specify the UV degrees of freedom for various gauge 
theories in our setup (see Fig 2.15 below). 



finite UV degrees of 
freedom added at r = r 




0123 <— > 



S^X 



UV cap from boundary at r= oo 




mm 



Figure 2.15: The equivalence between two different ways of viewing the boundary 
theory at zero temperature. 



For the left figure in 2.15 we add finite UV degrees of freedom aX r = of the 
deformed conifold geometry. Such a process is equivalent to the figure on the right 
where we cut-off the deformed conifold geometry at r = Tc and add a UV cap from 
r = Tc to r = oo. The boundary theory on the right has M^v degrees of freedom at 
r = oo and all physical quantities computed in either of these two pictures would 
only depend on A/'uv but not on r = r^- At non-zero temperature the UV descriptions 
remain unchanged. 

Once the UV descriptions are properly laid out, we can determine the form for 
anfi i^y writing Callan-Symanzik type equations for them. They tell 
us how A^°'\ B^°''> and E^"'' would behave with the scale Vc or equivalently fic- For A 
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the equation ic 



1 -e 



(a) 



Aa) I ?r(«) 

^mn{j) ^ ^mn{j) 



(2.87) 



with similar equations for i?*-"^ and E'^"\ These equations tell us that physical quan- 
tities are independent of scales. The parent cascading theory is defined as the scale- 
invariant limits of f l2.87p . i.e: 



A 



(a) 

mn{j) 



mn{j)'> mn{j) 



m(j) 



m{j) 



(2. 



The above relation gives us a hint how to express v4("), and E^") in terms of A4fr, 
the effective degrees of freedom at r = Tc and A/^iv, the effective degrees of freedom 
at r = oo i.e the boundary: 



mn{j) 

I g(°) g-i[A/'uv-(l-£{<,))A^off] 
' mn(j) 



■mn(j) 




~(a) 
mn{j) 


mn{]) 




mn(j) 









+ 



(2.89) 



where the actual boundary degrees of freedom are specified by knowing CLyYin{j)i (^mn(j) 
and bm{j) as well as Mnv Since j goes from to oo, there are infinite possible UV 
complete boundary theories possible^"!. For very large Afuv (i-e Afuv — ^ e~^,n » 1) 
the boundary theories are similar to the original cascading theory. The various choices 



~(t),b^rnji)C^)) ^^^^ degrees of freedom change from 



of («mn{i)l 

A/'uv to A/'cff under RG flow. The if' dependence of all the quantities will tell us how 
the UV degrees of freedom affect IR physics. This is to be expected: addition of 
irrelevant operators do change IR physics, but not the far Ilj"]. 



^The following equation is derived from the scale-invariance of [iSrcn] ^ • 
^"The connection of j with UV completions come from the coefficients a'"^^„^ , e^"'' 



and b 



(a) 



rnn{j) ' mn(j) mn{j) 

etc. that depend on j. For different choices of these coefficients we can have different UV completions. 
In this sense j and UV completions are related. 

^^The if dependences of the UV caps are also one-to-one correspondence to the changes in the 
local geometries near the cut-off radius Tc, as we discussed before. All in all this fits nicely with 
what one would have expected from UV degrees of freedom. Of course it still remains to verify the 
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Therefore with this understanding of the boundary theories we can express the 
energy-momentum tensor at the boundary with Afuv degrees of freedom at the bound- 
ary purely in terms of gauge theory variables, as: 



rprnm 



oo 

+ E (2.90) 

j=0 



where sum over a is again implied, and the first line is the universal property of the 
parent cascading theory inherited by our gauge theory. The second line specifies the 
precise degrees of freedom that we add at r = rc to describe the UV behavior of our 
theory at the boundary r oo. Using this procedure, the final results of any physica. 
quantities should be expressed only in terms of A/'uv i-e the UV degrees of freedon^f . 



2.4 Towards Large N QCD from Dual Gravity 

From the above discussions we see how one obtains gauge theories which are quite 
distinct in the UV but have common IR dynamics. With a clear understanding of 
how to UV complete a gauge theory with particular properties in the IR, we can 
construct a brane configuration that mimics large N QCD. 

First observe from fl3.33p and (12.81 1) that in OKS model the number of degrees 
of freedom keep on growing in the UV . Furthermore with B2 ~ log(r) and axio- 
dilaton behaving as in f l2.4ip we see that the gauge couplings run as Qi ~ log(A) 

story from an actual supergravity calculation. We need to analyze the metric near the junction by 
studying the continuity and differentiability of the metric and see how far below r = Tc we expect 
deformations from the UV caps. Various types of deformations will signal various sets of irrelevant 
operators. Needless to say, the far IR physics remain completely unaltered. In the following section, 
we will perform an exact calculation using appropriate sources in SUGRA action accounting the 
deformation in the region near the junction. 

^^Restoring back the 0{rh) contributions would mean that there should be an additional contri- 
bution to (I2J01) of the form Y1T=q G'(&i"\ sk"i)r'*^e-^'^"" where G is a function whose 
functional form could be inferred from the UV integral (|2.86p . 
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using fl2.2ip . In fact in the far UV, at least one of the two gauge couphngs blow 
up and we have Landau poles. Thus the UV of OKS model has runaway behavior 
while QCD is asymptotically a free theory [70] [71] . However the gauge couplings do 
run logarithmically and in the far IR using a cascade of Seiberg dualities one obtains 
SU{M) gauge group with fundamental matter. Hence the IR of OKS model has 
common features with large N QCD while the UV is quite different. But as discussed 
in the previous section, we can modify the UV of OKS model and build dual gravity 
of a gauge theory that resembles large N QCD more closely. 

The ideal scenario would be to build a gauge theory which becomes asymptotically 
free and for scales where coupling is large, it has a description in terms of weakly 
coupled classical gravity. Currently there are no brane setups that realizes such 
running of coupling and has dual gravity in the appropriate regime. What we have 
been able to construct in [62] is a gauge theory that becomes asymptotically conformal 
in the far UV and has logarithmic running of couplings in the IR. The regime where 
the gauge theory has a classical supergravity description is when gauge coupling is 
large and thus from UV to IR we have strong but running gauge coupling. This gauge 
theory has no Landau poles and the number of degrees of freedom do not diverge as 
A — )■ oo. In this section we will briefly describe the brane setup that gives rise to such 
a gauge theory and construct its dual gravity. Details of our model can be found in 

We start with the brane setup of Fig 2.11 i.e. the OKS model. The logarithmic 
running of the couplings come from the global logarithmic running of NS-NS two form 
i?2 ~ log(r) and the local logarithmic running of axio-dilaton near a seven brane as 
in (12.411) . On the other hand the three form fluxes run as 1/r, i.e. similar to the 
scalar potential for a single charge located at r = 0. By introducing anti charge at 
some location r ~ tq, we can get the total three form flux to behave like 1/r^ for 
large r ^ rg i.e. similar to scalar potential of a dipole. This would result in total 
NS-NS two form B2 ~ 1/r for r ^ tq. We can also arrange 7 branes in such a way 
that F theory gives r ~ (9(l/r") for large r ^ tq with n > 0. Details of the 7 brane 
embedding and the global behavior of r will be discussed in section 2.4.3. The anti- 
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charges we introduce are anti five branes in the neighborhood of r ~ ro . Of course 
anti branes want to shde down to the tip of the conifold and we need to introduce 
sufficient flux to keep them at r ~ tq. The resulting brane setup is depicted in Fig 



Figure 2.16: Brane configuration for asymptotic conformal gauge theory. 

The dual gravity for the brane setup in Fig 2.16 considering non-zero temperature 
is sketched in Fig 2. 17. The geometry splits into three regions of interest: Regions 1, 
2 and 3. The figure shows the various regions of interest. As should be clear, most 
of the seven branes lie in Region 3, except for a small number of coincident seven 
branes that dip till rmin i.e Region 1. The interpolating region is Region 2. Region 
1 is basically the one discussed in great details in [57] and the previous sections. In 
this region there is one (or a coincident set of) seven brane(s) along with B2 ~ log(r). 
The logarithmic dependences of the warp factor, fiuxes and subsequently of the gauge 
couplings originate from B2 and these coincident (or single) seven branes . 

The UV cap in the full F-theory framework is depicted as Region 3 in the above 
figure. In this region the seven branes are distributed so that axio-dilaton has the 
behavior i.e. r ~ (9(l/r"),n > while NS-NS two form B2 go as C(l/r™),m > 
. Using these precise forms in f l2.2ip one gets that the gauge couplings go as gi ~ 
C(l/r™) = (9(1/A™), i.e. we have asymptotic AdS space giving almost conformal 
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field theory in the far UV. 

One big advantage about our UV cap is related to the issues raised in [72]. Since 
the Hns and the axio-dilaton fields have well defined behaviors at large r, there would 
be no UV divergences of the Wilson loops in our picture! Therefore our configuration 
can not only boast of holographic renormalizability, but also of the absence of Landau 
poles and the associated UV divergences of the Wilson loops. We will have more to 
say about Wilson loops in section 3.3. 

The metric in all three regions can be written in the following form [62] 

1 



gidt^ + dx^ + dy^ + dz^ + Vh ^grrdr'^ + gmndx'^dx^ (2.91) 



with gi being the Black-Hole factors and h being the warp factor depends on all the 
internal coordinates (r, 6',, 0j, ip). To zeroth order in gsNf and gsM we have our usual 
relations: 

= ^, 9^'^ = 1 - ^, 9^^} = 1, 9^1dx-dx- = dsl.. (2.92) 

But in higher order in g^Nf, ggM, both the warp factor and the internal metric get 
modified because of the back-reactions from the seven-branes, three form fluxes and 
the localized sources we embed. We can write this as: 

h = ^ + \ grr = glr + 9rr^ 9mn = S'jnn + 5'mn (2.93) 
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where the superscripts denote the order of QsNf and QsM. Now observe that to hnear 
order in gsNf and QsM, gj/J = gl^^ = and we recover OKS metric (12.231) with warp 
factor given by fl2.3ip in Region 1. Thus only when we include higher order terms, 
the internal metric deforms from the metric of T^'^. Whereas for a given order in 
gsNf, gsM the warp factor h behaves logarithmically for small r and as l/r*" for large 
r. 

It is clear that one cannot jump from Region 1 to Region 3 abruptly. There should 
be an interpolating geometry where fluxes and the metric should have the necessary 
property of connecting the two solutions. This is Region 2 in our figure above. 

In the following, we briefly discuss the backgrounds for all the three regions . For 
a comprehensive analysis, please consult our paper 

2.4.1 Region 1: Fluxes, Metric and the Coupling Constants 
Flow 

As already mentioned, the metric of the entire geometry has the form given in fl2.9ip 
with warp factor h given by fl2.3ip in region 1. The internal space retains its resolved- 
deformed conifold form up to 0{gsNf). Beyond this order the internal space loses its 
simple form and becomes a complicated non-Kahler manifold. The background has 
a// the type IIB fluxes switched on, with the three- forms given by (12.39 p . (I2.40p while 
five-form and the axio-dilaton are given by (12. 41 p . 

2.4.2 Region 2: Interpolating Region and the Detailed Back- 
ground 

To attach a UV cap that allows conformal invariance i.e. vanishing beta function, 
we need at least a configuration of vanishing NS three-form and axio-dilaton that 
becomes constant as r — oo. This cannot be abruptly attached to Region 1: we need 
an interpolating region. This region, which we will call Region 2, should have the 
behavior that at the outermost boundary the three-forms vanish, while solving the 
equations of motion. The innermost boundary of Region 2 — that also forms the 
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outermost boundary of Region 1 — will be determined by the scale associated with 
the mass of the lightest quark, mo, in our system. In Fig 2.17, this is given by region 
in the local neighborhood of r^in = 'moTf^^ +rh, where Tq and are the string tension 
and the horizon radius respectively. 

For B2 to change its form from log(r) to behave as constant as r — )■ 00, we will 
need three form flux that behaves as 1/r for small r to run as 1/r", > 2 for large r. 
To get such behavior, it is useful to define two functions /(r) and M(r) as (see Fig 
2.18): 

o(r-ro) 

fir) ^ ^^^air-ro) ^ M{r) ^ M[l-f{r)], a»l (2.94) 

where M is as before related to the effective number of five-branes (or the RR three- 
form charge) and tq is the location of the sources, which we will elaborate shortly. 
Note that for r << Tq, f{r) ^ e''"'"", whereas for r > Tq, /(r) ^ 1. Thus for r smaller 
than the scale tq, f{r) is a very small quantity; whereas for r bigger than the scale 
ro, /(r) is identity. In terms of M(r) this means that for r < tq, M(r) ^ M whereas 
for r > Tq, M(r) 0. This will be useful below. 

Using these functions, we see that one way in which logarithmic behavior along 
the radial direction can go to inverse r behavior, is when the warp factor takes the 
following form: 



h 



E ^ (2-95) 



where Cj are constant numbers, and the denominator can be mapped to r(a) defined 
in f l2.56p with e(a) functions of ggNf.M.N and the resolution parameter a. L^s are 
functions of the angular coordinates {9i, (l)i,ip). For other details see [57]- The warp 
factor h has the required logarithmic behavior as long as the exponents of r are small 
and fractional numbers, and indeed switches to the inverse r behavior as soon as the 
exponents become integers. The question now is what are the background sources 
that give rise to a warp factor of the form (13.1961) . All the elements of OKS model 
only give rise to the logarithmic warp factor and it turns out, as we will discuss in 
some detail bellow, that we need to add sources (i.e. anti five branes) at the outermost 
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boundary of region 2 along with additional background fluxes. The specific point in 
the radial direction beyond which Region 3 starts tells us exactly where to add the 
sources and the AdS cap. 

The demarcation point can be found by looking at the behavior of H^s and i^/jR. 
For this we need to use the functions (12.941) to write the RR three-form. Our ansatz 
for F3 then is: 

F3 = (a. - ^-A^) E (sin de, Ad<f>,-Y.^ de, a 

^ 3gsM{r)Nfda , ^ . / ^ ^2 . ^ , , ^ Qa ^ ^1 • /, , 

A > : — dr A e^ A cot — sm 62 d(p2 — > cot — sm Ui dc 

2 ^ 47rr'M ^ \ 2 ^ ^ r^M 2 ^ 

- E ^^^f^'!^'"" sin e^ sin 6^ (cot ^ rf^^i + E 4^ T d^^) A (2 

where Oq = 1 + ^ and (c^, /i^) are constants. One may also notice three things: 
first, how the internal forms get deformed near the innermost boundary of the region, 
second, how the function /(r) appears for all the components, and finally, how Nf is 
not a constant but a delocalized functio The function /(r) becomes identity for 
r > Tq and therefore F3 — )■ for r > Tq. For r < Tq, the corrections coming from /(r) 
is exponentially small. Integrating F3 over the topologically non-trivial three-cycle: 

-ew, A I sin 9i dOi A rf^i - V ^ sin 62 d92 A #2 | (2.97) 
2 y r^(°) y 

we find that the number of units of RR flux vary in the following way with respect 
to the radial coordinate r: 

Mtot(r) = M(r) + — ^) V ^ (2.98) 



which is perfectly consistent with the RG flow, because for r < tq, and r — )■ re~3ilM^ 
Mtot decreases precisely as M — A^j as the correction factor e''~^'° coming from /(r) is 
negligible. For r > rg, Mtot shuts off completely. This also means that below tq, the 
total colors decrease by Mtot exactly as one would have expected for the RG flow 
with Nf flavors. Using similar deformed internal forms, one can also write down the 



13 



We will soon see that Nf in fact is the effective number of seven-branes. 
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--- For 0=0.5 
- - For 0=0.75 
— For o=l 
For o=1.25 



Figure 2.18: A plot of the /(r) function for tq = 5 in appropriate units, and various 
choices of a. Observe that for large a the function quickly approaches 1 for r > tq. 



ansatz for the NS three-form. The result is: 



^3 = E 



6gsM{r)K 



2tt 



(^cosec Y cosec 
27rr 2 



9s Nf 



dr 



(sin e, de, A #1 - E ^ sin 9, dO^ A ^0^) + E ^^'^lll^J^^ {-^H' \de^ 



r 



Afcot^^/^^-E- 
3gs dM{r 



% cot I de^] + g,^^ (h{r) cot | dO, + b,{r) cot | dO,^ 



Ae^ A dr + 



Att dr 



1 + gsNf 



1 ^ 9aVV/ 



sin^i dOi A rf^i A rfr - ^ [2 1^ + 9c/,A^/ - 



log I sm ^ sm — I + 63 (r) 

9a'gsNf' 



127r dr 



+ 



fozifr) dM(r) 

sin 6^2 ^6*2 A (i02 A dr - ^ ^ — sin 6^2 ^6*2 A (i02 A dr 



127r 



dr 



(2.99) 



with (ka, ■■■,qa) being constants and bn = J2m ""'^ where a„m = CLnmio? ■, g sN A and 
= €m{gsNf). The way we constructed the three-forms implies that H3 is closed. 



In fact the 0{df) terms that we added to f l2.99p ensures that. However F3 is not 
closed. We can use the non-closure of F3 to analyze sources that we need to add 
for consistency. These sources should in general be (p, q) five-branes, with (p, q) 



62 Chapter 2 The Gauge/ Gravity Correspondence 

negative, so that they could influence both the three-forms and since the Imaginary 
Self Duality (ISD) property of the three-forms is satisfied near r = r^i^ the sources 
should be close to the other boundary. A simplest choice could probably just be 
anti five-branes because adding anti D5-branes would change F3, and to preserve the 
ISD condition, would have to change accordingly. Furthermore, as we mentioned 
before, as r — )■ tq, both H3 = F3 ^ 0. Therefore r = tq is where Region 2 ends and 
Region 3 begins, and we can put the sources there. They could be oriented along the 
space time directions, located around the local neighborhood of r = tq and wrap the 
internal two-sphere (6*1,01) so that they are parallel to the seven-branes. However, 
putting in anti D5-branes near r = ro would imply non-trivial forces between the 
five-branes and seven-branes as well as five-branes themselves. Therefore if we keep, 
in general, the {p, q) five-branes close to say one of the seven-branes then they could 
get dissolved in the seven-brane as electric and magnetic gauge fluxes and F^^^ 

respectively. Here, as before, * is the hodge star operator and *F^^^ is the hodge 
dual F'^^\ In general a configuration of Dp branes is equivalent to a configuration of 
Dp + 2 brane with fluxes. This equivalence is due to M theory where a brane can be 
viewed as collapsed version of another higher dimensional brane. In type HA theory, 
this was first shown by [73] [71] where a DO brane is considered as a tubular D2 brane 
with electric and magnetic fields on the D2 brane world volume. Thus we may treat 
the five brane charges being soaked in by seven-branes, which in turn would mean 
that F3 in (12:961) and H3 in fICTD will satisfy the following EOMs: 

dFs = F(^) a A2{z) - d (Re r) A H3 
d*H3 = * a A2 (z) - d{Ci A F3) (2.100) 

where the tension of the seven-brane is absorbed in A2{z), which is the term that 
measures the delocalization of the seven branes (for localized seven branes this would 
be copies of the two-dimensional delta functions) and r is the axio-dilaton that we 
will determine below. In addition to that d * F^ will satisfy its usual EOM. 

However the above set of equations fl2.100p is still not the full story. Due to 
the anti GSO projections between anti-D5 and D7-brane, there should be tachyon 
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between them. It turns out that the tachyon can be removed (or made massless) by 
switching on additional electric and magnetic fluxes on D7 along, say, (r, ip) directions! 
This would at least kill the instability due to the tachyon, although susy may not be 
restored. For details on the precise mechanism on how brane-anti brane configurations 
are stabilized , one may refer to [73]- [TS]. But switching on gauge fiuxes on D7 would 
generate extra D5 charges and switching on gauge fiuxes on anti-D5s will generate 
extra D3 charges. This is one reason why we write {N,Nf,M) as effective charges. 
This way a stable system of anti-D5s and D7 could be constructed. 

To complete the rest of the story we need the axio-dilaton r and the five-form. 
The five- form is easy to determine from the warp factor h using fl2.4ip . The total 
five-form charge should have contribution from the gauge fiuxes also, which in turn 
would affect the warp factor. For regions close to r^i^ it is clear that r goes as z"^"'^! 
where z is the embedding fl2.16p . More generically and for the whole of Region 2, 
looking at the warp factor and the three-form fluxes, we expect the axio-dilaton to 

14 

go ao: 

r = [6o + fci/(r)]E^ (2.101) 

a ' 

where bi are constants and are functions of the internal coordinates and are com- 
plex. These Ca and the constants bi are determined from the dilaton equation of 
motion [55l [56] : 

r = — ^"^^ ' + p, q) sources 2.102 

«Im r \l~9 ^'^ 

where tilde denote the unwarped internal metric gmn-, and S'dt is the action for the 
delocalized seven branes. The /(r) term in the axio-dilaton come from the (p, g) 
sources that are absorbed as gauge fluxes on the seven-braneO. Because of this 



^^One may use this value of axio-dilaton and the three- form NS fluxes (|2.99p to determine the 
beta function from the relations ()2.21|) . To lowest order in gsNf we will reproduce the beta functions 
presented in section 2.2.1. Notice that for r > tq the beta function does not vanish and both the 
gauge groups flow at the same rate. This will be crucial for our discussion in the following subsection. 

^^The r"*^'") behavior stems from additional anti seven-branes that we need to add to the existing 
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behavior of axio-dilaton we don't expect the unwarped metric to remain Ricci-flat to 
the lowest order in QsNf. The Ricci tensor becomes: 

~ _ 2 d(rndn)T 9 ( fTll fm\ , 2 /f^(p,q)5-brane _ 1~ f^(p,g)5-braneA 

(2.103) 

where we see that Tlrr picks up terms proportional to e^^^ and derivatives of /(r), Nf{r), 
implying that to zeroth order in QsNj the interpolating region may not remain Ricci- 
flat. However since the coefficients are small, the deviation from Ricci- flatness is 
consequently small. 

Finally the warp factor can be obtained using the five-form equation of motion: 

d * dh'^ = ifa A Fs + kIo {f^'^^ A F^^^ - 7^ A 7^) A2{z) + 4^ tr F^^^A4{S) 

(2.104) 

where F^^'^ is the seven-brane gauge fields that we discussed earlier, F^^^ is the (p, q) 
five-brane gauge fields required for stabilization of five brane anti-five brane config- 
uration and they also give proper interpretation of the colors in the gauge theory 
sidcj. Here TZ is the pull-back of the Riemann two- form, and A4(5) is the term 
that measures the delocalization of the dissolved (p, q) five-branes over the space S 
embedded in the seven-brane (again for localized five-branes there would be copies 
of four-dimensional delta functions). The iJs A F-^ term in f l2.104p is proportional to 
^^^iJ^- This is precisely the form for the warp factor ansatz (13.1961) with the p{r) 
term there accounting for the M'^{r) term above. This way with the warp factor 
(I3.196P and the three-forms (I2.96P and (I2.99P we can satisfy (I2.104P by switching on 
small gauge fiuxes on the seven-branes and five-branes. 

Therefore combining (I3.196p . <^IM^ . (EHH]), ( I2.10ip and the five-form, we can 
pretty much determine the supergravity background for the interpolating region 



system to allow for the required UV behavior from the F-theory completion. The full picture will 
become clearer in the next sub-section when we analyze the system in Region 3. 

^^In fact one should view the gauge fluxes on the seven-branes and the five-branes as the total 
gauge fluxes that are needed to stabilize the system. We will see in the next subsection that the full 
stabilization would require additional fluxes, but the structure would remain the same. 
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'"min < r < tq. At the outermost boundary of Region 2 we therefore only have 
the metric and the axio-dilaton. Both the three-forms exponentially decay away fast, 
giving us a way to attach an AdS cap there. 

2.4.3 Region 3: Seven Branes, F-Theory and UV Comple- 
tions 

The interpolating region, Region 2, that we derived above can be interpreted alter- 
natively as the deformation of the neighboring geometry once we attach an AdS cap 
to the OKS-BH geometry. The OKS-BH geometry is the range < r < and 
the AdS cap is the range r > vq. The geometry in the range r^[j^ < r < vq is the 
deformation. Such deformations should be expected for all other UV caps advocated 
in [57]. In this section we will complete the rest of the picture by elucidating the 
background from r > tq in the AdS cap. But before that let us give a brief gauge 
theory interpretation of background. 

For the UV region r > ro we expect the dual gauge theory to be SU{N + M) x 
SU{N + M) with fundamental flavors coming from the seven-branes. This is because 
addition of {p, q) branes at the junction, or more appropriately anti five-branes at the 
junction with gauge fluxes on its world-volume, tell us that the number of three-branes 
degrees of freedom are -|- M, with the M factor coming from five-branes anti-five- 
branes pairs. As mentioned in the previous section, one comes to this conclusion 
using M theory, as five anti-five brane configuration with flux is equivalent to three 
brane charge. Furthermore, the dual gauge theory for AdS conifold geometry is the 
Klenabov-Witten theory with gauge group SU{N) x SU{N) for some A^. Thus the 
gauge theory dual to region 3 is indeed SU {N + M) x SU {N + M), but has RG flows 
because of the fundamental flavors (This RG flow is the remnant of the flow that we 
saw in the previous subsection. We will determine this in more details below). 

At the scale r = ro we expect one of the gauge groups to be Higgsed, so that we 
are left with SU{N + M) x SU{N). This Higgsing is justified as follows: the five 
and anti-five branes are separated by distance ro, which means the gauge bosons that 
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arise due to strings stretching between five and anti-five branes have length ~ tq. 
Thus these bosons are massive with mass O{ro = Aq). For scales less than ro = Aq, 
these bosons are not produced and we only have the gauge theory arising from D3, D5 
and seven branes i.e. OKS-BH type theory with gauge group SU{N + M) x SU{N). 
Now both the gauge fields fiow at different rates and give rise to the cascade that is 
slowed down by the Nf flavors. In the end, at far IR, we expect conflnement at zero 
temperature. 

The few tests that we did above, namely, (a) the flow of and M colors, (b) the 
RG flows, (c) the decay of the three-forms, and (d) the behavior of the dual gravity 
background, all point to the gauge theory interpretation that we gave above. What 
we haven't been able to demonstrate in is the precise Higgsing that takes us from 
Klebanov-Witten type gauge theory to Klebanov-Strassler type cascading picture. 
From the gravity side its clear how this could be interpreted. From the gauge theory 
side it would be interesting to demonstrate this. 

Coming back to the analysis of region 3, we see that in the region r > tq we do 
not expect three-forms but we do expect non-zero axio-dilaton. These non-zero axio- 
dilaton come from the the seven branes that are present in region 3. Of course the 
complete set of seven-branes should be determined from the F-theory picture [80] [81] 
to capture the full non-perturbative corrections. This is now subtle because the 
seven-branes are embedded non-trivially here (see fl2.16p ). A two-dimensional base, 
parametrized by a complex coordinate z, on which we can have a torus flbration: 

= x^ + xF{z)+G{z) (2.105) 

can be identifled with the z coordinate of (I2.16p . This way vanishing discriminant 
A of fl2.105p i.e A = AF^ + 276*^ = 0, will specify the positions of the seven-branes 
exactly as f l2.16p . Here we have taken F{z) as a degree eight polynomial in z and 
G{z) as a degree 12 polynomial in z. 

As is well known, embedding of seven-branes in F-theory also tells us that we can 
have SL{2, Z) jumps of the axio-dilaton [58] [H]. We can deflne the axio-dilaton r = 
Co -|- ie~'^ as the modular parameter of a torus flbered over the base parametrized 
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by the coordinate z. The holomorphic maj^ from the fundamental domain of the 
torus to the complex plane is given by the famous j-function: 

■^^ ' ~ r]^^{T) 27G^{z) + AF-\z) ^ ' ' 

where Qi,i = 1,2,3 are the well known Jacobi Theta-functions and rj is the Dedekind 
?7-function: 

V{r) = g^LKl-^'^)' <1 = e'^'^ (2-107) 

n 

For our purpose, we can write the discriminant A{z) and the polynomial F{z) gener- 
ically as: 

24 8 

A{z) = AF^ + 27G^ = aY[{z-Zj), F{z) = bl[{z - Zi) (2.108) 

j=i i=i 

where zj are (like n in (12.161) ) complex constants and the jth seven brane is located 
at z = Zj. When we have weak type IIB coupling i.e t = Gq + ioo, j{T) ^ g-2T«T 
using fl2.106p the modular parameter can be mapped to the embedding coordinate z 
as: 

! 24 



i i i °° 

- + — log (55926ari) - — E 
9s 2tt Its: 



1 



nz" 



= E tj^" (2.109) 

n=0 ^ 

where C„ = CniOi, (pi,^) and Vn = VniOi, (l)i,ip) are real functions and f = r^/^. To 
avoid cluttering of formulae, we will use r instead of r henceforth in this section 
unless mentioned otherwise. So the coordinate r will parameterize Region 3, and 

The above computation was done assuming that z > {zi,Zj), which at this stage 
can be guaranteed if we take 61^2 small. This gives rise to special set of configurations 
of seven-branes where they are distributed along other angular directions. However 
one might get a little worried if there exists some Zj = Zq related to the farthest seven- 
brane(s) where the above approximation fails to hold. This can potentially happen 



^^Holomorphic in r, the modular parameter. 
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when we try to compute the mass of the heaviest quark in our theory. The question 
is whether we can still use the r derived in (12. 1091) . or we need to modify the whole 
picture. 

Before we go into answering this question, the choice of z bigger than {zi,Zj) 
already needs more convincing elaboration because allowing 612 small is a rather 
naive argument. The situation at hand is more subtle than that and, as we will argue 
below, the picture that we have right now is incomplete. 

To get the full picture, observe first that z being defined by equation fl2.16p means 
that if we want fl2.109p to hold in Region 3, we need to specify the condition r > ro in 
f l2.16p . This way a given z will always imply points in Region 3 for varying choices of 
the angular coordinates {6i,(f)i,ilj). However a particular choice of {zi,Zj) may imply 
very large r with small angular choices or small r with large angular choices. Thus 
analyzing the system only in terms of the r coordinate is tricky. In terms of the full 
complex coordinates, z > {zi, Zj) would mean that we are always looking at points 
away from the surfaces given hj z = Zi and z = Zj. 

What happens when we touch the z = Zi surfaces? For these cases F{zi) — )■ 
and therefore we are no longer in the weak coupling regime. For all F{zi) = imply 
j(r) — > which in turn means r = exp (ivr/3) on these surfaces. These are the 
constant coupling regimes of [H2]-[B1] where the string couplings on these surfaces are 
not weak. On the other hand, near any one of the seven-branes z = Zj we are in the 
weak coupling regimes and so ( I2.109P will imply 

t{z) = — log {z - Zj) zoo (2.110) 

which of course is expected but nevertheless problematic for us. This is because we 
need logarithmic behavior of axio-dilaton in Region 2, but not in Region 3. For a 
good UV behavior, we need axio-dilaton to behave like f l2.109p everywhere in Region 
3. 

In addition to that there is also the issue of the heaviest quarks creating additional 
log divergences that we mentioned earlier. These seven branes are located at z = zj = 
Zo, and therefore if we can make the axio-dilaton independent of the coordinates Zo 
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then at least we won't get any divergences from these seven-branes. It turns out 
that there are configurations (or rearrangements) of seven-brane(s) that allow us to 
do exactly that. To see one such configuration, let us define F{z),G{z) and A{z) in 
(12.1081) in the following way: 

7 10 

F{z) = (z~Zo)l[{z-Zi), G{z) = {z-Zofl[{z-z,) 

i=l i=l 
21 

A(^) = (z-zoflliz-z,) (2.111) 

i=i 

which means that we are stacking a bunch of three seven-branes at the point z = Zo, 
and 

21 7 10 

Uiz-Zj) = Al[{z-Zif + 27{z-Zo)Y[{z-z,f (2.112) 

j=l 1=1 i=l 

implying that the axio-dilaton r becomes independent of z^ and behaves exactly as 
in (12.1091) with in (12.1091) varying up to (7, 21) respectively. 

The situation is now getting better. We have managed to control a subset of 
log divergences. To get rid of the other set of log divergences that appear on the 
remaining twenty-one surfaces, one possible way would be to modify the embedding 
(I2.16p . In fact a change in the embedding equation will also explain the axio-dilaton 
choice (12.1011) of Region 2. To change the embedding equation (I2.16P we will use 
similar trick that we used to kill off the three-form fluxes, namely, attach anti-branes. 
These anti seven-branea^^l are embedded via the following equation: 

^3/2gi(^-<^i-02)gi^ ^ sin y = roe'® (2.113) 

where is some angular parameter, and could vary for different anti seven-branes. 
The above embedding will imply that their overlaps with the corresponding seven- 



branes are only partial 







And since we require 



^^They involve both local and non-local anti seven-branes. 

""^^For example if we have a seven-brane at z = zi such that lowest point of the seven brane is 
r = < To, then the corresponding anti-brane has only partial overlap with this. 
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it will appear effectively that we can only have seven-branes in Regions 1 and 2, 
and bound states of seven-branes and anti seven-branes in Region 3o This way the 
axio-dilaton in Region 3 will indeed behave as (I2.109p . with the seven values of Zi in 
(I2.112P chosen to be in region 1 and 2. 

There are two loose ends that we need to tie up to complete this side of the story. 
The first one is the issue of Gauss' law, or more appropriately, charge conservation. 
The original configuration of 24 seven branes had zero global charge, but now with the 
addition of anti seven-branes charge conservation seems to be problematic. There are 
a few ways to resolve this issue. First, we can assume that that branes wrap topologi- 
cally trivial cycles, much like the ones of [38] . Then charge conservation is automatic. 
The second alternative is to isolate six seven-branes using some appropriate F and G 
functions, so that they are charge neutral. This is of course one part of the constant 
coupling scenario of [5S]. Now if we make the (6'2,02) directions non-compact then 
we can put in a configuration of 18 seven-branes and anti seven-branes pairs together 
using the embeddings ( 12.16P and ( 12.113P respectively. The system would look effec- 
tively like what we discussed above. Since the whole system is now charge neutral, 
compactification shouldn't be an issue here. 

The second loose end is the issue of tachyons between the seven-brane and anti 
seven-brane pairs. Again, as for the anti-D5 branes and D7-brane case [75]-|78]. 
switching on appropriate electric and magnetic fluxes will make the tachyon massless! 
Therefore the system will be stable and would behave exactly as we wanted, namely, 
the axio-dilaton will not have the log divergences over any slices in Region 3. 

This behavior of axio-dilaton justifies the r~'^(") in fl2.10ip in Region 2. So the full 
picture would be a set of seven-branes with electric and magnetic fluxes embedded 
via fl2.16p and another set of anti seven-branes embedded via (12.1130 lying completely 
in Region 3. 

Thus in Region 3 both the three-forms vanish and therefore gi = g2 = Oym with 



^"^Of course this effective description is only in terms of tlie axio-dilaton charges. In terms of the 
embedding equation for the seven-branes (|2.16l) this would imply that we can define z with r > Tq 
and ij with r < Tq. 
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91,92 being the couplings for SU{N + M), SU{N + M). From (12.211) we can compute 
the /9-function for gyu as: 

^^^^^^ = SbiA - rn^h^ ^ ^ ^ 

where A is the usual RG scale related to the radial coordinate in the supergravity 
approximation. For A — )■ oo, /3((7ym) implying a conformal theory in the far 
UV. We can fix the 't Hooft coupling to be strong to allow for the supergravity 
approximation to hold consistently at least for all points away from the z = Zi^i = 
1, 7 surfaces. 

Existence of axio-dilaton r of the form (I2.109P and the seven-brane sources will tell 
us, from (12.1031) . that the unwarped metric may not remain Ricci flat. For example 
it is easy to see that 

^ ^0 n,m=l ^ ^ 

where the last term should come from the seven-brane sources and, because of these 
sources, we don't expect TZrr to vanish to lowest order in 9sNfCJ The term Av is 
given by the following infinite series: 



^^ = i-E^;ir^+ E .Up+V +■■• (2.116) 

k,l=l k,l,p,q=l 



Similarly one can show that 

n T (^'^'^n + idaVn){dbCm - idbV^) I ^ / 1 \ .2 117) 

^0 n,m=0 ^ ^ 

for (a, h) 7^ r. For TZrb similar inverse r dependence can be worked out. In the far UV 
we expect the unwarped curvatures should be equal to the AdS curvatures. The warp 
factor h on the other hand can be determined from the following variant of (I2.104p : 

d * dh-^ = tr (f(^) a F(^) - 7^ a 7^) A2{z) + ... (2.118) 

because we expect no non-zero three- forms in Region 3. The dotted terms are the 
non-abelian corrections from the seven-branes. As r is increased i.e r >> tq, we 



^^Although, as discussed before, the deviation from Ricci flatness will be very small. 
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expect F^^'^ to fall-off (recall that they appear from the anti (1,1) five-branes located 
in the neighborhood of r = tq) and therefore can be absorbed in TZ. Once we embed 
the seven-brane gauge connection in some part of spin- connect ion, we expect 



°h-' = ^[-n) (2-119) 

where the box represents combinations of differential operators that arise from fl2.118p . 
Solving this will reproduce the generic form for h: 



4 



J.4 



1 + E 



(2.120) 



i=l 

with a constant and a/s are suppressed by powers of gsNf. More details on this 
is given in the Appendix A and B of [62j. At far UV we recover the AdS picture 
implying a strongly coupled conformal behavior in the dual gauge theory. 

To summarize, we have obtained the dual gravity of a thermal field theory with 
matter in the fundamental representation . The gauge theory becomes almost confor- 
mal in the UV with massive Higgs-like gauge bosons with matter transforming under 
SU {N + M) X SU {N + M) gauge symmetry. In the IR the massive gauge bosons 
are Higgsed away and we end up with a SU {N + M) x SU{N) gauge theory. The 
dual geometry has the metric of the form fl2.9ip with a warp factor given by fl2.3ip for 
r <^ To, by fl3.196p for r ~ rg, and of the form f l2.120p for r ^ tq. Thus in UV we have 
Klebanov-Witten type field theory while in the IR we have Klebanov-Strassler model 
with fundamental matter and temperature. Once we reach the Klebanov-Strassler 
type description of our field theory, we expect Seiberg duality cascade to occur as 
we go further down in energy scale. At the bottom of the cascade we end up with 
strongly coupled SU (M) gauge theory with Nf fundamental flavors and logarithmic 
running of coupling. The scenario is sketched in Fig 2.19. 
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SU(N+M) XSU(N+M) 
g. 1/A" 

Asymptotic AdS geometry 



Flow from UV to IR 



SU(N+M) X SU(N) 
g~log A 

Klebanov-Strassler type geometry 



Flow from UV to IR 



SU(M) 
g~log A 

Klebanov-Strassler type geometry 



Figure 2.19: Schematic depiction of flow in the gauge theory 



Chapter 3 

Application to Thermal QCD 



Having discussed in some detail the gauge theories that arise from various brane setups 
and their dual geometries, now we attempt to make connections with thermal QCD 
in the large N limit. The gauge theory arising from the brane setup in section 2.4 is 
obtained by UV completing OKS-BH geometry and is asymptotically conformal with 
logarithmic running of couplings in the IR. To our knowledge it is the brane setup 
that resembles large N QCD the most in the 'top down' approach where gauge/gravity 
duality is exactly derivable. Unlike the 'bottom up' approach where we do not know 
the exact brane configuration in higher dimensions or the precise origin of the gauge 
theory, our approach is based on open/closed string duality and in some sense the 
correspondence is more theoretically complete. 

As sketched in Fig 2.19, although the gauge theory has a very rich structure, it is of 
course not exactly QCD. However, the benefit of this construction is at large coupling 
it has an exact classical supergravity description. Certain gauge theory observables 
which are extremely difficult to obtain using conventional field theory techniques 
become rather simple to calculate using dual gravity. The field theory lives in four 
flat spatial dimensions and incorporates matter in fundamental representation. We 
will refer to the fundamental matter as 'quarks' of our theory. The theory has no UV 
Landau poles while effective degrees of freedom converge in the far UV. As the theory 
reaches conformal fixed point at large energy scales, our construction maybe viewed 
as the UV completion of Klebanov-Strassler model which is most relevant for QCD. 
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Although the field content is somewhat different from large N QCD, our theory 
has several common features to it and in fact becomes almost large N QCD in the 
IR. Thus the physics extracted from our gauge theory could be very similar to that 
of QCD and it is worth analyzing the thermodynamic properties of the field theory 
plasma. As this may turn out to be crucial in understanding the collective excitations 
of QCD plasma, in this chapter we will examine the dynamics of 'quarks', transport 
coefficients of the medium and the confining nature of the our gauge theory. 



3.1 Quark Dynamics 

Strings in ten dimensional dual gravity with endpoint on a D7 brane corresponds to 
fundamental matter i.e. quarks in four dimensional Minkowski space time. As the 
D7 brane fills the four dimensional Minkowski space, the endpoint of the string is a 
point there and the quark is localized at the endpoint. Energy of the string gives 
energy of the quark and minimal energy of the static string configuration in the bulk 
geometry gives the mass of the quark. If we consider a string traveling through the 
bulk geometry, it corresponds to a quark moving through gauge theory plasma. In 
the following subsections we will use this key concept to compute thermal mass of a 
quark and the drag it experiences as it traverses a medium. We will also quantify the 
wake it leaves behind in the medium along with transverse momentum broadening of 
a fast moving quark. 



3.1.1 Thermal Mass and Drag of a Quark 

We start with the action of a string in ten dimensional geometry with metric ( I2.23p . If 
X*((T, r) is a map from world sheet coordinates a, t to 10 dimensional space time, then 
string action or fundamental string Born Infeld action is (see for example [85] |86j): 
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where J{(f)) is the additional couphng of the dilaton (f) to the string world-sheet, Tq is 
the string tension, X"' are the ten bosonic coordinates, G is a 32 component spinor, 
Fabc... = [dC]abc.. with Cabc. being the background RR form potentials, and f^s is 
world sheet metric, given by the standard pull-back of the space time metric on the 
world-sheet: 



where we have taken the background metric of the form f l2.9ip ignoring ggNj terms 
which means Qmn is the metric of T^'^ in fl2.9ip . Here 7, 5 = 0, 1 with parametrization 
rp = T = t and r]^ = a = r. 

In the ensuing analysis we will keep Bab-, J{<P) £^s well as (9^0 zero and the jus- 
tification will be given shortly. The interesting thing however is to do with the 
background RR forms. Note that the RR forms always couple to the 32 component 
spinor. Therefore once we switch-off the fermionic parts in (13. ip . the fundamental 
string is completely unaffected by the background RR formi^. Thus in the following 
analysis, for the mass and drag of the quark, we can safely ignore the RR fields. We 
have also defined: 



where Gij is more generic than the background metric, and could involve the back 
reaction of the fundamental string on the geometry. The analysis is very similar to 
the AdS case discussed in |87][88], however, since our background geometry in [57] 
involves running couplings, the results will differ from the ones of [BZ] [HB] ■ 

As mentioned earlier, a fundamental quark will be a string starting from the D7 
brane and ending on a D3 or fractional D3 (which is a D5 brane) brane, giving the 
quark color. At non zero temperature, the dual geometry has a black hole and quark 

^This is of course the familiar statement that the RR fields do not couple in a simple way to the 
fundamental string. 





det / = ~G,,PX'^ + (G,,X'*X'^)(GhX'=X') 



(3.1) 
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in the thermal medium is represented by a string stretching between seven brane and 
the horizon rh of the black hole. For simplicity of the calculation, we will then restrict 
to the case when 

X° = t, X^ = r, X^ = x(a,r), X'^ = (A; = 2, 3, 5, 6, 7) 

(x^ X') = {6,, = TT, e = e = (3.2) 

and we choose parametrization t = t,a = r also known as the static gauge. Thus 
we are only considering the case when the string extends in the r direction, does 
not interact with the RR fields, and moves in the x direction of our manifold. More 
general string profile, while being computationally challenging, does not introduce 
any new physics and hence our simplification is a reasonable one. 

Before moving further, let us consider two points. First is the effect of the black 
hole on the shape of the D7 brane. We expect due to gravitational effects the D7 
brane will sag towards the black hole and eventually the string would come very 
close to the horizon. In fact putting a point charge on the D-brane tends to create 
a long thin tube on the D-brane that in general extends to infinity. The end point 
of the string being a source of point charge should show similar effects (see [90] for a 
discussion of a somewhat similar scenario)!^. For our analysis here we will ignore this 
effect altogether and we hope to address the issue in our future work. 

The second point is to see how the background varying dilaton and NS-NS two 
form effect the string. With the axion and dilaton behaving as f l2.4ip near the location 
of D7 brane and as 0{gsNf /r"") globally, da(f) ~ 0{gsN f /r"^) for all r at the location 
of the string. Note that QsNf ^ 1 and as the string extends from r = tqIo r = rh., for 
r/i large i.e. high temperature, 0{gsNf/r^) is negligible at the location of the string. 
This means we can ignore the contribution from 9^0 in the string action. Similarly 
we can ignore effects of B2 as it gives 0{gsM\ogr) for small r and 0{gsM/r"^) for 
large r, both of which can be ignored at the location of the string. This is because 

^Even in the supersymmetric case, putting a point charge on a D-brane tends to create a long 
tube that extends to infinity. One can then view an open string to lie at the end of the thin tube. 
This effect is somewhat similar to the one discussed in [89l . 
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for reasonably large temperature, the IR divergence of log(r) is absent (temperature 
which is of 0{rh) is the IR cutoff) and at the UV we can ignore l/r" terms. On the 
other hand, the dilaton additionally couples to the string world sheet through J(0) 
term which is proportional to the Ricci scalar i?(2) of the world sheet metric, f.ys and 
under a reparametrization of the world sheet metric, we can make i?(2) small enough 
as there are no divergences in the background metric. This way we can ignore 
near the location of the string. However in section 3.3 we will consider the effects of 
i?2 and dilaton on the string world sheet and we will see that these effects although 
rather technical to track, do not change the overall physics. Thus we can simply set 
J(0) = B2 = df^(j) = in the following analysis. 

Getting back to the string action, with our choice of parametrization and string 
profile we get: 

-det/=»iM + -»i(rL."-9.W-'i^ (3.3) 

g2{r) h{r, 71,71) 

where the warp factor h{r, 61, 62) = h{r, vr, tt) , evaluated at the location of the string. 

With this, the rest of the analysis is a straightforward extension of [HZ|[HH]- The 

Euler-Lagrangian equation for = x{t,r) derived from the action (13. ip and the 

associated canonical momenta are: 

1 / X \ d / gix' \ ^ 
52C^t V-det c/rUv-det 

n, - -ToG., ^=== (3.4) 

If we consider a static string configuration, i.e. x(cr, r) = b = constant, then energy 
can be interpreted as the thermal mass of the quark in the dual gauge theory. Using 
the static solution in (13. 4p . we obtain the thermal mass m(T) using E = — J da 11°, 

as 

m(r) = Toiro - r,) = Tq (I/x^/^ - t) (3.5) 



Observe that the thermal mass decreases as temperature is increased. In our 
analysis, /x^/^ = ro > T and /x^/^ is proportional to zero temperature mass. This 
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means, we are dealing with heavy quarks where the temperature is always less than 
the mass of the quark. On the other hand, our predictions for the thermal mass gets 
exact in the limit where number of colors is infinite. Thus our results suggest that at 
large N and small temperatures, the heaviest quarks get less massive as temperature 
is increased. 

When a probe particle moves through a plasma, it interacts with the medium 
through collisions with the constituents of the medium and if it is charged, it also 
radiates. Overall because of the interaction with the medium, the probe experiences 
drag force and this drag is a key characteristic of the plasma. As a moving string 
in the bulk corresponds to a moving quark in the medium, we can compute the 
drag coefficient by considering a string moving with some velocity with the endpoint 
representing a quark which traverses the medium. For the computation of drag, 
that is to analyze the response of the medium to a moving probe, it is enough to 
consider a constant velocity quark and calculate the drag it experiences. Of course 
drag will try to slow down the quark and we need to apply force to keep it at constant 
speed. If there is no external force applied, the probe will slow down but the drag 
coefficient we extract considering constant speed remains unchanged - as we can 
consider instantaneous velocity to be the constant velocity and the analysis stays the 
same. 

We consider the string profile fl3.2p with 

x(t,r) = x{r) + vt (3.6) 

Then from (13. 4p . noting that / is independent of time, we can solve the equation 
of motion to get: 

9i92 91 - hC^v^ 

where C is a constant of integration that can be determined by demanding that 
— det / is always positive. Using the value of x'^ from (13.71) we can give an explicit 
expression for the determinant of / as: 
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For — det / to remain positive for all r, we need both numerator and denominator 
to change sign at same value of r. This is the same argument as in [87] [88]. The 
numerator changes sign ajf] 



+ 0{g,Nf,gsM) 



(3.9) 



where we use h near horizon to be of the form in (12.311) . □ Requiring that denominator 
also change sign at that value fixes C to be: 

1 



C 



1 + log 



(1-1-2)1/4 



(1 + 2^ {log 



(l-i;2)i/4 



(3.10) 



where M and Nf differs from M,Nf due to the 0{gsNf, QsM) terms in (13. 9p . The 
first part of C is the one derived in [87] [88]. The next part is new. Now the rate 
at which momentum is lost to the black hole is given by the momentum density at 
horizon 



(3.11) 



while the force quark experiences due to friction with the plasma is ^ = —vp with 
p = mv/\^l — v"^. To keep the quark moving at constant velocity, an external field 
£i does work and the equivalent energy is dumped into the medium [87] [88]. Thus 
the rate at which a quark dumps energy and momentum into the thermal medium is 
precisely the rate at which the string loses energy and momentum to the black hole. 
Thus up to 0{gsNf,gsM) we have = -Uf{r = rh) and 



TqCVT^ 



m 



(3.12) 



Tn 



1 + log 



r 



(l_^2)l/4 



(1 + 2^ {log 



r 



(1-1,2)1/4 



^Note that by 0{gsNf,gsM) we will always mean 0{gsN f , glM N f , QsM"^ / N) unless mentioned 
otherwise. 

"^This is justified as for all the geometries considered in section 2.2-2.4, the near horizon r ^ rh 
warp factor is always of the form (|2.3ip . We have various choices for the large r behavior of the 
geometry and subsequently for h, but the IR behavior remains the same as the OKS model. 
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which should now be compared with the AdS result [87] [88]. In the AdS case the 
drag coefficient u is proportional to T^. For our case, when we incorporate RG flow 
in the gravity dual, we obtain O A log T + B log^T^ correction to the drag 
coefficient computed using AdS/CFT correspondence [87] [88] . 

3.1.2 Transverse Momentum Broadening 

With the computation of thermal mass of the quark and the drag it experiences in 
the medium, we will now study the diffusion process through which a probe parti- 
cle transfers momentum with the plasma. The analysis is of particular interest as 
it relates to the formation of the quark gluon plasma at the relativistic heavy ion 
colliders. To be more precise, at the earliest stages of a central collision, energy den- 
sities are expected to be high enough to form the quark gluon plasma and several 
observables have been proposed to probe the plasma. Strangeness enhancement [91], 
J/ip suppression [92], electromagnetic radiation [93] are among the key candidates 
while the ones that are produced at the earliest stages of the collision gather special 
attention as they directly interact with the plasma. Hard scatterings where the par- 
tons transfer large momentum take place right after the collision as only then there 
is enough energy and the resulting partons fragment into jets of hadrons with large 
transverse momentum Pt- The partons produced in the heavy ion collisions are ex- 
pected to lose energy in the medium [91] which should result in a suppression of high 
Pt hadrons [95j when compared to the high Pt hadrons produced in proton-proton 
collision. Experiments at RHIC have indeed observed this suppression [96]- [98] and 
this phenomenon is known as 'jet quenching'. 

There has been a lot of effort to model the energy loss mechanism [99]- |109j 
and gluon bremsstrahlung with Landau-Pomeranchuk-Migdal (LPM) effect |110] is 
thought to be the dominant process in jet quenching. As the medium formed in the 
coUison expands, one needs to model the time evolution of the system accounting the 
dynamics of the fluid and then compare with experimental data of suppression. As 
the dual geometry we constructed is time independent, we cannot address the effects 
due to the evolution of medium. Rather we will quantify the diffusion process by 
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computing the mean square transverse momentum transfer between the non expand- 
ing medium and a fast moving parton. In principle, one needs to construct a time 
dependent dual geometry by considering the collision of strings ending on D7 branes 
and computing their back reaction, similar in the spirit of Then using this 

time dependent background geometry to compute Wilson loops, one can calculate 
momentum distribution and finally quantify 'jet quenching'. Constructing the dual 
geometry of a heavy ion collision is rather challenging and we hope to address the 
issue in our future work. 

For our current analysis, consider a parton moving through a plasma in four 
dimensional Minkowski space time with the following world line 



x{t) = vt 

z{t) = y{t)=5y{t) (3.13) 

For a fast moving parton, we can always choose coordinates such that f l3.13p is 
the world line. Now if the plasma has matter in fundamental representation, has 
logarithmic running of coupling in the IR but becomes asymptotically conformal, we 
can treat the ten dimensional geometry with metric of the form fl2.24p to be the dual 
gravity of this gauge theory which lives in four dimensional flat space time. 

To obtain the momentum broadening of the parton we shall use the Wigner dis- 
tribution function / as defined in QCD kinetic theory |112] 

/(X,r^) = </,,(X,rx)>=Tr[pQl(X^)Uab(X^,Xi)Qa(Xi)] (3.14) 

where = X — r_|_/2, X^ = X + r^jl^ X = (t,x) is the world line of the parton 
field Qa without any fluctuation, Uab is the link and p is the density matrix. Of 
course the index ab refers to color and for our choice of the world line for the parton 
we have r'j_ = 25y^. Now Fourier transforming the distribution functions, we can get 
the average transverse momentum to be 

<pI>=J d'x J ^,klf{X, k^) (3.15) 
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X=(to,Xo,y(l,)) 

X=(t,x +vAt,y(t)) 

V J 

Xf=(to-is,Xo-is,y(to-is)) 

Figure 3.1: The contour C with X denoting the coordinates of a point on the contour. 
The upper hne with real time coordinate is the world line of the heavy parton while 
the lower line has complex time coordinate. Here s — )■ is real. 



Now if we assume that the initial transverse momentum distribution is narrow, 
then we have 

<pI> = 2iiTT 

(3.16) 

with T some large time interval, and is the diffusion coefficients. 

Following the arguments in |114j . one can write the diffusion coefficients solely in 
terms of functional derivative of Wilson loops. The final result is 

KT = lim - / dt e*'^* (iGUt, 0) + iGUt, 0) + iGUt, 0) + iGUt, 0)) (3.17) 
4 J 

where the Greens functions are 

Gy(tt') - 1 /tr.° ^'^^[^y^'°^ 

^^^ '^ " trpOWc[0,0]^''^ 5y,(t)5y,(tO^ 
^^^ '^ " trpOWc[0,0]^*'^ 5yi(t)5y2(tO ^ 

trp'JWc[0,OJ 5yi(t)5y2(t') 

where t,t' are the real part of complex time tct'^ on contour C (Fig 3.1) . Here 
we have introduced type '1' and '2' fields {Syi,i = 1,2) of thermal field theory in 
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real time formalism |115] |116] evaluated with real and complex time coordinates - 
hence evaluated at the upper and lower horizontal line of contour C . We denote by 
Wc[Syi,Sy2] the Wilson loop with deformation 6yi and 6y2 on the upper and lower 
line of C. 

We will now compute the Wilson loop at strong coupling by using holography, 
that is we identify 

where p° is the density matrix |114] . S'ng being the Nambu-Goto action, with the 
boundary of the string world sheet being the curve C. That is the string world sheet 
ends on the world line fl3.13p of the heavy parton. If : (a, r) — )■ (t, (, x, y, z, ip, 0i, 02, Oi, 
is a mapping from string world sheet to ten-dimensional geometry given by f l2.23p with 
( = 1/r, then with parametrization a = (,t = t, we have 

x{tX) = vt + xiC) 
z{t,0 = y{t,0 = Sy{t,0 (3.20) 

where x is the unperturbed solution for the mapping. Now using background metric 
of the form fl2.9ip with grr = H, the Nambu-Goto action up to quadratic order in the 
perturbation 6y gives 



where we approximated Qi = g and prime means derivative with respect to ( while 
dot means derivative with respect to 't', A = ^,A' = Minimizing this action 
with respect to x with ignoring 0{6y) gives the solution 

~ 9K\9 - ^'■''> 
where C is a constant of integration. 

Now to solve for the transverse fluctuation 6y, note that the puUback metric 
/a/3 = daX^dpX^G^y is off diagonal with components 



7 
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with 5y = d5y/dt, 5y' = d5y/dC, g{a) = 1 — for any a and 7 — l/\/T^^v^ . This 
puUback metric can be diagonahzed at zeroth order in y with the reparametrization 



C 

dF 



The resulting puUback metric is 



vx' 



9(0 - 



(3.24) 



fti 

k 



= 0{Sy'); 



-9(0 + 



+ 



1 



-./2 



7VC%(C)-CW(C) ^/,(c) 

Sy\M(?C'H{0 



(3.25) 



5(C)^VMC)C^(5(C)-t'^)(^(C)-C^VMC))/ 

Here 5y = ■yJl^V, — ^ = Observe that the puUback metric has 

a horizon at C = Cft which means at C = Ch/ ^/l < Ch and thus the puUback has 
larger horizon radius than the space-time metric. Thus the world sheet only extends 
in the radial direction from ( — to ( — Ch/-\/l- Also observe that at the horizon 
g — C'^v^h = by our choice of C and thus indeed /^^ = 00 at the horizon. Now 
with this reparametrization , the terms quadratic order in 5y in Nambu-Goto action 
become 



-[2] 

^NG,5y 



A 



27ra' J 

g(C)g(C) 

g-Wh{Q 



9{C)C' 



2^''h{C)VT+A 
- 1 



1 + -f^hipv^C^ 
2^9(0 



-Sy 



(3.26) 
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If we denote Sy{t,() = J du e*'^* Sy{u})Y{(), then the equation of motion can be 
written as 



D^-^!il±iMO^ (3.27) 



We try solution to fl3.27p of the form 

Y = g{lf7 

J^={l + f3n) (3.28) 

with /3 being a constant and we have written J-" only up to linear order in /3 which is 
sufficient for what is to follow. Now from equation (I3.27P near horizon C, ^ C,h and 
only considering up to quadratic order in /3, we can isolate the most divergent terms 
to obtain 



= (3.29) 

where ^ = ^(C = Ch), S is of 0{gsM''^^/N,gl{Nf)f and h = h{C = ChI ^) while 
T = ^^'(^(^^-^ is the temperature associated with the space time black-hole (not the 
world sheet black hole). This determines our constant /3 and now we can solve fl3.27p 
order by order in Cj. For the purpose of calculating kt-, we will eventually take the 
zero frequency — )■ limit (albeit after dividing by w), thus it is sufficient to solve 
(I3.27P only up to linear order in u. With the form of the solution as in fl3.28p . the 
equation of motion fl3.27p gives the following equation for 



(3.30) 
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above form of the equation, we observe with ly = "H' we can write it as 

'gjcr ~h'9{cy \ 

. 9(0 hgiO J 



where h = 2-f^h{C)^^^ and prime denotes a derivative with respect to (. From the 



[WB]' = -B 



= B(^M-'m-] (3.31) 

We will impose the boundary condition 1^(0) = 1 which implies 'H(O) = 0. Now 
using the solution for Y and taking appropriate linear combinations to build the 
type '1' and '2' fields Syi,6y2 as in |114j . we can write the boundary action after 
integrating the Nambu-Goto action and the result is equation (3.51) of |114] but Y 
is replaced with our solution, u,u replaced by cj/ttT, w/vrT and R= 1. Finally from 
the boundary action we can obtain the Greens function Gij and the result for the 
diffusion coefficient is 



f^T = V7^7siVcfrvrT^(l + B) (3.32) 

where B is of 0{gsM^g/N, g'^{N'f^)y and higher, while N^s is the number of effective 
degrees of freedom for the boundary gauge theory 



logro 



where ro is as in section 2.4 i.e. it is the scale where warp factor changes to inverse 
power series from logarithm. Note for duality to hold, must be quite large, making 
iVeff quite large. For Nj^ = M^s = 0, we get back the value of kt as computed in 
[114] . However for a non conformal field theory with fundamental matter - which 
is more relevant for QCD - Mcs ^ 0, N^^ ^ 0, and our analysis thus generates a 
correction to the AdS/CFT result. 



3.1.3 Wake 

In the previous sections we computed the drag force on the quark along with the 
broadening of its momentum. Clearly a moving quark should leave some disturbance 
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in the surrounding media. This disturbance is called the wake of the quark. In order 
to quantify the wake left behind by a fast moving quark in the Quark Gluon Plasma, 
we need to compute the stress tensor T^'^, p,q = 0, 1, 2, 3 of the entire system. We 
expect a cone like disturbance in the medium with the quark located at the tip of the 
cone and this cone becomes apparent when one plots the stress tensor as a function 
of location of the fast quark. We would like to compute the energy momentum tensor 
the medium including the fast quark, then subtract the contribution from the quark 
to obtain the disturbance left behind in the medium. Our goal therefore would be to 
compute: 

rppq _ rppq /o r,A\ 

medium+quark quark yj-'j-i:) 

where the first term is basically the energy-momentum tensor obtained from dual 
geometry considering back reaction of a moving string i.e ^background+string ■ Similarly 
the second term is the energy momentum tensor of the string i.e T^tring restricted 
to four- dimensional space-time. This is similar to the analysis done in [69] for the 
AdS case. For our case the above idea, although very simple to state, will be rather 
technical because of the underlying RG flow in the dual gauge theory side. Our second 
goal would then be to see how much we differ from the AdS results once we go from 
CFT to theories with running coupling constants. 

For a strongly coupled QGP, we will apply the gauge/gravity duality to compute 
T^'^ of QGP using the supergravity action. In the ten dimensional bulk geometry, we 
introduce an additional string moving with some velocity and this string is dual to 
the fast parton creating the wake. The total metric takes the following form: 

kj = lijir,x,y,z,t) (3.35) 

where g^j is the background metric and Uj = 0, .., 9) denote the perturbation from 
the moving string source (with k — )■ 0). In order to compute T^"^, we need to write 
the supergravity action as a functional of the perturbation Ipg. We have O = T^'^ in 
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f l2.5ip and thus 0o = i^^pq is the source in the partition function (12.491) . It follows that 

1 6S 



K Slpq 



(3.36) 



where 5" = S'totai + •S'gh + •S'counterterm as discussed in section 2.3.1 where now 5'totai 
includes the DBI action for the string. 

Minimizing 5'totai , we obtain the equation of motion for Gij : 

R^J {gap + Kl^p) - \ (g.j + Kk,) R {g^fs + K^p) = T^f + + T,^^''^' (3.37) 

where the Tf^'"'^ come from the five-form fluxes F(5) (that give rise to the AdS^^ part) 
and the remnant of the Hnsi Hrr and the axio-dilaton along the radial r direction 
(that give rise to the deformation of the AdS^ part). . The effect of t/J'^''^ will not 
be substantial if we take it as a probe in this background whereas the strings stress 
tensor is given explicitly by 



5S. 



^string(^) — (3.38) 

, , (2X ■ X'X'X'i - X'^X'^X'^ - PxiX''^\ 

'"^"i V^HlT ^""'"^ 

where X is the mapping from string world sheet to space time, and we can consider 
the string profile given by (13. 2p . 

The Einstein equations can be worked out if one considers the effects of all the 
background fluxes in our theory. The result of such an analysis can be presented in 
powers of k. For our case we are only interested in back reactions that are linear in 
K. To this order the equation of motion satisfied by lap is determined by expanding 
(I3.37P in the following way: 

- Btf - A%%, = (3.39) 

where A'^^^,a,P = 0,..,9 is an operator whereas B'^f and A'^f are functions of r, 
the radial coordinatqj. We have been able to determine the form for the operator 



■^There will be another contribution from the (p, q) seven branes in the background, although for 
small gsNf these are sub leading. 
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for any generic perturbation /a/3 in five dimensions that is by setting lij = for 
i,j = 5, .., 9. The resulting equations are rather long and involved; and we give them 
in the Appendix B of [57]. For the functions A'^f and B°f we have worked out a toy 
example in Appendix C of [57] with only diagonal perturbations. For off diagonal 
perturbations we need to take an inverse of a 5 x 5 matrix to determine the functional 
form. We shall provide details of this in the following. The variables defined in fl3.39l) 
are given as: 

/ 6R.ij \ 1 I 5R \ 1 



•^ij - ^ r{b)tibcdar(b)ub'c'd>a'g g g g g 



b,c,d,. 



\ ^ 17 TP „aa'bb'cc'dd'nn'{^ ^nf} 

% - ~g ^{5)nabcdi'i5)n'a'b'c'd'g g g g g {gijg -0^0'^) 

a,b,c,d,.. 

- I E g^,F^'^F|;'^g'^"g'(' + ^ (3.40) 

i=l i=l 

where we have given the most generic form in f l3.40p above for a five dimensional 
perturbation. Furthermore, 1^ and are operators and not functions. 

Once (I3.39P is solved, with the exact solutions lij, we can write write the action as 
a functional of kj. We can then integrate the radial direction along with the internal 
directions and add appropriate Gibbons-Hawkings terms along with counterterms to 
renormalize the action if necessary. The renormalization and calculation of the stress 
tensor was described in some detail in section 2.3.1 and 2.3.2. The final result for 
the stress tensor T^"^ with ulij = (pm is given by (I2.90p . Knowing all the solutions 
lij, one can easily evaluate the stress tensor and finally plot the wake. Here we have 
outlined the procedure and an exact calculation will be done in our future work. 



3.2 Transport Coefficients 

At low energies, an effective theory of fluids is hydrodynamics which describes the 
kinematics of the fluid and uniquely determines its stress energy tensor. As dissipation 
is allowed, the theory is not described in terms of action but in terms of equation 
of motion, namely the equation which guarantees conservation of energy-momentum 
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tensor 



V^Tf"" = 



(3.41) 



where iJ^iy = 0, 1, 2, 3. There are four equations above and they can be solved by four 
independent variables |116j |117]|118] . If we assume local thermal equilibrium then the 
state of the system can be uniquely determined by the local fluid velocity u^^{x^) and 
local temperature T{x^). There are actually three independent components of the four 
velocity as u^u^ = — 1 fixes one of the components and the fourth independent 
variable is the temperature. Using these four independent variables, we can write 
the expression for stress tensor that satisfies equation 03.411) and at zeroth order in 
derivatives of we obtain the stress tensor of ideal fluid: 



where e is the energy, P is the pressure and g^y is the metric of four dimensional 
space time where the fluid lies. 

Now at linear order in derivatives, the allowed terms are restricted by rotational 
symmetry and only nonzero components are t|^', i, j = 1, 2, 3 



where Pia = Qia — UiUa and the coefficients of the two terms determine the transport 
coefficients shear viscosity t] and bulk viscosity ( tllZj- Thus f l3.43p is the defining 
equation for viscosity. 

Physically shear viscosity measures the mixing between two layers of a fluid [57]. 
More viscous the fluid, the faster momentum can be transferred from a layer to the 
next. Somewhat counter intuitive is the fact that the stronger the coupling in the 
fluid, the less the shear viscosity. This is because the rate of mixing is controlled 
by the mean free path. When the mean free path is small compared to the flow 
velocity variation of the two laminas, the layers cannot easily mix since the exchange 
of particles is limited to the small volume near the interface of the two laminas: Most 
particles in the fluid just flows along as if there is no other layers nearby. On the 



rPj = (e + P)u,uy + Pg, 



(3.42) 




(3.43) 
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other hand, if the mean free path is comparable to the typical size of the flow velocity 
variation, then mixing between different layers can proceed relatively quickly. For an 
ideal fluid, laminas of fluid do not interact at all and we have zero viscosity. 

Transport coefficient such as shear viscosity becomes crucial in understanding 
physics of quark gluon plasma. In particular, as already mentioned in the introduc- 
tion, in a heavy ion collision the overlapping region of the two nuclei is of elliptical 
shape with different short and long axes. Thus the plasma formed undergoes elliptic 
flow due to difference in pressure along the long and short axes. On the other hand, 
a large viscosity would mean greater interaction between the layers of fluid which 
would quickly equilibrate the system and one would observe very little elliptic flow. 
But the data from RHIC is well described by ideal hydrodynamics with zero shear 
viscosity and shows strong elliptic flow P!n]-|22]. This means the fluid created in 
heavy ion collisions has small viscosity and thus is strongly coupled. Thus one can- 
not apply conventional perturbative field theory techniques to compute the transport 
coefficients of the plasma. 

However in the regime of strong coupling, certain gauge theories can be described 
by dual supergravity. In particular for the gauge theories described in section 2.2- 
2.4, which have several features common to QCD, one can compute the transport 
coefficients at strong coupling using dual geometry. The result we obtain can then 
be suggestive of the viscosities of QGP. In the following sections we compute shear 
viscosity r] for gauge theories using dual supergravity and the ratio rj/s which appears 
in various experimental observations. 

3.2.1 Shear Viscosity 

In this section we will present our calculation [57] of shear viscosity of the four di- 
mensional theory following some of the recent works |146]|12U] . The shear viscosity 
described earlier can be obtained from correlation functions using the Kubo formula 



V = lim — / dtd^x e'^'{[T23{x),T23m) = - lim 



Im G 



(3.44) 
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where G^{uj, if) is the momentum space retarded propagator for the operator T23 at 
finite temperature, defined by 

G\u,t) = -^J dtd'x e^(-*-^-^)e(t)([T23(a;),T23(0)]) (3.45) 

In the following, we will compute the Minkowski propagator following the conjecture 
made in |121j [122] using dual gravity. Note that our prescription fl2.49p computes a 
path integral with a classical action Ssugra, unaware of the ordering of the operators 
whose expectation value is being computed. Therefore computing any commutator is 
subtle here. Hence, we compute only the correlator (T23(r, af )T23(0, af )) and using 
the conjecture in |121] . relate this to the retarded Greens function. However before 
we compute this explicitly, let us evaluate the higher order corrections to the effective 
action from the wrapped D7 brane in our theory. 

In the case of a single D7 brane, the disc level action contains the term [123] [124] : 

= — ^ . —^—^ I \Ci At! (RAR)- e-hi (R A *R)] (3.46) 
°^ 1927r^, (47ra')2 Jm^^ ^ ' ^ ^ ' 

where C4 is the four-form, R is the curvature two-form, (f) is the dilaton and is 
a non-trivial eight manifold which is the world-volume of the D7 brane. The action 
is SL{2, Z) invariant which was shown by doing an explicit analysis [123] |124] . Since 
for our case the D7 wrap a non-trivial four-cycle, we can dimensionally reduce it over 
the four-cycle and obtain the following action: 

= I Re [log r/(r) tr {RA*R-iRA R)] (3.47) 

where 77 (r) is the Dedekind function, and r is the modular parameter defined as 
follows: 

with V4 being the volume of the four-cycle on which we have the wrapped D7 brane. 
In fact the above action can be derived from the following action that has two parts, 
CP-even and CP-odd [125] : 

/ log|r^(r)|2tr {RA*R)--^ [ log ^ tr (i? A R) (3.49) 
327r^ JAf4 327r^ Jm* //(r) 
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where the first part is CP-even and the second part is CP-odd. To compare f l3.49p with 
fl3.47p note that the Dedekind 77 function has the following expansion with q = e^'^*'^: 

.2 A ^3 



log|r?(r)p 



6 



r2 



3q' 4g' 

— + — + ... + C.C 



log 



r]{T) 



3g2 4g3 



c.c 



(3.50) 



Combining everything we see that, up to powers of q f l3.49p and f l3.47p are equivalent. 
However writing the action in terms of f l3.49p instead of fl3.47p has the following 
advantage: from D7 point of view fl3.49p captures the D3 instanton corrections in the 
system |125j - [T29] . But a deeper reason for writing the action as f l3.49p is that the 
CP-even and CP-odd terms can be expanded further to account Gauss-Bonet type 
interactions |125] : 



^CP-even = -«1 / 6 "^^GB - ^7 / e" 



(3.51) 



where ai is a constant, and £gb and Cr are respectively the Gauss-Bonnet and the 
curvature terms defined in the following way: 



GB 



R 



327r2 



p Tya/S-fS r) p pa/3 p jyairfS i o p p^ 

^aS^S-ti — ZKaprC — Kab-yS-ti + ^riabrC 



ab 



(3.52) 



In the above note that the three curvature terms RapR"^ , RatR"^^ and R^ are not the 
pull-backs of the bulk Ricci tensor. We have also used the notations {a,/3) to denote 
the world- volume coordinates, and (a, b) to denote the normal bundle. 
From the CP-even terms, the coefficient of Ra/s-ysR"'^'^^ is given by [57]: 



C3 



327r2 



(3.53) 



which has an overall plus sign because ai in many cases is zero (see |125j for a 
discussion on this). However in general for certain exotic compactifications we can 



have ai « 



■. If we now compare this to |146j we see that C3, which is the 



coefficient of Raj3-ysR°'^'^^ in [146] . is indeed positive. This would clearly mean that 
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adding fundamental flavors lowers the viscosity to entropy boundl^ 

The CP-odd term on the other hand has a standard expansion [131j - |134] . |123] - 

(2 '2 \ 
Cs + ^^C,AtT ear] (3.54) 

where the first term gives the dual axionic charge of the D7 brane. Combining (13.511) 
and (I3.54P we get the full back reactions of the D7 brane up to 0{a"^). 

Having computed the back reactions of the embedded D7 brane, we can use this 
result to compute the shear viscosity. To start the analysis, we need the correlation 
function of T23(x) and T23(0) to use it in the Kubo formula (13.441) . From gauge/gravity 
duality we know that switching on T23 in the gauge theory is equivalent to considering 
graviton modes along = x and x^ = y directions. In the ten dimensional dual 
geometry with the metric of the form (12.911) . we introduce graviton perturbations 
in the xy direction. We start with the ten dimensional SUGRA action with the 
perturbed metric and then integrate over the five internal compact directions to obtain 
an effective five dimensional action . The resulting five dimensional metric dsf = 
g^udx^^dx" which minimizes the five dimensional action takes the following form: 
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where we have ignored 0{ggNj), k > 2 terms (which is equivalent to setting g^r = 1 
in (12.911) and thus in five dimensions the effective grr is 1) and h{r) is now only 
a function the the radial coordinate. Note that h is obtained from h appearing in 
(I2.9ip by integrating over the internal coordinates on which h depends. When h is 

^The analysis here was motivated by discussions with Aninda Sinha. His paper [130] dealing 
with the violation of viscosity to entropy bound appeared recently and has some overlap with our 
analysis. 
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independent of any of the angles of internal space, like the case in AdS space, then 
h = h. We have also taken the approximation gi = g2 = g{r) for simplicity and this 
approximation becomes exact if we ignore 0{gsNf,gsM'^/N) corrections to the black 
hole factors gi. 

Since our goal is to compute the Fourier transform of (T23(t, a^)T23(0, a^)), we can 
do this by first writing the supergravity action in momentum space, treating it as a 
functional of Fourier modes for 0(r, t) where: 



(j){r,t) = (p{r,t)(t){t) = J du e-^'^"0(r,w) = J du e-''^^'(j){r,u) 

(f){r,u) = ^{r,\uMu) (3.55) 

where as before, we defined the Fourier transform using the curved space time r = 



g{rc) t and not simply t. Although this definition is precise for the theory at the 
cut-off r = Tc only, we will use it also for any r because in the end we will only provide 
description at the boundary (i.e r — t- oo) where the results would be independent of 
the choice of the cut-off. 

The way we proceed now is the following^. We consider the metric fluctuation 
as in fl3.2.ip and plug this in the five dimensional effective action. Finally, we will 

(2) 

call this resulting action as S^q where the subscript (2) involves writing the action 
in terms of quadratic (j){r,u). We do this as there exists a very useful relation for 
computing the shear viscosity (see for example |121j|122] ): 

lim Im Gf^{u, t) = lim — Im G^{u, t) (3.56) 

where Gfj^ is the Schwinger-Keldysh propagator [121] - [T35] . Comparing this with our 
earlier Kubo formula (I3.44p . we get the following expression for shear viscosity: 

V = -^lim lrciGff{u,t) (3.57) 
Thus if we can write our effective supergravity action in the following way: 

si'Mn,uj)] = lJ^Mn,u;)Gf{\u\,t)<PAn,-uj) (3.58) 



^This is similar to the procedure of |146| . Notice however that the theory considered by [146) has 
no running but contains higher curvature-squared corrections. 
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where ri, is a specified point on the boundary. Then taking the Gff{u, if) part and 
using ( I3.57P we can easily obtain the shear viscosit}o. In other words, we will be 
taking two functional derivatives of S'gQ[0(rb, w)] with respect to (j){rb,ui) and thus 
are interested about terms quadratic in (f){rb, oj) in the action. Of course, in real 
time formalism, we are concerned with the Schwinger-Keldysh propagator Gfj^ of the 
doublet fields (j)i{r,t), (f)j{r,t). In the context of gauge/gravity duality, we follow the 
procedure outlined by [135j for AdS/CFT correspondence and treat 4>i{r,t), (j)2{r,t) 
as the perturbation (f){r,t) and its doublet in the four dimensional Minkowski space. 
In ten dimensional gravity theory, 0i(r) = (j){r) is the field in the R quadrant of 
the Penrose diagram and 02('") is the field in the L quadrant. For more details see 
[T35]-[TIOj. 

To be more precise, our aim is to get the effective action in the form (13.581) . To 
this effect we take our metric fl3.2.ip and plug it in the five dimensional effective 
action with net result: 



4?1 = /^/ dr 



dud^q 



+ C{r)(f){r, —u)(f)'{r, u) + D{r)(f){r, —u)(j){r, u 



A{r)(j){r, -uj)(f)"{r, oj) + B{r)(j)'{r, -uj)(p'{r, u) 



(3.59) 



where prime denotes derivative with respect to r and the explicit expressions for 
A, B, C, D are given in Appendix E of [57j. The five dimensional Newton's constant 
is given by: 

Gn - (3.60) 

where volume of T^'^ i.e Vt^,! is dimensionful and is proportional to ten dimen- 
sional Newton's constant. 



^Notice that there would be an overall volume factor of T^'^ that would appear with the effective 
action. This factor just modifies the Newton's constant in five dimensions and does will not effect 
dimensionless ratios like 77/s, as we shall find out. 

^ Although our background is a deformation of the AdS space, the arguments of [135j also apply 
here. We can still consider the 0i perturbations to compute the Schwinger-Keldysh propagator 
because by definition a propagator is what appears sandwiched between the fields. See also [140] for 
a generic approach. 
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The fluctuation (/)(r, u) satisfies the following Euler-Lagrange equation of motion: 

,/// N ^'W-^'W,,/ X 2D(r)-C'(r) + A"(r),, , 

which we can derive from 13.59] by minimizing it. Once we plug in the values of A, B 
etc., the above Euler-Lagrange equation takes the following form: 



^ ^ + -+M{r 



9{r) 



r4 f 

W = -T 1 



2-kN 



^gsNf 



u'^g{rc)h{r) 



9[r) 



2tx 



logr + - 



Ait 



logr 







(3.62) 



where J7'(r) and Ai{r) appear due to seven branes and fluxes in the geometry 
As before, primes in fl3.6ip and f l3.62p denote derivatives with respect to the five 
dimensional radial coordinate r. 

Now as we mentioned above in f l3.55p . (j){r,u) can be decomposed in terms of 
0(r, \u\) and Then as a trial solution, just like in [146J . we first try 0(r, cj) = 

g{ry and look at f l3.62p for r near the horizon where g{r) — )■ 0. Plugging this in 
( I3.62P with g{r) =0 we obtain : 



7 



±?|ci;| 



'h{rh)g{rc) 



16 



±i- 



\UJ\ 



(3.63) 



where in the last step we have used the definition of temperature Tc as in fl2.45p . 

To get the solution with for general r, where g{r) ^ 0, we propose the following 
ansatz for the solution to (I3.62p : 

(t){r,uj) = g{rf'^cF{r,\u\)(f){uj) (3.64) 

Plugging this in fl3.62p we see that the equation satisfied by F(r, |c<;|) can be expressed 
in terms of 7 and 7^ in the following way: 

'\^?9{rc)h 



F"{r, \co\) + ( ^ + - + Mir)] F'{r, \co\) 
\g{r) r J 



9'^{r) 



+ J{r)\F{r,\u\) (3.65) 



-"^^In special cases we expect J'(r) and A4 to vanish (see for example |141| V However when this 
is not the case ,as possible for non-trivial UV completions of our model, we could expect a non- 
minimally coupled scalar field. 
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+ o4^nr.i-i) 



9"{r) 
9{r) 



5 , A g'(r) 
r / g[r) 



F(r,|u;|)|+7'^F(r,|u;|) = 



where the 7^ terms come from both the last term in the above equation as well 
as the |a;p term above. Furthermore, note that the source J{r) ~ 0{gs) + 0{g1), 
so in the limit — we find that (13.651) has a solution of the form F(r, |a;|) = 
Ci + C2g{r)~'^'^ with ci, C2 constants. Then we expect the complete solution for gs 7^ 
to be F(r, = Ci + C2g{r)~^'^ + f{r, \uj\). Demanding that F(r, be regular at the 
horizon r = r^ forces C2 = as g{rh) = 0. We choose Ci = 1 and / = ^+7^+7^/0 + ... 
as a series solution in 7. Then our ansatz for the solution to (13.651) becomes 



F(r, = 1 + Qir) + -^n{r) + 72/C(r) + ... 



(3.66) 



Once we plug in the ansatz (I3.66P in (13.651) we see that the resulting equation can be 
expressed as a series in 7: 



g r 



9 



9" /5 , A g 

r 



9 



y \9 r ) g 
+ (-o + ^)(l + ^)} 



31 



+ 7' ^ic 

I 9 



12' 



9 



+ 0(7^ 



(3.67) 



where we have avoided showing the explicit r dependences of the various parameters 
to avoid clutter. We have also defined kq in terms of the variables of (13.631) in the 
following way: 

16 



T'g' 



(3.68) 



Although the above equation (I3.67P may look formidable there is one immediate 
simplification that could be imposed, namely, putting the coefficients of 7^,7,7^,... 
individually to zero. This is possible because one can view 7 to be an arbitrary 



100 Chapter 3 Application to Thermal QCD 



parameter that can be tuned by choosing the graviton energy cu or the temperature 
Tc. This means that the zeroth order in 7 we will have the following equation: 

g'{r) 5 



+ M(r) 



g'{r) + J{r)[l + g{r)] = 



(3.69) 



,gir) r 

In the above equation observe that the source J{r), for cases where it is non-zero, has 
a complicated structure with logarithms and powers of r. To simplify the subsequent 
expressions, let us choose to work near the cut-off r = Tc. This is similar to the spirit 
of the previous section where we eventually analyzed the system from the boundary 
point of view. Then to solve (13.69^ near r ~ Tc we can switch to following coordinate 
system 



r = rc(l - C) 



(3.70) 



Taylor expanding all the terms J{r),g{r), ^n(^i_Qn in f l3.69p about C = 0, we obtain a 
power series solution for Q as: 



g{r) 



EE 



a i=0 "^cla) ( ^ 



E «^c^ 



(3.71) 



Since f l3.69p is a second order differential equation, we can fix two coefficients and we 
choose = 0. Then the rest of a'^s are determined by equating coefficients of 

on both sides of equation f l3.69p . The exact solutions are listed in Appendix E of 
\57\ . Note that all are proportional to Qs and in the limit Qs 0, Q ^ 0. 

To next order in 7 we have an equation for "H that also depends on the solution 
that we got for Q. The equation for 'H(r) can be taken from fl3.67p as: 

2g'{r) 



n"(r) + 



,g[r) r 



H'ir) + J{r)H{r) 



g{r) 



g"{r) 
g{r) 



+ 



+ M{r) 



gin 
g{r) 



[i + g{r) 



(3.72) 



To solve this we make the coordinate transformation f l3.70p and plug in the series 
solution for g{r) given above. The final result for "H can again be expressed as a 
series solution in ( in the following way: 



nr) = EE 



4i 



a i=0 c{a) 



:i-c) 



4i 



E^^c 



(3.73) 



i=0 
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We again set 60 = = and following similar ideas used to solve for Q, we determine 
all b[s by equating coefficients inf l3.72p . The exact solution is given in Appendix E of 
[57] . Again note that all bi are of at least 0{gs) and thus with — )■ 0, "H — )■ 0. 

Finally the second order in 7 is a much more involved equation that uses results 
of the previous two equations to determine /C. This is given by: 

2g'{r) 



IC"(r) 



+ - +Mir)] /C'(r) + :7(r)/C(r) 



9[r) 



gij) 



-H'{r) (3.74) 



9"ir) 
.9{r) 



+ 



9'{r) 
9{r) 



9'\ry 
9^ir) _ 



[1 + g{r)] 



which could also be solved using another series expansion in (we haven't attempted 
it here). Therefore combining f l3.7ip and f l3.73p we finally have the solution for the 
metric perturbation: 



(/>(r, \u\)± = ^(r)^*4.T, 



l + g{r)±i- 



\U!\ 



■Uir) - 



-Mr) + 



(3.75) 



We can analyze this in the regime where the gravitons have very small energy, i.e 
— 7- or equivalently 7 — > 0. In this limit we can Taylor expand (f){r, \u\) about 
7 = to give us the two possible solutions: 

Iwl 



1 + g{r) ± ^^{nr) + [1 + g{r)]\og g{r)} 



(3.76) 



\UJ\ 



;{lC{r) + H{r) log g{r) + [1 + g{r)]\og^g{r)] + 0{\uj\ 



167r2T2 

which consequently means that to the first order in u the off diagonal gravitational 
perturbation at low energy is given by two possible solutions corresponding to positive 
and negative frequencies as: 



[1 + g{r)] 0(u;) ± ^^{H{r) + [1 + 6;(r)]log g{r] 



(3.77) 



As is well known, following |142j |135j . we can define field on the right R and left L 
quadrant of the Kruskal plane in terms of 0+(r, uj) and 0-(r, u) in the following way: 



0R,±(c^,r) 



(j)±{uj,r) in R 
in L 
(t>±{uj,r) in L 
in R 



(3.78) 
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Now </)R,±, (/)L,± contain positive and negative frequency modes but a certain linear 
combination of 0r,±,0l,± gives purely positive or purely negative frequency modes 
in the entire Kruskal plane |142] |135] . Furthermore imposing that positive frequency 
modes are in falling at the horizon in R quadrant and negative frequency modes are 
outgoing at the horizon in R fixes two combinations : 



•'pos 



■'neg 



<j)R,4u,r) + e^/^''^<f)i.,4u,r) (3.79) 



With fl3.79p we see that we can define fields in R(L) quadrant as linear combination 
of positive and negative frequency modes 

0R(w,r) = ao[0R,+ (w,r) - e'^^^'^^R r)] = 0i 

M^,r) = aoe'^/2^^[0L,+ (c^,r)-0L,-(c^,r)] =02 (3.80) 

where we have identified 0r(0l) with the thermal field 0i(02) defined on the complex 
time contour which familiarly appears in the Schwinger-Keldysh propagators of real 
time thermal field theory. Here Sq is a constant. The final physical quantity that we 
will extract from here will only depend on T, as we will show soon. 

Having got the graviton fluctuations (f){r,u) = </)R(a;,r), we are almost there 
to compute the viscosity rj using fl3.57p . Our next step would be to compute the 
Schwinger-Keldysh propagator Gff{0,~(f'). All we now need is to write the action 
fl3.59p as fl3.58p and from there extract the Schwinger-Keldysh propagator. This 
analysis is similar to the one that we did in the previous section, so we could be brief 
(see also |120j ). The action (13.590 can be used to get the boundary action once we 
shift (f){r, u) to (j){r, u) + S(f){r, oj) in the following way: 



4S( 



+C(r)0(r, -a;)0'(r, u) + D{r)<p{r, -a;)0(r, uj) + [2A(r)0"(r, uj) (3.81) 
-25(r)0"(r, uj) - 2B\r)4>'{r, u) - C"(r)0(r, u) + 2D(r)0(r, co) 
+A"(r)0(r, to) + 2A'(r)0'(r, u)] S^{r, -co) + dr [25(r)0'(r, w)50(r, -tu) 
+C(r)0(r, uj)6(f){r, —oo) + A{r)(j){r, uj)6(j)'{r, —oo) — dr {A{r)(j){r, oo)) 6(j){r, -co) | 
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Plugging in the background value of 0(r, u) will tell us that only the boundary term 
survives. And as before, to cancel the A{r)(f){r, uj)6(f)' {r, —u) we will have to add the 
Gibbons-Hawking term to the action ^61j. The net result is the following boundary 
action: 



(2) 
SG 



1 



(2vr)4 
dujd?q 



r, -u) 



^2 

J'(a;,r) 



C{r)-A'{r) + B{r)-A{r) 



(r, -u) 



(j){r, u 
(3.82) 



rh 



rh 



Now comparing f l3.58p with 13.821 we see that the terms between the braces combine 
to give us the required Schwinger-Keldysh propagator: 

1 J^{cj,r) 



lim 



lim - 



(3.83) 



rh 



N\ 9[rc) 



G{r) - A'{r)] + \B{r) - A{r) 



01 (r, -w) 



where we assum^iU that 0i(r/j, w) = ao[0R,+('^/i, i^) — e'^^"^''0R,-('^fe, i^)]- Now to eval- 
uate the shear viscosity from the above result we need to perform two more steps: 

• Evaluate the contributions from the UV cap that we attach from r = to r = oo. 

• Take the imaginary part of the resulting total Schwinger-Keldysh propagator. This 
should give us result independent of the cut-off. 

To evaluate the first step i.e contributions from the UV cap, we need to see precisely 

(2) 

the singularity structure of S^q. The second step would then be to extract the 
imaginary part of SK propagator from there. Since the imaginary part can only come 
from the second term of (I3.82p . we only need to evaluate: 

1 r„ i0'i(r,-u;) 



lim - 



N\ 9[rc) 



B{r) - A{r) 



(r, -uj) 



(3.84) 



^^At this point one might worry that the solution for 0i is only known around Tc- That this is 
not the case can be seen in the following way: Integration by parts gives p.83p which says one only 



needs to know the value of the field 0i at Vc and rh. The solution for 0i = ^ is given in (I3.75P 
from which it is clear that 4>i{rh) = as g{rh) — 0. Furthermore to know rj we only need to know 
the imaginary part of p.83p . which is evaluated using (|3.85l) in (|3.83p and using boundary values of 

4>i[rc) and (t)i{rh). 
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with 0(r, —u) being the graviton fluctuation in the regime r > Tc. To analyze this let 
us flrst consider a case where gs ^ and (^(r), 'H(r), /C(r), ...) — > 0. In this limit we 
expect for < r < Tc- 

1 



B{r) - A{r) 



'^91 



g{r)r'' + 0{gsN^ 



f) 



(3.85) 



ao 
g'{r) 



+ + i^iog gii + G) 



4nTr 



2'KTr 



2nT, 



27r 



0i(r, -w) g{r) \A7r'^ + \og'^g{r) J 

The above considerations would mean that the contribution to the viscosity, t]i, for 
this simple case without incorporating the UV cap will be: 



2nT,g^,GNJg(^ V47r + ^ log^g{r,)J 2glGM \Ati + i \og^g{ 



(3.86) 



rhJh{rh) 



and h{rh) ~ ^ in this 



where we have used the relations t^Tc\j gir^ 
limit. This helps us to write everything in terms of T and not the scale dependent 
temperature T^.. In fact as we show below, once we incorporate the contributions from 
the UV cap, the Tc dependence of the above formula will also go away and the final 
result will be completely independent of the cut-off. Note that in the limit Tc — )■ oo 
we recover the result for the cascading theory. 

Combining all the ingredients together, the contribution to the viscosity in the 
limit where (^(r), 'H(r), /C(r), ...) etc are non-zero can now be presented succinctly as 
(although rji below doesn't have any real meaning on the gauge theory side as this is 
an intermediate quantity): 

" m' 

nx/nyrh) I -^^^ T+g~'{ 



(3.87) 



2g^,GN |^47r + i[log(?(r,) + ^ 

Note that the above expression is exact for our background at least in the limit where 
we take the leading order singularity of the background. This is motivated from 
our detailed discussion that we gave in the previous section. Note that the second 
term in the action (13.821) is exactly the second equation of the set (12.691) whose 
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singularity structure has been shown to be renormahzable. Thus taking the leading 

(a) 



order singularity instead of the actual rf^ will not change anything if we carefully 



compensate the coefficients with appropriate gsNf,gsM'^/N factors! 

But this is still not the complete expression as we haven't added the contributions 
from the UV cap. Before we do that, we want to re-address the singularity structure 
of the above expression. The worrisome aspect is the existence of factor in fl3.87p . 
Does that create a problem for our case? 

The answer turns out to be miraculously no, because of the form of H and Q given 
in fl3.73p and (13.71 p . This, taking only the leading powers of rc, yields: 

n' = -^-^ + ...., g' = -^-^ + .... (3.88) 

c c , , c c 



killing the dependence in This would make r/i completely finite and all 

the Tc dependences would go as 0{l/rc). Therefore we expect the contribution to the 
viscosity from the UV cap to go like: 

m^v\-=T.^. (3-89) 

where the total viscosity will be defined as r/ = r/i + 772- As this is a physical quantity 
we expect it to be independent of the scale. Therefore 

which will give us similar Callan-Symanzik type equations, as discussed in the previous 
section, from where we could derive the precise forms for Gi in 13.891 Finally when 
the dust settles, the result for shear viscosity can be expressed as: 

1 + Er=i«fce-^'-^- 



?7 = 



N 



(3.91) 



_47r + i log^(l -r^e-^"-), 
where Uk are functions of T that can be easily determined from the coefficients 
{ai,bi) in fl3.7ip and f l3.73p or {ai,bi) worked out in Appendix E of [57]; and hiT) = 



-"^^It is now easy to see why 61 ai = is consistent. For non-zero 61, ai there would have been 
additional log r terms from r^^j . These would have made the theory non-renormalizable. Thus 
holographic renormalizability would demand 61 = oi = from the very beginning — consistent with 
what we choose earlier. 
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Yi + 0{gs, Nf, M). Observe that the final result for shear viscosity is completely 
independent of Vc and Tc, and only depend on T and the degrees of freedom at the 
UV i.e through e~-^"". Needless to say, for large enough A/'^^ (which is always the case 
for our case because Afuv — > e~", ?2 > 1), the shear viscosity is only sensitive to the 
characteristic temperature T of the cascading theory. The interesting thing however 
is that the shear viscosity with finite but large enough Afuv can be smaller than or 
equal to the shear viscosity with Afuv e^",n >> 1 i.e for the parent cascading 
theory provided: 

in the limit of small characteristic temperature T. This will have effect on the viscosity 
to entropy ratio, to which we turn next. 

3.2.2 r]/s 

Going back to the stress tensor of fiuid described by hydrodynamics, one can obtain 

m 

m,) = -j^ (V.(T,°) + V,{T^) - ^5,,VKT'°)) - ^%V,(T'°) (3.93) 

where 6T^^ is the deviation from the ideal fiuid stress tensor T*-' in the fiuid rest frame 
and e and P are the local energy density and the pressure, respectively. Using the 
thermodynamic identity Ts = e+P where s is the entropy density, the two coefficients 
can be also written as rj/Ts and C/Ts. Since the temperature is the only relevant 
energy scale in the highly relativistic fiuid, one can easily see that the importance of 
the viscous terms depends on the size of the dimensionless ratios r]/s and (/s. Thus 
not the absolute value of shear viscosity rather the ratio rj/s is what is relevant in 
viscous hydrodynamics. 

Furthermore, the dimensionless number rj/s has been shown to be universal for a 
large class of gauge theories where the dual graviton perturbation is a scalar field |143] 
and even more strikingly r]/s = I/Att has been conjectured to be its lowest possible 
value |144] for any fiuid using arguments of ^145j . Although it has been shown by 
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|57j|146j - |154j that the lower bound can be violated for gauge theories which have 
higher derivative terms (i.e. curvature square terms and beyond) in their dual gravity, 
these theories may violate other constraints and may not be most relevant for physics 
of fluids. In any case, the ratio t]/s becomes crucial in describing fluid dynamics and 
in this section we will compute this ratio using dual gravity. 

We will calculate the ratio for two cases i.e one with only RG flow, and the other 
with both RG flow and curvature squared corrections. As usual the former is easier 
to handle so we discuss this first. 

Starting with the type JIB supergravity action (I2.19P in ten dimension, the entropy 
is given by Wald's formula [155] . [156] . [T57] ■ [158] 

dxdydzd^MVV—-^e''''e'"^ (3.94) 

(J^abcd 

where the integral is over the eight dimensional surface of the horizon at r = rh,Cio 
is the Lagrangian density of the action in fl2.19p . Vat, a,b = 1..8 is the induced 8x8 
metric at horizon, eab is the bi normal normalized to eab^"''' = —2. Finally using explicit 
expression for the metric f l2.9ip with warp factor given by fl2.3ip near horizon, we have 



TF = 2 Uxdydz d'M sin 9, sin 9^ JhiJ^Zo^) ^^e^'e^'^ 

^3 108V3fi;io oRabcd 



2g^GN 1 2nN 



log rh + - 



1/2 



log (3.95) 



27r V 2/ 47r 

where V3 is the infinite three dimensional volume and we have used the definition of 
five dimensional Newton's constant Gn introduced in f l3.60p . 

Once we replace rh by the characteristic temperature T, we see that the entropy 
is only sensitive to the temperature and is independent of any other scale of the 
theory. Since the above result is also independent of Nuv it would seem that the 
Wald formula only gives the entropy for the theory with Muv = 00 i.e for the parent 
cascading theorjo. The interesting question now would be to ask what is the entropy 
for the theory whose UV description is different from the parent cascading theory? 
In other words, what is the effect of the UV cap attached ai r = on the entropy? 



""^^This can be argued by observing that fact that in a renormahzable theory, hke ours, the depen- 
dences on degrees of freedom go hke O (e~^"") corrections as we saw in the previous sections. 
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To evaluate this, observe first that in finite temperature gauge theory, entropy 
density of a thermalized medium having stress tensor (T^^) = diagonal(e, P, P, P) is 
given by 

^"^^ (3.96) 



T 

where e is the energy density, P = Px = Py = Pz is the pressure of the medium and 
T being the temperature. With our gravity dual we can compute the stress tensor 
earned) (^^"^ ^h^s the energy e = (Tj^ed) the pressure P = {T^^^)) of the medium 
through equation of the form (13.361) . i.e 

(TZd) = ^ (3.97) 

Obgpq 

where again p,q = 0, 1, 2, 3 and gpg is the four dimensional metric obtained from the 
ten dimensional OKS-BH metric gij,i, j = 0,1. .,9; and 6f, operation has been defined 
earlier. There are two ways by which we could get a four- dimensional metric from 
the corresponding ten-dimensional one. The first way is to integrate out the 6i,(j)i 
directions to get the four-dimensional effective theory. This is because the warp factor 
for our case is dependent on the 6i directions. The second way is to work on a slice 
in the internal space. The slice is coordinated by choosing some specific values for 
the internal angular coordinates. Such a choice is of course ambiguous, and we can 
only rely on it if the physical quantities that we want to extract from our theory is 
not very sensitive to the choice of the slice. Clearly the first way is much more robust 
but unfortunately not very easy to implement. We will therefore follow the second 
way by choosing the the five dimensional slice a.s 6i = 02 = 7^,1^ = 4'i = 4'2 = and 
thus obtaining 



= g^.u{9^ = ^l,i^ = (pi = 0) (3.98) 



with /i, z/ = 0,1,2,3,4. The next step would be to evaluate all the fiuxes and the 
axio-dilaton on the slice. To do this we define: 

IH3I' = \H,\\e, = -K,^ = (f), = oy, 

IF5I' = |F5|2(e, = 7r,V^ = 0, = 0); 
$ = $(^i = 7r,7/; = 0, = 0) (3.99) 



= 




9i = TT,tp = 


(P^ = 0) 




\m 


9i = 7!-,^p = > 


^i = 0) 
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Once the fluxes have been defined, we need the description for Stotai which gives rise 
to fl3.98p . This is easily obtained from f l2.19p as: 



^ J (fx e-2*v^(R - 49,$9^$ - ^|H3 



'total n 9 



1 f / / 1 1 O 1 1 o 1 



d'x^l \F,\' + \F,\' + -iFsl^ (3.100) 



with R being the Ricci-scalar for g^^, and g = det g^j^u- Note that in the definition 
for the slice sources H3, Fi, F3, F5 and R, we still have gij, i,j>5 which we evaluate 
at 6i = 7t,iIj = (pi = 0, treating them simply as functions and not metric degrees of 
freedom. 

To complete the background we need the line element. Here we will encounter 
some subtleties regarding the choice of the black-hole factors and the corresponding 
QsNf type corrections to them. With the definition of g^u the line element is: 
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^4 00 
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^ i,j=0 

where aij,/3ij are all of 0{gsNf, ggM) and only involve the parameters of the theory 
namely, r^, L and /i from the embedding equation f l2.16p and guarantees that ^ 
are dimensionless. On the other hand crj"-*, k^"^ can incorporate zeroth orders in gsNf. 
However note that so far we have been assuming gi{r) ^ g2{r) = gi'f'), ignoring their 
inherent 6i dependences, and also the inequality stemming from the choices of aij and 
Pij. This will be crucial in what follows, so we will try to keep the black hole factors 
unequal. These considerations do not change any of our previous results of course. 

Now looking at the form of the metric, knowing the warp factor h{r, vr, vr) and 
gi{r), just like before we can expand the line element as AdS^ line element plus 
0{gsNf, g^M) corrections. We can then rewrite the line element fl3.10ip as: 

^2 \ hiv Ti vr) 

ds^ = -L^[g[r) + h]dt'+ ^ ^' ' U r^ + — {l + l2)d-^^ 

g2{r) 
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Mr) ^ t 0.'-'<^' 



i,j=0 



oo 



Hr) = E (3-102) 

i,j=0 



where again dj are of 0{gsNf, QgM) and we are taking h{r, vr, vr) = ^+0{gsNf, QsM). 
Such a way of writing the local line element tells us that there are two induced four- 
dimensional metrics at any point r along the radial direction: 

g(;) = diagonal (-(7(r), 1,1,1), g« = di^gonsX {-hMMM) (3.103) 

where we haven't shown the r"^ / L'^ dependences. The reason for specifically isolating 
the four- dimensional part is to show that we can study the system from boundary 
point of view where the dynamics will be governed by our choice of the boundary 
degrees of freedom. It should also be clear, from four-dimensional point of view, the 
metric choice g^^-* is directly related to the AdS geometry whereas the other choice 
g^q^ is the deformation due to extra fluxes and seven branes. 

The above decomposition also has the effect of simplifying our calculations of the 
energy momentum tensor (Tj^ed)- "^^^ rewrite the total energy momentum tensor 
as the sum of two parts, one coming from the AdS space and the other coming from 
the deformations, in the following way: 



jrppq \ ^'O" total 

\ mcd/ r _.n 



'total 



.pq \ _ ^fcS[°l total _^ '^bS[Jtai 
^b^pq ^bg.lq 
(^mcd)AdS + (rmed)def 

sill + Sll, (3.104) 



where is zeroth order in gsNf,gsM and S[li is higher order in gsNf,gsM. 

Note that (Tj^cd)Ads = r is the well known AdS/CFT result obtained from the 
analysis of [B3]-|nZ] in the limit — t- oo. With the 0{l/r) series expansion of our 
metric ggg = 1 — r'^/r^, = = §33 = 1; ^^e result at the boundary is 

/rjlOO \ _ ' h _ I 

\ med/AdS 2^2^^^ 2g^^GN 

(^rncd)Ads = m,n = 1,2,3 (3.105) 
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This only gives the CFT stress tensor as we evaluate the tensor on the AdS boundary 
at infinity, reproducing the expected first term of f l3.95p . How do we then evaluate the 
0{gsNf, QsM) contributions from the deformed AdS part i.e the energy momentum 
tensor (r^ed)def at any r = Tc cut-off in the geometry? 

In fact the procedure to evaluate exactly such a result has already been discussed 
in section 2.3. Therefore without going into any details, the final answer after in- 
tegrating by parts, adding appropriate Gibbons-Hawkings terms and then using the 
equation of motion for gW, we have 



'total 



Cr(r, q) - Xr{r, q) $W(r, g)$W(r, -q) 



+ 



5r(r, q) - Arir, q)] M\r, g)<|.W(r, -q) + <|.W(r, g)<|.'W(r, -q) 



+ (E^-FrWj;}{r,q) 



(3.106) 



The values of the coefficients are given in Appendix F of [57]. The above form is 
exactly as we had before, and so all we now need is to get the mode expansion for 
Note however that the subscript m can take only two values, namely m = 0, 1 



as there are only two distinct fields gQQ and g[i 
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mode expansion is: 
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. Therefore our proposed 



(3.107) 



Just like our analysis in section 2.4, the action in f l3.106p is divergent due to terms 
of 0{r'^),0{r^) and hence we need to renormalize the action. The equations for 
renormalization are identical to the set of equations (I2.7ip — fl2.76p . and therefore we 

analogously subtract the counter terms to obtain the following renormalized action: 
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45(4)[/3]$ _ 45(4)[/3]$ , s(5)[/3]$' , s(5)[/3]$' 1 + I V "^^^^^^ I $ 



(3.108) 
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where the radial coordinate is measured at the two boundaries Vh and Tc and $m are 
independent of r as before. Note that X[r(Q,)] is a function independent of and 
appears for generic renormahzed action. 

Now the generic form for the energy momentum tensor is evident from looking at 
the linear terms in the above action fl3.108p . This is again the same as before. However 
now we also need the entropy from the energy- momentum tensor as in fl3.96l) . The 
result for the energy-momentum tensor at r = Tc is given by: 



irpmm\ _ [ ^ 1 ^ \ ^ 



TTmn[l] rjnm[l]. U)[l3] T^m.n[l] 



1=0 '"c(a) 



(3.109) 



The explicit expressions for the coefficients listed above, namely, Hl™^^^\ K'!™^^^\bl^^^-^^j^^ 
and sll}j^^ are given in Appendix F of [57] . 

To complete the story we need the contribution from the UV cap. This is similar 
to our earlier results. The final expression for the ratio of the energy-momentum 
tensor to the temperature takes the simple form: 



(rZd)def _ T^T^hjT) f d^q ^ r/rr'n«[l] , it"'"[1]x (4)[/3] 4/7^™"[l] 

oo 

+ (^H + K7^'^>nF^ + E b^m^nmC-^"'^^] (3.110) 

j=0 

We would like to make a few comments here: First, observe that the final result is 
independent of our choice of cut-off. Secondly, in the string frame there should be a 
l/g^ dependence. Finally, we can pull out a term because the coefficients have 
an explicit dependences (see Appendix F of [57J). This means that both from the 
AdS and the deformed calculations performed above we can show that the entropy is 
of the form: 



l + 0(gsNf,gsM,e-^-^] (3.111) 

where the first part is from fl3.105p and the second part is from (13.1101) . The result 
for the parent cascading theory is (I3.95p . and so we should regard (13.1111) as the 
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entropy for the theory with Af^v degrees of freedom at the boundary. Of course in the 
hmit Muv n » 1 we should recover the entropy formula f l3.95p for the parent 

theory. All in all we see that the correction due to f/uv degrees of freedom only goes 
as e"-'^"", so in practice this is always small for the type of Afuv that we consider here. 
This means that we can use the entropy for the parent cascading theory to estimate 
the viscosity by entropy ratio for a system with degrees of freedom at the UV as: 

1 + Er=i «fce-^'^- 



V 

s 



(3.112) 



471 + i log^ (1 - r^e-^""] ^ 

where we see that the boundary entropy term (13.951) neatly cancels the coefficient 
in the viscosity (I3.9ip to give us the precise bound of when Af^^ — i- oo. Of course 
from our other analysis (13.1 lip we might expect a O (^ggNj, ggM, e~-^""^ contribution 
that would make (I3.112p saturate the celebrated bound ^ if the total entropy density 
factors compensate the factors coming from the viscosity. This would seem consistent 
with, for example, |143j^ . In fact our conjecture would be for non-zero M,Nf and 



J^uv — > e >> 1, the bound is exactly saturatecj^fl i-e 



V _ 1 



in' 



Our second and final step would be to incorporate both the RG fiow as well as 
curvature square corrections. As we discussed before the curvature squared correc- 
tions are typically of the form c^R^upaR'^'^^'^ with cs being the coefficient ( I3.53P that 
we computed before. 

The crucial point here is that (see |146j where this has also been recently empha- 
sized) in the presence of curvature squared corrections the five dimensional metric 
itself changes to: 

ds'^ -«■''•> dt- W'''T;"' dr-+ (3.113) 

'h{r,7i,7i) 92{r) Jh{r,n,n) 



where the black hole factors Qi are no longer given by (I2.32p . They take the following 



^^Provided of course if we assume that a^'s are more general now, being functions of T, gsM, gsNf. 
This way even for non-zero M, Nf , whenever we have Muv oo the bound is exactly . 

^^Note that any possible deviations from ^ due to p.88p in (|3.87|) cannot happen because the 
underlying holographic renormalizability will make di = bi = 0, as discussed earlier. Thus the 
bound in itself is a rather strong result. 



114 Chapter 3 Application to Thermal QCD 



forms: 



gi{r) = l-^ + a + 7^ + a, 



mn _ 



g.{r) = l-^ + « + 7^+^^.^ (3.114) 

where arran,/3mn are all of 0{gsM, ggNf) and can be worked out with some effort (we 
will not derive their explicit forms here). Similarly we could also express f l3.114p in 
terms of inverse powers of r to have good asymptotic behavior. Observe that we can 
still impose gi ~ g2 because the corrections are to 0{gsNf, ggM), although all our 
previous analysis have to be changed in the presence of curvature corrections because 
the explicit values of gi{r) have changed. We will address these issues in our future 
work. Finally {a, 7) are given by 

At this point one might get worried that the metric perturbation on this background 
would become very complicated. On the contrary our analysis becomes rather simple 
once we ignore terms of 0{c3gsM, c^ggNf) (which is a valid approximation with C3 << 
1). In this limit the metric perturbation can be written simply as a /mear combination 
of the terms proportional to C3 in $ which appears in |146j and our solution (13. 75 p . 
(I3.76P derived for RG flow. The final result is: 

4>{r,\oo\)±,R^ = 1 ± ^J^lnr) + [1 + g(r)]log g{r) + + - a + a/M 

+ Q{r) - Y^^{^(0 + ^(^) log 9{r) + [1 + Q{r)]\oig{r)] 

+ 0(|a;n + 0{c^gsNf) + 0{c^gsM) (3.116) 

where we have written 12.321 as; 

gi{r) = g{r) + 0{c:,gsNf) + 0{c^gsM) 

92{r) = g{r) + 0{c:,gsNj) + 0{c^gsM) (3.117) 

with g[r) = 1 — ^ being the usual black hole factor. Of course as emphasized above, 
this is valid only in the limit C3 << 1, which at least for our background seems to be 
the case (see 13.530 . 
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The above corrections are not the only changes. The entropy computed earher also 
gets corrected and therefore the horizon can no longer be at r = r/^. To evaluate the 
correction to entropy we again ignore the terms of 0{c^gsM, csQsNf). In this limit the 
correction terms are precisely given by the analysis of [146] and are proportional to 
the C3 factor fl3.53p as expected. This means that the final result for r]/s including all 
the ingredients i.e RG flows, Riemann square corrections as well as the contributions 
from the UV caps; is given by: 
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s 



1 + Er=i«fce"''-^"^ 
47r+i log^ (1 - r%-^"\ 



C3K 



3^^2(1 _7-4g-4A/;„)3/2 



5,(47r2 - log' Co) + inAo log Co 



(3.118) 



[Att^ - log' Co) + 167r2 log' Co 
where we see two things: one, the bound is completely independent of the cut-off 
r = Tc in the geometry, and two, the bound decreases in the presence of curvature 
square corrections even when A/"™ — )■ e~" with n = (9(1) O The constants appearing 
in 13.1181 are defined as: 



Co = i_rv^"" 

Ao = -isr^e"^-^"" + {sT 

+ 24 (1 + rV^"") log Co 

Bo = -SSnT^e-^^-^ + ASnT^e-^-^ 



47r^e 



+ 48 



-4J\fuv 



log Co + 2Qr^e 



(3.119) 



This is consistent with [146j, and the only violation oi rj/s may be entirely from the 
C3 factor provided the increase in bound from the first term of 03.1181) is negligible, 
as we discussed earlier for fl3.112p . This means in particular: 
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s 



(3.120) 



""^^In [159] ■ non-relativistic systems that appear to have no lower bound were constructed. However, 
these are systems which necessarily require large chemical potentials and low temperature. In highly 
relativistic system created at high energy colliders such as the RHIC or the LHC, chemical potentials 
are small and temperature is high. Our discussion here assumes that the system under discussion 
has such properties so that the use of thermodynamic identity e + P = Ts is valid. Hence, our 
discussion here is in no direct conflict with the models constructed in |159] . 



116 Chapter 3 Application to Thermal QCD 



where Ub can be extracted from (13.1181) . Here we will not study the subsequent 
implication of this result, for example whether there exists a causality violation in 
our theory due to the curvature corrections as in |147] - |154] . We hope to address this 
in our future work. 

3.3 Confinement in QCD 

At low temperatures only color neutral states exist and there are no free color charges 
i.e. quarks in QCD at temperatures below some critical value. The mechanism which 
may explain this confinement of quarks was first proposed by Wilson back in the 
early 70 's |16U] where the only requirement is the existence of abelian or non-abelian 
gauge fields which are strongly coupled to the matter fields that confine. The strong 
coupling calculation can be done using lattice gauge theory and it turns out that at 
low temperatures one can show that the free energy of a couple of charges grow linearly 
with the distance between them |161]|162] . This means it takes infinite amount of 
energy to separate two color charges which are strongly coupled to gauge fields and 
thus only color singlet combinations of charges have finite energy at low temperatures. 
The linear behavior of free energy as a function of inter charge separation is referred as 
linear confinement. At high temperatures the free energy takes the form of Coulomb 
potential |161j |162j and thus the interaction between the charges is negligible for large 
distances. This results in a deconfined phase of matter and gauge fields at sufficiently 
large temperatures. 

The linear confinement of quarks at large separation and low temperatures is a 
strong coupling phenomenon and one uses lattice QCD to compute the free energy 
of the bound state of quarks. On the other hand at high enough temperatures, if the 
gauge coupling is weak, one can compute the free energy using perturbative QCD 
and obtains Coulombic interactions between the quarks. Whether using perturba- 
tive QCD at weak coupling or lattice QCD and effective field theory techniques at 
strong couplings, one finds that at high temperatures, the Coulomb potential is Debye 
screened |163] |ltilj . 
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The study of heavy quark potential gathers special attention as it is linked to a 
possible signal from Quark Gluon Plasma. In particular, heavy quarks are formed 
during early stages of a heavy ion collision as only then there exists enough energy for 
their formation . On the other hand right after the collision QGP is formed as only 
then the temperature is high enough for quarks to be deconfined. As temperature 
goes down, heavy quarks form bound states and as these states are very massive, 
they are formed at temperatures higher than deconfinement value. This means heavy 
quark bound states like J/ip can coexist with QGP and can act as probes to the 
medium. 

In particular one can study the J/ip bound state formed in proton-proton colli- 
sions and compare with heavy ion collision where a medium is formed. The medium 
will screen the cc interaction making them less bound and eventually resulting in a 
suppression of J/ip production. This phenomenon is known as the J/ip suppression 
and considered as a signal of QGP formation [92] . 

In order to quantify quarkonium suppression and analyze features of plasma that 
causes the suppression, one has to compute free energy of J/iJ. At large couplings, 
lattice QCD calculations are most reliable and there has been extensive studies quan- 
tifying the potential for heavy quarks |164] - [T73] . The effective potential between 
the quark anti-quark pairs separated by a distance d at temperature T can then be 
expressed succinctly in terms of the free energy F{d,T), which generically takes the 
following form: 

F{d, T) = ad fs{d, T) - ^Md, T) (3.121) 

d 

where a is the string tension, a is the gauge coupling and fc and fs are the screening 
functionO (see for example |164] - PTD] and references therein). At zero temperature 
free energy is the potential energy of the pair. On the other hand free energy or 
potential energy of quarkonium is related to the Wilson loop which we will elaborate 
here. 



"'^''We expect the screening functions fs,fc to equal identity when the temperature goes to zero. 
This gives the zero temperature Corneh potential 
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Consider the Wilson loop of a rectangular path C with space like width d and time 
like length T. The time like paths can be thought of as world lines of pair of quarks 
QQ separated by a spatial distance d. Studying the expectation value of the Wilson 
loop in the limit T — > oo, one can show that it behaves as 

{WiC)) ~ exp{-TEQQ) (3.122) 

where Eqq is the energy of the QQ pair which we can identify with their potential 
energy YgQ^d) as the quarks are static. At this point we can use the principle of 
holography [23] [SS] |174] and identify the expectation value of the Wilson loop with 
the exponential of the renormalized Nambu-Goto action, 

{W{C)) ~ exp(-^^-) (3-123) 

with the understanding that C is now the boundary of string world sheet. Note that 
we are computing Wilson loop of gauge theory living on flat four dimensional space- 
time Whereas the string world sheet is embedded in curved five-dimensional 
manifold with coordinates and r. We will identify the five- dimensional manifold 
with Region 3 that we discussed in section 2.4. For the correspondence in (13.1231) to 
be valid, we need the t'Hooft coupling which is the gauge coupling in the theory to be 
large. On the other hand as discussed before, it is in this regime of strong coupling 
that linear confinement is realized in gauge theories. Thus using gauge/gravity duality 
is most appropriate in computing quarkonium potential. 

To be consistent with the recipe in [59], we need to make sure that the induced 
four dimensional metric at the boundary of the string world sheet C is fiat. For an 
AdS space, this is guaranteed as long as the world sheet ends on boundary of AdS 
space where the induced four dimensional metric can indeed be written as 1]^^. Using 
the geometry constructed in section 2.4, with metric of the form (I2.9ip and warp 
factor given by (I2.120p for large r, we see that the metric is asymptotically AdS and 
therefore induces a fiat Minkowski metric at the boundary via: 

lira u^g^^ = 7]^y (3.124) 
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where u = and g^y is the full metric (including the warp factor) in Region 3. Thus 
we can make the identification f l3.123p . Once this subtlety is resolved, comparing 
fl3.122p and f l3.123p we can read off the potential 

Cr 
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^ren 



(3.125) 



Thus knowing the renormalized string world sheet action, we can compute Vqq for a 
strongly coupled gauge theory. 

For non-zero temperature, the free energy is related to the Wilson lines W (±f ) 

via: 
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(3.126) 



In terms of Wilson loop, the free energy fl3.12ip is now related to the renormalized 
Nambu-Goto action for the string on a background with a black-holc£|. One may 
also note that the theory we get is a four- dimensional theory compactified on a circle 
in Euclideanised version and not a three-dimensional theory. 



3.3.1 Computing the Nambu-Goto Action: Zero Tempera- 
ture 

Our first attempt to compute the NG action would be to consider the zero temperature 
for the field theory. This means that we take the black hole factors Qi in f l2.9ip to 
be identity. The string configuration we take to compute the action is shown in Fig 
3.2. Note that we are considering the case when the string is exclusively in region 
3 of geometry shown in Fig 2.17 and the justification for this will be given shortly. 
Even if the string enters region 2 and 1, the arguments for linear confinement still 

-•^^There is a big literature on the subject where quark anti-quark potential has been computed 
using various different approaches like pNRQCD [Ill]-[ll3], hard wah AdS/CFT [175] -[179] and 
other techniques |180j — [52] , Its reassuring to note that the results that we get using our newly 
constructed background matches very well with the results presented in the above references. This 
tells us that despite the large N nature there is an underlying universal behavior of the confining 
potential. 
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holds provided the warp factor is monotonic and satisfies a certain equation. We will 
discuss it later in this section. For the time being observe that the configuration in Fig 
3.2 has one distinct advantage over all other configurations studied in the literature, 
namely, that because of the absence of three- forms in region 3, we will not have the 
UV divergence of the Wilson loop due to the logarithmically varying B field [72] . 



B 



C D 




r =co 



A 



-0123 



l"min r - a 



Figure 3.2: The string configuration that we will use to evaluate the Wilson loop in 
the dual gauge theory. The line A determines the actual boundary, with the line B 
denoting the extent of the seven brane. We will assume that line B is very close to 
the line A. The line C at r = denotes the boundary between Region 3 and Region 
2. Region 2 is the interpolating region that ends at r = r^[^. At the far IR the 
geometry is cut-off at r = a from the blown- up S^. 



As the system is not dynamical, the world line for the static QQ can be chosen to be 

a;2 = = (3.127) 
and using u = 1/r we can rewrite the metric in region 3 aa^^l: 



X 



d 

±2' 



—g{u)dt^ + dut"^ 
1 



(3.128) 



^^We use the Einstein summation convention henceforth unless mentioned otherwise. 
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where An are the coefficients that can be extracted from the in fl2.120p . the black 
hole factor g{u) = 1 for the zero-temperature case, and ds'j^^ is the metric of the 
internal space that includes the corrections given in f l2.117p . This can be made precise 
as 



giving ^0 = f^)-^! = —ji7,A.2 = ji: — and so on. Observe that since a^, 
i > 1 are of 0{gsNf) and oc yfgj^, all Ai are very small. The r"" corrections along 
the radial direction given in (12.1 ISp are accommodated above through Biv} series. 

Now suppose : (cr, r) — )■ u, ^z^, 0j, 6'j) is a mapping from string world 

sheet to space-time. Choosing a parametrization r = x° = t, a = a;^ = x with the 
boundary of the world sheet overlapping with the world line of the QQ pair, we see 
that we can have 

X° = t, = X, = = 0, = X^ = ^ = 

(X^X^) = (^1,00 = (7r/2,0), (X^X9) = (^2,02) = (vr/2,0) (3.130) 



a^u / 3af ao . 



(3.129) 



which is almost like the slice that we chose in ^57j. The advantage of such a choice is 
to get rid of the awkward angular variables that appear for our background geometry 
so that we will have only a r (or u) dependent background. We also impose the 
boundary condition 

u{±d/2) = ^ (3.131) 

where u-^ denote the position of the seven brane closest to the boundary u = 0. The 
Nambu-Goto action for the string connected to this seven brane is: 

5'stnng = ^ / d^d-^ [v^-det [{g^, + 9^09,0) S^Xa'S.X-] + \e'''B,, + J(0) 

+ e»^9aX™9,x" 6 r^r^^^-T, e Fa^,,,,, + o(e^)] (3.132) 

where a, 6 = 1, 2, di = 82 = The other fields appearing in the action are the 
pull backs of the NS B field Bab, the dilaton coupling J(0) and the RR field strengths 
Fabc- Its clear that if we switch off the fermions i.e = = the RR fields 
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decouple. The B^s field do couple to the fundamental string but as we discussed 
before, in region 3 we don't expect to see any three-form field strengths. This is 
because the amount of i?NS that could leak out from region 2 to region 3 is: 

Sns = MS[l-f{r)] = M5e-"(^-''°\ r > tq (3.133) 

where S is the two-form: 

S^9s (h{r) cot ^ dOi + &2(r) cot ^ rf^2 ) A e^, - sin 02 d02 A #2 (3.134) 



An 



^ + 9sNf-^ + log sm ^ sm I + 63 (r) 



sin ^1 dOi A d0i 



and 6n have been defined in section 2.4.2. We see that not only 5ns has an inverse r 
fall off, but also has a strong exponential decay as a >> 1. This is the main reason 
why there are no NS or RR three- forms in region 3, making our computation of the 
Wilson loop relatively easier compared to the pure Klebanov-Strassler model. 

On the other hand the dilaton will couple additionally via the J(0) term. Although 
this coupling of 4> is not directly to the X'^, we can still control this coupling by 
arranging the other seven-branes such that: 

E^-E^ <e forO<n<m„ (3.135) 

with e very small and nio a sufficiently big number. Under this condition the dilaton 
will be essentially constant and the axio-dilaton r would behave as: 

- = -o + E^ + iE^ (3.136) 

n=l n>mo 

SO that its contribution to NG action can be ignored although the Biu^ contribution 
still remains, because the seven-branes continue to affect the geometry from their 
energy-momentum tensors and the axion charges. In this limit both string and Ein- 
stein frame metrics are identical and the background dilaton is 

= log gs - + 0{gf) (3.137) 
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which, in the hmit Qs 0, will be dominated by the constant term (note that ruo is 
fixed). On the other hand, ~ 4'R{2), where i?(2) is the Ricci scalar for the world 
sheet metric. As the background space time metric run as (9(l/r"'), we can make i?(2) 
arbitrarily small by a reparametrization of the world sheet. This means we can ignore 
J(0) in the string action while the NG string will see a slightly different background 
metric as evident from fl3.132p . 

Thus once all the effects have been accounted for, using the metric f l3.128p with 
the embedding X^^ given by fl3.130p . one can easily show that at zero temperature 
the NG action is given by: 

I dx 



NG 



T 
2^ 



A 



du 
dx 



(3.138) 



where we have used J dt = T,T = T/a', Vn+mo = {n + mo)T>n+mo'i and An,Bn and 
T^n+nio are now defined for choices of the angular coordinates given in (I3.130p . The 
above action can be condensed by redefining: 



(3.139) 



which would mean that the constraint equation i.e diTl = 0, being the stress- 
tensor, for u{x) derived from the action fl3.138p using fl3.139p can be written as 

f n 2 \ 

"""^ = (3.140) 



d_ 

dx 



implying that: 



Co 



(3.141) 



where Co is a constant, and u'{x) = This constant Co can be determined in the 
following way: as we have the endpoints of the string at x = ±d/2, by symmetry the 
string will be U shaped and if Umax is the maximum value of u, we can choose the x 
coordinate such that u(0) = Umax and u'{x = 0) = 0. Plugging this in (13. 141 p we get: 



Co 



•^n "max 
ul„„ 



(3.142) 
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Once we have Co, we can use fl3.140p to get the following simple differential equation: 

1/2 



du ^ 1 
dx CoVO^- 



-ci{A.,u^y 



(3.143) 



which in turn can be used to write x{u] as: 

Co 



X[U] 



dw 



(3.144) 



where we have used x(ur 



0. Now using the boundary condition given in f l3.13ip 



i.e x{u = Uy) = d/2, and defining w 



Ur 



<v,eo 



d 



dv V 



~lAJcJf^ 



— we have 

max 

A 7/" 



V 



A n 



m n,m 



-1/2 



(3.145) 



At this stage we can assume all An > 0. This is because for An > we can clearly 
have degrees of freedom in the gauge theory growing towards UV, which is an expected 
property of models with RG flows. Of course this is done to simplify the subsequent 
analysis. Keeping An arbitrary will also allow us to derive the linear confinement 
behavior, but this case will require a more careful analysis. We will get back to this 
in the section 3.3.3. Note also that similar behavior is seen for the the Klebanov- 
Strassler model, and we have already discussed how degrees of freedom run in regions 
2 and 3. Another obvious condition is that d, which is the distance between the 
quarks, cannot be imaginary. From f l3.145p we can see that the integral becomes 
complex for 

- ^^f Y > 1 (3.146) 



whereas for J-'{v) = 1 the integral becomes singular. Then for d to be always real we 
must have 



J^{v) < 1 



(3.147) 



We will now set, without loss of generality, = 1 and = . Such a choice is 
of course consistent with supergravity solution for our background (as evident from 
fl3.129p ). This choice also defines our units in that we are setting the AdS throat 
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radius L = 1 and this is only to make subsequent computations convenient. . Now 
analyzing the condition f l3.147p . one easily finds that we must have 



1 



-^m + 1)^^+3 <if < 1 



(3.148) 



for d to be real. This condition puts an upper bound on Umax and we can use this 
to constrain the fundamental string to lie completely in region 3. Observe that for 
AdS spaces, An = for n > and hence there is no upper bound for Wmax- This is 
also the main reason why we see confinement using our background but not from the 
AdS backgrounds. Furthermore one might mistakenly think that generic Klebanov- 
Strassler background should show confinement because the space is physically cut-off 
due to the presence of a blown- up S^. Although such a scenario may imply an 
upper bound for Umax, this doesn't naturally lead to confinement because due to the 
presence of logarithmically varying 5ns fields there are UV divergences of the Wilson 
loop. These divergences cannot be removed by simple regularization schemes [72] and 
creates problems in the interpretation of Wilson loop which describe the potential. 

Coming back to (13.1381) we see that it can be further simplified. Using (I3.14ip . 
(13.1421) and (13.1431) in (I3.138p . we can write it as an integral over u: 



5 



NG 



T 

TT 



T 



1 dv 



-1/2 



^ .Aim "max ^ , 



-1/2 



(3.149) 



where in the second equality we have taken v = tt/wmax- 

This simplified action (I3.149P is however not the full story. It is also divergent in 
the limit Cq — 0. But we can isolate the divergent part of the above integral (13.1490 
by first computing it as a function of eo- The result is 



S- 



NG 



+ 



^NG 
f 1 



NG 



T 1 



TC U„ 



TT Ur, 



1 dv 

ymU 



m"'max'^ 



A 11^' 



A-m.^ 



m n,m 



-1/2 



- 1 ks.iso) 
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Now by expanding y QmUmax^^ = we can compute the first integral to be 

S],o = l-^(-^o + ^ + T.Al + <^i^o) + .) (3.151) 
where Qq = Qq, Qi = ^Giu^ax and so on. The second integral becomes 

-1/2 \ 

- 1 + 0{el) 

(3.152) 

where the eo dependence here appears to C(ej^); and we have set Qi = without loss 
of generality. Now combining the result in fl3.15ip and f l3.152p . we can obtain the 
renormalized action by subtracting the divergent term 0{l/e) in the limit e^, — )■ and 
we obtain the following result 

^+0(^?2) 

+ 0{eo)\ (3.153) 

where the third term in (13.1531) . including the 0{g1) correction, is related to the action 
for a straight string in this background in the limit ^ Q. Our subtraction scheme 
is more involved because the straight string sees a complicated metric due to the 
background dilaton and non-Ricci fiat unwarped metric. This effect is independent 
of any choice of the warp factor and we expect this action to be finite in the limit 

Once we have the finite action, we should use this to compute the QQ potential 
through fl3.125p . Looking at fl3.145p we observe that the relation between d and Umax 
is parametric and can be quite involved depending on the coefficients An of the warp 
factor. If we have An = 0, Qn = for n > 0, we recover the well known AdS result, 
namely: d ~ Mmax and Vqq ~ ^. But in general fl3.145p and f l3.153p should be solved 
together to obtain the potential. 

As it stands, fl3.145p and fl3.153p are both rather involved. So to find some corre- 
lation between them, we will observe the limiting behavior of Wmax. Therefore in the 



oil 



T 1 



Vr Mr, 



A v"' 



A. 



cren 
^NG 



+ 



T 1 

^ ^ma; 

1 dv 
f ^ 



1=2 



I - 1 



1-t;^ 



1 dv 
v'^ 

AnU 

. AmU 



,m 



m n.m 



m^-max' 
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following, we will study the behavior of d and S^q for the cases where Wmax is large 
and small. 



Quark-Anti quark potential for small Wmax 

Let us first consider the case where "Umax is small. In this limit we can ignore higher 
order terms in Mmax and approximate 



^0 + •^2U 



2 

max 



1 + rj 



(3.154) 



where = 1 and ^2iimax — V- Using this we can write both f l3.145p and fl3.153p as 
Taylor series in rj around 1] = 0. The result is 



d = ^/rj ao + aiT] + 0{r]'^ 
bo + biT] + 0{ri'^ 



Qren 



T 



(3.155) 



with oo, ai, 6o; are defined in the following way: 



ao 
ai 

bo 
bi 



dv- 



/A2J0 



dv 



V 



1 / 



-1 + 



dv 



VA2J0 



1.1981 



l-V 



4A, 



2A2 



1-VT 







O.62JA2 



t;2(l + 1,2)^/13^ 



(3.156) 



where we have taken Qq = 1 and Qi = without loss of generality. In this limit clearly 
increasing 77 increases d, the distance between the quarks. For small t], d = aQ^yrj, 
and therefore the Nambu-Goto action will become: 

/ ao\bo\\ 1/61 



oren 
'-'NG 



T 



1 

d + 



d + Oid' 



(3.157) 



TT /a yTfCio, 

where all the constants have been defined in f l3.156p . Using f l3.125p we can determine 
the short- distance potential to be : 



Vr 



QQ 



1 

a' 



9sN 



' ao\bo\ \ 1 
^ Ti ) d 
0.236 



d + 0{d^) 



+ (0.174^2 + 0.095^2) + 0{d^ 



(3.158) 
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where N is the number of D3 branes in the gauge theory and we have used string 
tension 1/a' = ^/gJ^ / L"^ = \JgsN as = 1 according to our choice of units. 
The potential is dominated by the inverse d behavior, i.e the expected Coulombic 
behavior. Note that the coefficient of the Coulomb term which is a dimensionless 
number only depends on the number of D3 branes. Thus it is independent of the 
warp factor and hence should be universal. This result, in appropriate units, is of 
the same order of magnitude as the real Coulombic term obtained by comparing with 
Charmonium spectra as first modelled by |184j and subsequently by several authors 
|164] - [T7D] . |171] - 1173] . |176] - 1179] . This prediction, along with the overall minus sign, 
should be regarded as a success of our model (see also |185] where somewhat similar 
results have been derived in a string theory inspired model). The second term on the 
other hand is model dependent, and vanishes in the pure AdS background. 

Note also that the above computations are valid for infinitely massive quark-anti 
quark pair. For lighter quarks, we expect the results to differ. It would be interesting 
to compare these results with the ones where quarks are much lighter. 

Quark-Anti quark potential for large Wmax 

Lets analyze the integrals (13.1 45 p and (I3.153P in the limit Wmax is close to its upper 
bound set by (I3.148P (see also |185j ). In particular if Umax is the upper bound of Umax, 
then it is found by solving 



We observe that both the integrals (I3.145P and (13.1531) are dominated by w ~ 1 
behavior of the integrands. Near v = 1 and Umax ~^ Umax the distance d between the 
quark and the anti quark can be written as: 



(m + l)^^+3u; 



m+3 



(3.159) 



max 



d 




■max 



■max 



L 



log A - log (2^B(A + B) + 2B + a) 



y'A(l -f) + B(1 -u)2 



dv 



Vb 



(3.160) 
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where note that we have taken the lower hmit to 0. This will not change any of our 
conclusion as we would soon see. On the other hand, the renormalized Nambu-Goto 
action for the string now becomes: 

dv 



Qren 



TT 



1 



Afl 



B(l-^;)2 



T 



7ru„ 



+ 0{VL 



maxV 



f \ GmK 



TT 



f 



logA - log (2y^B(A + B) + 2B + A 



- 1 



+ 0{u 



' 1 

max./ 



where A and B are defined as: 



4-2- 



■^m "max 



(3.161^ 



(3.162) 



B 



j ^-^""^max 
•^rn "max 



^■^""^max 



(n^ — n)AnU 



n 

max 



'^m^'max 



6 



Observe that in the integral f l3.160p and f l3.16ip we have to take the limit Ur 
So A, B should be evaluated in the same limit. Interestingly, comparing (13.1621) to 
(13.159P we see that 



lim 

Umax Climax 



thus vanishes when computed exactly at u 
at that point and in fact behaves as: 



A^O (3.163) 
The other quantity B remains finite 



B 



'^rn "-max 



(3.164) 



Our above computation would mean that the distance d between the quark and the 
anti quark, and the Nambu-Goto action will have the following dominant behavior: 



d 



oren 



2A/^mU™ax^max log € 



A ii" 



lim 



T v ^mii^ax log e 
lim — j=- 

£^0 TT Umax a/B 



B 



(3.165) 



which means both of them have identical logarithmic divergences. Thus the finite 
quantity is the ratio between the two terms in 13.1651 This gives us: 



Qren 
^NG 

d 



± wA^n^^max 



T X constant 



(3.166) 
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Now using the identity f l3.125p and the above relation fl3.166p we get our final result: 

Vqq = i^^) d/a' (3.167) 

\ max / 

which is the required linear potential between the quark and the anti quark. 

Before we end this section one comment is in order. The result for linear con- 
finement only depends on the existence of u^ax which comes from the constraint 
equation fl3.159p . We have constructed the background such that Umax lies in region 
3, although a more generic case is essentially doable albeit technically challenging 
without necessarily revealing new physics. For example when u~^^^ is equal to the 
size of the blown up at the IR will require us to consider a Wilson loop that goes 
all the way to region 1. The analysis remains similar to what we did before except 
that in regions 2 and 1 we have to additionally consider fields of the form 
and log u respectively. Of course both the metric and the dilaton will also have non- 
trivial M-dependences in these regions. One good thing however is that the Wilson 
loop computation have no UV or IR divergences whatsoever despite the fact that now 
the analysis is technically more challenging. Our expectation would be to get similar 
linear behavior as fl3.167p here too. 



3.3.2 Computing the Nambu-Goto Action: Non-Zero Tem- 
perature 

After studying the zero temperature behavior we will now discuss the case when 
we switch on a non-zero temperature i.e make g{u) < 1 or equivalently the inverse 
horizon radius, Uh finite in fl3.128p . where 



9[u) 



u 



uj 



(3.168) 



Choosing the same quark world line (13.1270 and the string embedding fl3.130p with 
the same boundary condition (I3.13ip but now in Euclidean space with compact time 
direction, the string action at finite temperature can be written as 



5'ng — — I ^ 



27r 



^^^(m)(Aw"J + 



' du\' 
dx I 



(3.169) 
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where QmU"^ is defined in (13.1391) and the correction to QmU^ is suppressed by g1 
as well as u^/u\ because the background dilaton and non-zero temperature induces 
a slightly different world-sheet metric than what one would have naively taken. To 



avoid clutter, we will further redefine these corrections as: 



4:+n+l+k+2mo 



(3.170) 



Minimizing this action gives the equation of motion for u{x) and using the exact 
same procedure as for zero temperature, the corresponding equation for the distance 
between the quarks can be written as: 



d 



1 _ "max A „.n 
h 



1-^ 



'1/2 



1 



^ .Aim "max ^ 



(3.171) 



Once we have d, the renormalized Nambu-Goto action can also be written following 
similar procedure. The result is 
T 1 



Qren 



71 U„ 



j - ^0 + E 7^ - /' + 0{gl) (3.172) 

L 1=2 I- i- JO V 
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^ dv 





1-v' 



1 ^"^^iax 

i a — 



A 7,m ,,m 
^-^mU'max'^ , 



-1/2 



which is somewhat similar in form with (I3.153p . which we reproduce in the limit 
■u/i — i- oo. Also as in (I3.153p . we have defined '\J i'mU^ax^"^ = ^z"^'- 

Now just like the zero temperature case, requiring that d be real, sets an upper 
bound to Mmax, that we denote again by Umax, and is found by solving the following 
equation: 



1 



-(m + l)^m+3urx + -n ^-o 
^ J- k=0 ^ 



u: 



J r 



ul 



+ 2j - 1 



1 



(3.173) 



Once we fix Uh and the coefficients of the warp factor A^-, Umax will be known. We 
will assume that Umax lies in region 3. 
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Rest of the analysis is very similar to the zero temperature case, although the 
final conclusions would be quite different. We proceed further by defining certain 
new variables in the following way: 



A, 



EA-n 



1 







[ i-m \ I n 



k-i 

2 



/ - m > 4 

/ - m < 4 
I — m = 



(3.174) 



As before, we observe that for Wmax — ^ Umax, both the integrals (I3.17ip . fl3.172p are 
dominated by the behavior of the integrand near f ~ 1, where we can write 

dv 



d 



m "-max^max 



1 ^max A 

J- ::^>*-nU 



" max 



-2 



'-^m '-'max "-'max 
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logA - log 2VB(A + B) + 2B + A 



(3.175) 



where taking the lower limit of the integral to again do not change any of our 
conclusion. On the other hand, the renormalized Nambu-Goto action for the string 
now becomes: 

/■I dv 
^0 \lMl~v) + Bil-vY 
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-t \/ ^m"r 
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TT 
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max/ 
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vru„ 



(3.176) 



where A and B are defined exactly as in fl3.162p but with An replaced by An given 
by fl3.174p above. It is also clear that: 



lim 



A ^ 



(3.177) 

and so both (13.1750 as well as (I3.176P have identical logarithmic divergences. This 
would imply that the finite quantity is the ratio between (I3.176P and fl3.175p : 
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Now using the identity f l3.126p and the above relation fl3.166p we get our final result: 



\ 



1 - """'"^ d/a' (3.179) 



This gives linear potential at large distances and from the form of the potential above, 
we can see that the higher the temperature, i.e. lower the Uh, the lower the slope of 
the linear term. In fact we will explicitly plot the potential at various temperature 
for certain choice of warp factors in section 3.3.4 and analyze this melting of the 
potential. 

3.3.3 Linear Confinement from generic dual geometries 

Having argued for linear confinement for geometries which are like region 3 of Fig 
2.17, we will now consider a more general warp factor choice of h in dual geometry 
with metric f l2.9ip . Much of the discussion here closely follows our recent work |186] 
and further details can be found there. We will restrict to cascading gauge theories 
where the effective number of colors grows as scale grows. This property of a gauge 
theory is most relevant for physical theories as new degrees of freedom emerge at UV 
and effective degrees of freedom shrink in IR to form condensates at low energy [46j . 
The number of colors at any scale m = 1/r is given by ( I3.33P and for the analysis 
given here, it is simpler to define 



instead of A^'eff(M). The coefficients An in the previous section are related to l-Liu) by 
l-iiu) = AnU"", and h is the warp factor. In terms of ?/(«), the condition that Nes{u) 
is a decreasing function of m = 1/r becomes 

n\u) > (3.181) 

Combining Eqs. fl3.180p and fl3.18ip yields the following condition 

niu) > ^ (3.182) 

From now on, the value of L is set to 1 for the rest of this subsection, so that T-L{u) > 1. 
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Zero temperature 

Let Mmax be the maximum value of u for the string between the quark and the anti- 
quark. Then the relationship between "Umax and the distance between the quark and 
the anti-quark is given by [62] 



2'Urna,x'^ ( 



1 V 

dv 



QmU 



m^-rnax' 



CO 



('H(Mmax^))^ 



T-i{Ur, 



'1/2 



(3.183) 



For the above expression to represent the physical distance between a quark and 
an anti-quark in vacuum, the integral must be real. This is guaranteed if for all 
< < 1: 



< 1 



(3.184) 



For AdS space, fl3.184p is automatic, as "H = 1 and then d is proportional to Wmax 
which results in only Coulomb potential. But for a generic warp factor, (I3.184p gives 
rise to an upper bound for Wmax as already discussed. 

To show confinement at large distances the potential between the quark and the 
anti-quark must be long ranged. That is, diu^aa.^ must range from to oo as Wmax 
varies from to its upper bound, say Wmax = a;max- Since l-i{u) > 1, the only 
way to satisfy these conditions is via sufficiently fast vanishing of the square-root in 

Eq. fl3.183P as f 1 at Umax = a^max- 

For most Mmax, 1 — W iyluuiay)"^ vauishcs only linearly as v approaches 1. In this 
case, (i(Mjnax) is finite as the singularity in the integrand behaves like l/Vl — v and 
hence it is integrable. To make (i(Mmax) diverge at Wmax = a^max, 1 — W^('y|a;max)^ must 
vanish quadratically as v approaches 1 to make the integrand sufficiently singular, 

;|. Therefore, the function iy(f|xii 
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maximum at f = 1. 
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n'jxn 
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(3.185) 
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For this to vanish at f = 1, Xmax must be the smallest positive solution of 

xn'ix) - 21-i{x) = (3.186) 

With the definition Tilu) = AnU"", one can easily show that this is equivalent to the 
condition f l3.159p which was originally derived in |62j. The allowed range of Wmax is 
then 

< Umax < a^max (3.187) 

and within this range, d{uma.x} varies from to oo. How it varies will depend on the 
values of Qm as well as l-i{u). 

Finite temperature 

At finite temperature, the relation between Wmax and the distance between the quark 



and the anti-quark is obtained by replacing 1-L{u) with y 1 — u^/u\l-i{u) in Eq.( l3.183|l : 

rfT(Mmax) = 2Mmax\/l " U^^J ulH{u^s^^) f dv — 

Jen i 



X 



^ _ ^4 (1 - <ax/4) f ^(Wmax) ^ ^ 



1/2 

(3.188) 



;i - ^^<ax/0 \^{Um^^v)J 

The explicit factor of "Umax makes ^^(Mmax) vanish at Wmax = as in the T = case. 
As Umax approaches Uh, the integral near v = 1 behaves like 



dAu^..) ~ / dv ^ (3.189) 

y^(l - U)(l - WUmax/Uh) 

which indicates that (iT(wmax) goes to as Umax approaches Uh. Hence, at both 
Mmax = and Umax = Uh, dxiumax) vauishcs. Siuce dxiumax) is positive in general, 
there has to be a maximum between Umax = and Umax = Uh- Whether the maxi- 
mum value of (iT(umax) is infinite as in the T = case depends on the temperature 
(equivalently, uj^'^) as we now show. 

The fact that the physical distance needs to be real yields the following condition. 
For all < u < 1, 



Wr(u|Umax) = 



K(Ur 



1 '^max/'^fe 



< 1 (3.190) 
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Taking the derivative gives 



X 



(3.191^ 



Similarly to the T = case, (irlwrnax) can have an infinite range if the derivative 
vanishes at f = 1 for a certain value of Mmax, say Mmax = Z/max- This value of ?/max is 
determined by the smallest positive solution of the following equation 

yU\y) - 2U{y) = {yjuuf yn\y) (3.192) 

which then forces W^(l|?/max) to vanish. Note that the left hand side is the same as 
the zero temperature condition, Eg. (13. 1861) . The right hand side is the temperature 
{uh) dependent part. Using the facts that: 

i-i oo 

n (A: - 1/2) = -(2j - 3)!!/2^ ^a;^' (2j - 3)!!/2^j! = 1 - v^T^, (3.193) 
fc=o i=i 

it can be readily shown that Eg. (13. 192]) is eguivalent to Eg. (13. 173]) as long as Mmax < 

Uh- It is also clear that y = cannot be a solution of Eg. (l3.192p because at y = u^, 

the eguation reduces to ^{uh) = which is inconsistent with the fact that Hiy) > 1. 

Recall that we are considering gauge theories for which "H^/) > 1 and 'H'(y) > 0, 
and we assume that the eguation yT-L'^y) — 2'H{y) = has a real positive solution Xmax 
which gives confinement at zero temperature. Hence as y increases from towards 
ajjnax, the left hand side of Eg. (l3.192p increases from —2 while the right hand side 
increases from 0. The left hand side reaches when y = x^ax which is the point 
where the distance (i(umax) at T = becomes infinite. At this point the right hand 
side of Eg. (l3.192p is positive and has the value {xraax/uhYxmaxH,' {xj^ax) ■ Hence the 
solution of Eg. (13. 192 p . if it exists, must be larger than Xmax- 

Consider first low enough temperatures so that Uh ^ a; max- For these low tem- 
peratures. Eg. (13. 192]) will have a solution, as the right hand side will be still small 
around y = Xmax- This then implies that the linear potential at low temperature will 
have an infinite range if the zero temperature potential has an infinite range. 
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Now we show that the infinite range potential cannot be maintained at all tem- 
peratures. We can have a black hole such that Uh = Xmax- When the left hand 
side vanishes at y = Xmax, the right hand side is Xmax'H'(xmax) = 2'H(xmax) which is 
positive and finite. For y > Xmax, the left hand side (yTi'^y) — 2'H{y)) may become 
positive, but it is always smaller than yTi'^y) since T-Li^y) is always positive. But for 
the same y, the right hand side {{y I UhYyl-i' {y)) is always positive and necessarily 
larger than yTi'^y) since (y/uh) > 1. Hence, Eq.f l3.192|) cannot have a real and posi- 
tive solution when Uh = Xmax- Therefore between Uh = oo and Uh = Xmax, there must 
be a point when Eq. fl3.192p cease to have a positive solution. 

When Eq.f l3.192|) has no solution, then the expression for driumax), (I3.188P will 
not diverge for any Mmax within (0, Uh)- Furthermore, since the expression vanishes at 
both ends, there must be a maximum c?j^(Mmax) at a non-zero Wmax- When the distance 
between the quark and the anti-quark is greater than this maximum distance, there 
can no longer be a string connecting the quark and the anti-quark. 

To summarize, we have just shown that if we start with a dual geometry that 
allows infinite range linear potential at zero temperature, there exists some critical 
temperature above which the string connecting the quarks breaks. This shows that 
at high enough temperatures quarkonium state melts and gives rise to 'free quarks'. 
In the following subsection, we will try to quantify the melting temperature using 
geometries with exponential warp factors. 



3.3.4 Numerical analysis of melting temperatures 

After discussing the most general choice for warp factors that give rise to ?/max and 
consequently linear potential, we will now give specific examples of geometries that 
may arise as solutions to Einstein's equation. We start with the following ansatz for 
the metric: 

ds'^ = — ^dt"^ H — -^{dx^ + dy"^ + dz^) + { — -du^ + ds]. 




(3.194) 



138 Chapter 3 Application to Thermal QCD 



where h = h{u, 9i, (pi, ip),H = H{u, 9i, (pi, ip), and g = 1 — u'^/uf; Ai^ is the compact 
five dimensional manifold parametrized by coordinates {6i, (pi, ip) and can be thought 
of as a perturbation over T^'^. Here m = is the boundary and m = m/i is the horizon. 
As discussed in [62], the above metric arises in region 3 of [62] when one considers the 
running of axio-dilaton r, D7 brane local action and fiuxes due to anti five-branes on a 
geometry that deviates from the IR OKS-BH geometry from the back reactions of the 
above sources. The three-form fiuxes sourced by (p, q) anti-branes are proportional 
to r^'^fijr) for some positive i (see [62] for details about f{r)), where the function 
/(r) — 7- 1 as r — )■ oo and /(r) — )■ as r — )■ 0. With the coordinate u = 1/r, there is 
another function: k(u) = exp(— m-^),^ > 0, that also has somewhat similar behavior 
as f{u) and may allow us to have a better analytic control on the background. With 
such a choice of k{u), the total three form fiux is proportional to w^Miu) with 

M{u) = M[l - k{u)] = M [l - exp(-u-^)l (3.195) 

where M is the number of bi-fundamental fiavors. Thus three-form fiuxes are decay- 
ing fast as Mu^ 1 — exp {^—u-^^ and, as shown in [52], the seven-branes could be 
arranged such that the axio-dilaton r behaves typically as r ~ u'^. This means that 
from the behavior of the internal Riemann tensor one may conclude that the internal 
metric (jmn behaves as (jmn ~ u'"exp(coM'') where A,C,A and C are all positive and 
Co could be positive or negative depending on the precise background informations. 

From the above discussions it should be clear that taking the three-forms and 
world-volume gauge fiuxes to be exponentially decaying in the IR (but axio-dilaton 
to be suppressed only as u^) should solve all the equations of motion, giving the 
following behavior for the warp factor h and the internal metric H in fl3.194p r1 : 

h = L\%xp{-av^), H = exp{/3u^) (3.196) 

where we are taking a, a, /3, /3 to be all positives with a, /3 to be functions of internal 
coordinates {6i,(pi,ip) and = QgNa'"^ to be the asymptotic AdS throat radiucl. 



^"See also the interesting works of |187| where exponential warp factors have been chosen. 
^^Note that /? in p.l96|) could be considered negative so that H would be decaying to zero in 
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Motivated by the above arguments, we will consider Nambu-Goto action of the 
string in the geometry with {a, j3) = (3,3) and {a, (3) = (4,4) at temperatures T^^^ 
and T^"^^ respectively in Eq.f l3.196|) . As in [57] [62] we consider mappings X^((T, r), 
which are points in the internal space, to lie on the slice: 



0, 







(3.197) 



so that on this slice a, (3 are fixed and we set it to (a, (3) = (0.1, 0.05) for both choices 
(a,/3). (Such a choice of slice will also help us to ignore the three-form contributions 
to the Wilson loop.) With these fixed choices for the warp factors, we plot the inter 
quark separation as a function of Umax in Figures 3.3 and 3.4 for various values of 
T = IK. 



T=0.01 

T=0.2 

T=0.25 

T=0.286 

T=0.33 

T=0.357 



1.5 



2.5 



3.5 



Figure 3.3: Inter quark distance as a function of Umax for various temperatures and 
warp factor with {a, a, /3, j3) = (0.1, 3, 0.05, 3) in the warp factor equation. 



Note that for both choices of warp factors, for low enough temperatures, there 
exist Umax = Umax wherc d — 7- oo. As the temperature is increased, ?/max increases 
modestly. On the other hand from figure 3.3, one sees that when T > TJ^^ ~ 0.28 
there exists a (imax which is finite. This means for inter quark distance d > (imax, there 



the IR. However since region 3 doesn't extend to the IR we don't have to worry about the far IR 
behavior of Eq. (|3J96| . 
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Figure 3.4: Quark-anti quark distance as a function of Umax for various temperatures 
and warp factor with {a, a, /3, /3) = (0.1, 4, 0.05, 4) in the warp factor equation. 



is no string configuration with boundary condition x(0) = ±d/2 implying that the 
string attaching the quarks breaks and we have two free partons for d > (imax- Thus 
we can interpret (imax to be a "screening length". From Fig 3.4 we observe similar 
behavior but now dmax exists for T > T^^^ ~ 0.399. 

In figure 3.5, (imax as a function of T is plotted. We note that for a small change 
in the temperature near T^^^ (or near Tj^-* equivalently) there is a sharp decrease in 
screening length (imax, but for T >> Tj*\ i = 1,2, the screening length does not 
change much. In fact (imax behaves as C + exp(— 7T) (where C and 7 are constants) 
which in turn could be an indicative of a phase transition near 

rp(i) for i = 1,2 i.e the 

two choices of warp factor. 

Finally we plot the potential energy Vqq as a function of d in Figures 3.6 and 3.7 
for the two choices of warp factor. For T < in Fig 3.6 and T < in Fig 3.7, 
we have energies linearly increasing with an arbitrarily large increment of the inter 
quark separations. Thus we have linear confinement of quarks for large distances 
and small enough temperatures. For T > TJ*\ i = 1 or 2, there exists a (imax and 
for all distances d > (imax there are no Nambu-Goto actions, S'ng, for the string 
attaching both the quarks. This means that we have free quarks and Vqq is constant 
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Figure 3.5: Maximum inter quark separation dmax as a function oiT = l/uu for both 
cubic and quartic warp factors. 



for d > (iinax- Of course looking at Fig 3.6 and 3.7 one shouldn't conclude that the 
free energy stops abruptly. What happens for those two cases is that the string joining 
the quarks breaks, and then the free energy is given by the sum of the energies of 
the two strings (from the tips of the seven-branes to the black-hole horizon) and the 
total energies of the small fluctuations on the world-volume of the strings. The latter 
contributions are non-trivial to compute and we will not address these in any more 
detail here, but energy conservation should tell us how to extrapolate the curves in 
Fig 3.6 and 3.7, beyond the points where the string breaks, for all T > T^. Of course 
after sufficiently long time the two strings would dissipate their energies associated 
with their world-volume fluctuations and settle down to their lowest energy states. 

To compare with lattice QCD calculation of free energy |188j . in Fig 3.8, we plot 
side by side the lattice results and our calculation with = (0.1,3,0.05,3) 

where the potential energy has been extrapolated to account for the energy of the 
disjoint strings i.e. the free quarks for d > dmax- We observe the striking similarity 
between the shape of curves in the two plots. Two completely different approaches 
yield very similar results which only strengthens the validity of applying gauge/gravity 
correspondence in the study of strongly coupled QCD. Qualitatively the curves show 
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Figure 3.6: Heavy quark potential Vqq as a function of quark separation d with cubic 
warp factor, or equivalently, (a, a, /3, (3) = (0.1, 3, 0.05, 3) in the warp factor equation 
for various temperatures. 



how hnear potential melts at high temperatures and there is also evidence of screening 
of the Coulomb potential at small separation of the quarks. 

Observe that for a wide range of temperatures < T < TJ*\ the potential and thus 
the free energy hardly changes. But near a narrow range of temperatures T^*-* — e < 
T < TP + e (where e ~ 0.05), free energy changes significantly. For Fig 3.6 the 
change is more abrupt than Fig 3.7. This means as we go for bigger values of a, the 
change in free energy is sharper. 

In Fig 3.9 , we plot the slope of the linear potential as a function of T. Again 
for a wide range of temperatures, there are no significant changes in the slope but 
near Tj*\ the change is more dramatic: the slope decreases sharply, indicating again 
the possibility of a phase transition near TJ*^. As we noticed before, here too bigger 
exponent a gives a sharper decline in the slope hinting that when a » 1, the 
transition would be more manifest. 

To conclude, the above numerical analysis suggest the presence of a deconfinement 
transition, where for a narrow range of temperatures 0.28 < T^ < 0.39 the free 
energy of QQ pair shows a sharp decline. Interestingly, changing the powers of u 
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Figure 3.7: Heavy quark potential Vqq as a function of quark separation d with 
quartic warp factor, or equivalently, {a,a,/3,/3) = (0.1,4,0.05,4) in the warp factor 
equation for various temperatures. 

in the exponential changes the range of Tc only by a small amount. So effectively 
Tc lies in the range 0.2 < < 0.4. Putting back units, and defining the boundary 
temperaturj^ T as T 



S-^^ — , our analysis reveal: 
1.06 



iwy/ h{uh) 

0.91 



L2 



< Tc < 



L2 



(3.198) 



which is the range of the melting temperatures in these class of theories for heavy 
quarkonium states. Since the temperatures at both ends do not differ very much, 
this tells us that the melting temperature is inversely related to the asymptotic AdS 
radius in large thermal QCD. 



22 



See sec. (3.1) of [57] for details. 
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Figure 3.8: Comparison between lattice QCD results [188 ] and our analysis. The left 
figure is the lattice plot whereas the right figure is our calculation for the potential. 
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Figure 3.9: Slope of linear potential as a function of T for both cubic and quartic 
warp factors. Note that in the figure the slopes have been computed as ^ with a 
range of d from 1.6 to 1.7 in appropriate units. 
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Conclusions 



In this thesis we have proposed the gravity dual for a non conformal finite temperature 
field theory with matter in fundamental representation. The gauge couplings run 
logarithmically in the IR while in the UV they become almost constant and the theory 
approaches conformal fixed point. To our knowledge, the brane construction and the 
dual geometry in section 2.4 is the first attempt to UV complete a Klebanov-Strassler 
type gauge theory with asymptotically conformal field theory. 

Although our construction is rather technical with the gauge group being of the 
form SU {N + M) x SU (N) in the IR, one can perform a cascade of Seiberg dualities 
to obtain the group SU{M) and identify this with strongly coupled QCD. One may 
even interpret that the gauge groups depicted in Fig 2.19 contain strongly coupled 
large N QCD. This is indeed consistent as the coupling Qn+m of SU{N + M) factor 
in SU{N^M) X SU{N + M) or in SU{N^M) x SU{N), always decreases as scale is 
increased. To see this, observe that in the IR, the coupling Qn+m runs logarithmically 
with scale and gets stronger as scale is decreased. At the UV, we can arrange the 
sources in the dual geometry such that Qn+m decreases as the scale grows and runs as 
gM+N ~ O/k/A'^. By demanding —ka^/K^^^ < 0, we see that Qn+m indeed decreases 
as scale is increased. Thus from UV to IR coupling always increases, just like QCD. 

Of course the 't Hooft coupling \n+m — {N -\-M)gN^M for the group SU{N + M) 
is still large in the limit N ^ oo, even if Qn+m has decreased to very small value. 
This allows us to use the classical dual gravity description for the gauge theory which 
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has a coupling that shrinks in the UV, mimicking QCD. For even higher energies, 't 
Hooft couphng will eventually become too small for finite N + M and supergravity 
description will no longer hold. But we can use perturbative methods to analyze 
the field theory in the highest energies. For finite N + M and very high energies, 
the gauge coupling may even vanish giving rise to asymptotically free theory. These 
arguments lead us to conclude that, in principle, the brane configuration we proposed 
can incorporate QCD and in the large N limit, the gravity dual we constructed can 
capture features of QCD. 

Using the dual geometry, we have studied the dynamics of 'quarks' in the gauge 
theory and computed shear viscosity rj and its ratio to entropy rj/s. The key to 
most of our analysis was the calculation of the stress tensor of gauge theory and 
we showed how different UV completions contribute to its expectation value. Using 
the correlation function for the stress tensor, we computed the shear viscosity of the 
medium while the computation of pressure and energy density allowed us to calculate 
the entropy of the system. Using a similar procedure to calculate correlators of stress 
energy tensors, with introducing diagonal perturbations in the background metric, we 
can easily evaluate the bulk viscosity ( of the non-conformal fluid. One can consider 
vector and tensor fields of higher rank to couple to the graviton perturbations in 
the five dimensional effective theory and study how this coupling affects the bulk 
viscosity. The calculation is underway and we hope to report on it in the near future. 

All the calculations we performed regarding properties of the plasma did not 
account the effect of expansion of the medium which is crucial in analyzing fiuid 
dynamics. One possible improvement would be to construct a time dependent dual 
gravity which can describe the expansion of the QGP formed in heavy ion collisions. A 
first attempt would be to consider collisions of open strings ending on D7 branes and 
then compute their back-reactions on the geometry. The gravity waves associated 
with the collisions evolve with time and from the induced boundary metric, one 
can compute the energy momentum tensor of the field theory. Analyzing the time 
dependence of this stress tensor, one can learn about the evolution of the medium 
and subsequently account for the effects it has on the quark dynamics. 
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We have computed the fluxes and the form of the warp factor for our static dual 
geometry, but did not give explicit expressions for the deformation of the internal 
five dimensional metric to all orders in QsNj. However, we have explicitly shown the 
Einstein equations that determine the form of the internal metric and using our ansatz 
in [62] , one can in principle compute coefficients in the expansion of the internal metric 
to all orders in QgNf. Even without a precise knowledge of the internal geometry, we 
were able to extract crucial information about the dual gauge theory and formulated 
how the higher order corrections may enter into our analysis. Most of the calculation 
only relied on the warp factor and with the knowledge of its precise form, we were able 
to calculate thermal mass, drag and diffusion coefficients, rj/s and finally free energy 
of QQ pair. For completeness of the supergravity analysis, we hope to compute the 
exact solution for the internal metric using the ansatz of [62j in our future work. 

In our computation of the heavy quark potential, we classified the most general 
dual gravity that allows linear confinement of quarks at large separation and small 
temperatures. We showed that if a gauge theory has dual gravity description and its 
effective degrees of freedom grows monotonically in the UV, it always shows linear 
confinement at large distances, as long as the dual warp factor satisfies a very simple 
relation given by fl3.192p . Thus fl3.192p can be regarded as a sufficient condition for 
linear confinement of gauge theories with dual gravity. It would be interesting to 
study what are the general brane configurations that allow warp factors which satisfy 
f l3.192p and thus give rise to confining gauge theories. We leave it as a future direction 
to be explored. 

As the dual geometry incorporates features of Seiberg duality cascade, our con- 
struction is ideal for studying phase transitions. The gauge theories we studied have 
description in terms of gauge groups of lower and lower rank. From the gravity dual 
analysis, by cutting the geometry at certain radial location and attaching another 
geometry up to infinity, we can construct gravity description for various phases of a 
gauge theory. Each phase will have different dual geometries attached in the large 
r region while the small r region will be common to all the theories. A flow from 
large r to small r geometry can be interpreted as 'flow' from an effective theory in 
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the UV to another in the IR. Prom UV to IR, the different effective theories will 
describe different phases of the gauge theory and this can allow one to study the 
various phases of dense matter. Thus our construction is not only useful to analyze 
strongly coupled gauge theory, but also has potential for studying phase transitions 
in ultra dense medium and we hope to address this issue in the future. 
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